AO-A055  699 


UNCLASSIFIED 


GRUMMAN  aerospace  CORP  BETHPAGE  N Y F/6  20/9 

AN  automated  procedure  FOR  COMPUTING  THE  THREE-DIMENSIONAL  TRAN— ETC (U) 
FEB  78  W H MAS0N»  C MACKENZIE*  M STERN  F33615-75-C-3073 

AFFDL-TR-77-122-V0L-1  NL 


DOC  FILE  COPY  AOA055899 


AFFDLTR-77122 
VOLUME  I 


AN  AUTOMATED  PROCEDURE  FOR  COMPUTING  THE 
THREE  DIMENSIONAL  TRANSONIC  FLOW  OVER 
WING-BODY  COMBINATIONS,  INCLUDING 
VISCOUS  EFFECTS 


VOLUME  I 

DESCRIPTION  OF  ANALYSIS  METHODS  AND  APPLICATIONS 


GR UMMAN  AEROSPACE  CORPORA TION 
BETHPAGE,  NEW  YORK  11714 


FEBRUARY  1977 

TECHNICAL  REPORT  AFFDL-TR-77-122,  VOLUME  I 
FINAL  REPORl  FOR  PERIOD  1 MAY  1975-1  OCTOBER  1977 


Approved  for  public  release ; distribution  unlimited 


I , 

AIR  FORCE  FLIGHT  DYNAMICS  LABORATORY  N 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES  il 

AIR  FORCE  SYSTEMS  COMMAND  3 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433  j 


78  06  27  101 


V 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever/  Md  the  fact  that  the  government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications , or  other  data,  is  not  to  be 
regarded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture, 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 


Approved  for  public  release;  distribution  unlimited. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


J.  Kenneth  Johnson 
Project  Engineer 


FOR  THE  COMMANDER 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations,  contractual  obligations,  or  notice  on  a specific  document. 


AIR  FORCE/S67tO/l  May  1*78—  ISO 


CLASSI FICATION  OF  THIS  PAGE  flWi«ii  Dmim  Enlmrttl) 

(7m  REPORT  DOCUMENTATION  PAGE 


12.  GOVT  ACCESSION  NO. 


\FFDlfTR-77-122,  Volume  I 


^ ^TOMATED 

Jim;  ^SIGNAL  ' 

■ftONS,  ^CLlf 
_CRIPTION  OF  . 


OCEDURE  FOR 


PUTING  THE 


THE  WRE 
-BODir^ 


SIGNAL  D^NSONIC^OW  OVER  ttlNG-MDYJ^BI 
, pCLlirriNG  3tISC(K^  EF’^ECTsJ^LOTlE  t^S- 
lON  OF  analysis  methods  and  Af^LICATfoffs  t . 


jilLLIAM  £. 
“TU  ANITA 


.LLHAUS* 


aiLLIAM  Ji.fMSON 
^NALDf^CraNZ  IE 
lARK/sT?kNv 


GRUMMAN  AEROSPACE  CORPORATION 
BETHPAGE,  NEW  YORK  11714 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  ^ 

AIR  FORCE  FLIGHT  DYNAMICS  LABORATORY  (J 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 

U MONITORING  AGENCY  NAME  S ADORESSfK  jfom  ConlrolllnA  Of»c»J 


READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

3.  RECIPIENT’S  CATALOG  NUMBER 


EINAL  *EP«IT/ 

|\Y  OCTWW  »77^ 

^ ■■WRMIMnTTftOi  iRBIMBT  imllBI 

TRACT  OR  GRANT  NUMBERC*) 


-75-C-3(^73 


rF33615 


to.  PROCR>frm>EMENT.  PROJECT^ 
AREA^IO^  UNIT  NUMBERS  f * 

PROGRAllft^ENT  622  m 
PROJECTlUm  TASK  147t  FTl 


■nr  number  of  pages 
218 

IS.  security  class,  (of  (him  reparf) 


JNCLASSTFIED 

a.  OECL  ASS^IF'i'cATION  DOWNGRADING'' 

schedule 


16.  distribution  STATEMENT  (of  thia  Raport) 

Approved  for  public  release;  distribution  unlimited. 


DISTRIBUTION  STATEMENT  (of  tha  ahstract  antarad  in  flioc*  20.  ft  diffarant  ftom  Report) 


18.  supplementary  notes 

*AMES  RESEARCH  CENTER,  NASA  and  **  INEOIOIATICR-PMI , INC. 

AEROMECHANICS  LAB,  U.S.  ARMY  AVIATION  R&D  COM>tAND  PALO  ALTO,  CALIFORNIA 
MOFFETT  FIELD,  CALIFORNIA 

19.  key  WORDS  (Corrtfnua  on  r#v#f*«  aide  if  rrecaaamry  and  identify  by  block  number) 

AIRCRAFT  AERODYNAMICS 

JETS,  WAKES  AND  VISCID-INVISCID  FLOW  INTERACTIONS 
SSUBSONIC  AND  TRANSONIC  FLOW 

20.  lU^^RACT  fConllnu.  on  fgvaf  «■  aide  it  n.c.aary  mnd  Identity  by  black  ntimb.r) 

Thl.s  report  describes  a numerical  method  to  predict  the  detailed  pressure 
distribution  and  force  and  moment  results  for  wing-body  combinations  at  tran- 
sonic Mach  numbers  less  than  one.  The  resulting  computer  code  has  been  devel- 
oped with  the  Intent  of  providing  the  user  with  an  easy  to  use  and  reliable 
tool  that  produces  the  most  accurate  possible  engineering  prod  let  ions  .,3  The 
basic  inviscld  prediction  method  is  the  modified  transoni^'Timnl 1 disturbance 
theory  progr.im  developed  by  Ballhaus,  Bailey  and  Frick,  ^n  order  to  provide 


provide 


DD  1473  EDITION  OF  1 NOV  65  IS  OBSOLETE 


, !> 


SECURITY  CLASSIFICATION  OF  T^iS  g^GE  nmtm  Fntered) 

78  06  2i' 


StCUWiTY  CLASSIFICATION  OF  THIS  PAOtflWlMi  Dmim  Enlmn4) 


20.  ABSTRACT  (Continued) 

accurate  surface  pressure  predictions  on  the  wing,  several  additional  features 
of  the  typical  transonic  flow  field  have  been  Incorporated.  These  consist  of 
the  viscous  displacement  effect,  local  strong  viscous  Interaction  at  the  shock 
wave  foot  and  at  the  trailing  edge  (Including  an  aporoxlmate  treatment  of  local 
shallow  separations),  and  finally,  the  Interaction  effect  of  the  fuselage. 


FOREWORO 


This  final  report  was  prepared  by  the  Aerodynajnics  Section  of  the  Grumman 
Aerospace  Corporation,  Dotiipage,  New  York  for  the  P’light  Mechanics  Division, 

Air  Force  Flight  Dyna.’nics  Laboratory,  Wriglit-Patterson  Air  Force  Base,  Ohio. 

The  work  was  performed  imder  Contract  No.  F33f'15-75-C-3073,  wliich  was  initiated 
under  Project  No.  1476,  "Advanced  Wing-Body  Aerodynamic  Analysis  and  Design." 

Mr.  J.  Kenneth  Johnson  (FXiM)  was  the  Project  Monitor  of  this  contract. 

The  report  consists  of  throe  volujiies.  Volume  I,  entitled  "Description  of 
Analysis  Methoiis  and  Applicatio;is , " describes  the  methods  used  to  predict  sur- 
face pressure  distributions  and  aerodynamic  forces  on  three-dimensional  wing- 
body  combinations  at  transonic  speeds,  including  viscous  effects.  Volu-me  I 
also  contains  an  extensive  sot  of  comparisons  between  nuitierical  predictions 
tmd  ex^jerimental  results.  Detailed  instructions  I'equired  to  use  the  prograai 
are  provided  in  Volume  II,  "User's  Manual  and  Cotie  Description. " Volune  III 
was  written  at  Sybucon,  Inc.,  Atlanta,  Georgia  and  contains  a complete  descrip- 
tion of  the  tlieory  and  progrrun  that  computes  tlie  Ihill  tiirce -dimensional  boundary 
layer  over  tlie  wing.  This  work  was  peri'ormetl  by  Sybucon  under  subcontract  to 
Grumman  Aerospace.  Although  this  program  operates  independently  of  the  program 
described  in  Volume  II,  the  input  data  set  reciulred  for  the  ihill  three-dimensiona 
botmdary  layer  computation  is  generated  by  the  code  documented  in  Volune  II. 

Mr.  F.  Berger  was  the  Progr;un  Manager;  Ihu  W.  Mason  and  Mr.  D.  MacKenzie 
served  as  Project  Engineei-s.  'Die  work  was  performed  in  close  cooperation  with 
the  co-authors  from  the  NASA  Ames  Research  Center,  Dr.  W.  F.  Balliiaus  and 
Ms.  J.  Frick.  Additional  contributors  to  the  project  included  G.  Simpers, 

A.  Vachris,  D.  Raila,  P.  Aldala,  M.  Sturm  and  A.  Bunnell  of  Grumman,  and 
Drs.  F.  R.  Bailey  and  V.  Holst  of  NASA  Ames,  Moreover,  contributions  have  been 
made  by  A.  Clien  of  Boeing,  Drs.  R.  Melnik,  B.  Grossman  and  G.  Volpe  of  the 
Grunmfui  Reseai'c}!  Department  ;mii  Gmimian  Consultants  Prof.  A.  Jamieson,  Prof.  J. 
Werner  and  Dr,  £1.  Murman.  As  noted  above,  the  three-dimensional  bouidary  layer 
program  was  written  by  Dr.  J.  Nasii  and  Dr.  K.  Scruggs  of  Sybucon,  Inc. 


TABLE  OF  CONTENTS 


SECTION  PAGE 

I.  ANALYSIS  METHODS  1 

1.  Introduction  1 

2.  Inviscid  Calculation  Method  5 

3.  Viscous  Calculation  Method  14 

4.  Coupling  of  Inviscid  and  Viscous  Programs  38 

a)  General  38 

b)  Shock  Boundary-Layer  Interaction  Effects 47 

c)  Separation  Treatment  62 

5.  Body  Model 64 

6.  Transonic  Wave  Drag  on  Swept  Wings 72 

a)  General 72 

b)  Exact  Expression  for  Lift  and  Drag 73 

c)  First  and  Higher  Order  Expansions  75 

d)  The  Drag  Integral 79 

e)  Wave  Drag 81 

7.  Fully  Three-Dimensional  Boundary  Layer  Calculation  ...  87 

II.  APPLICATIONS 89 

1.  Introduction 89 

2.  Experiences  and  Limitations  of  the  Code  - How  To 

Successfully  Run  the  Code 90 

3.  Examples  of  Configurations  Computed  and  Typical  Agreement 

With  Experiment 95 

a)  ONERA  M6 97 

b)  RAE  "A" 98 

c)  L54H18 99 

d)  TN  D712 100 

e)  L51F07  100 

f)  A55B21 101 

g)  Forward  Swept  Wing 102 

h)  Low  Aspect  Ratio,  Highly  Tapered  102 

1)  Transonic  Transport  103 

j)  Advanced  Fighter  104 

k)  Summary 

I 


V 


ffiECXDlMO  PiOl 


1 


r 

I 


4.  Detailed  Data  Comparisons  - TACT  and  F-8 138 

a)  TACT 138 

b)  F-8 154 

5.  Common  Input/Output  168 

a)  Introduction 168 

b)  Common  I/O  - - A Systems  Viewpoint 168 

c)  Common  I/O  - - The  Aerodynaraicist 169 

d)  Common  I/O NASF 170 

e)  QUICK  - - A Baseline  Code  for  Generating  Body  Definition  170 

f)  Wing  Definition 172 

g)  Feasibility  of  Including  Graphical  Results  172 

h)  Common  Output 172 

i)  Recommendations  For  Common  I/O 173 

6.  Design  Feasibility  177 

a)  Introduction 177 

b)  Review  of  Design  Process 178 

c)  Review  of  2-D  Transonic  Design  Methods  180 

d)  Review  of  3-D  Transonic  Design  Methods  183 

e)  Application  of  the  3-D  Transonic  Program  to  Wing 

Design  Problems 185 

f)  A Practical  Approach  to  an  Intermediate  Design  Code.  . . 188 

g)  Conclusions  on  Design  Feasibility  190 

III.  CONCLUSIONS 197 

IV.  RECOMMENDATIONS  FOR  FUTURE  WORK 199 

REFERENCES 201 


H 


1 


itj 

1 


vi 


11 


PI 


LIST  OF  ILLUSTOATIONS 


FIGURE  PAGE 

1.  Setnispan  Wing I3 

2.  Wing  Planform  Transformation  I3 

3.  Infinite  Yawed  Wing  Notation  25 

U.  Effect  of  Wing  Sweep  on  Roundary  Layer  Development  26 

5.  Comparison  with  the  Nash  Test  Case  C = .I5 27 

6.  Coa.parison  with  the  Nash  Test  Case  C.  = .50 28 

7.  Comparison  with  Adams  3“P  Test  Case  Zero  Lift 29 

8.  Comparison  with  Adams  3~D  Test  Case  C = 0.7I+  Case 3O 

Li 

9.  Comparison  with  the  Compressible  Adams  Test  Case 

^*/c,  e/c 31 


10.  Comparison  with  the  Compressible  Adanis  Test  Case  Cf^^ 32 

11.  Compai’ison  of  Bradshaw  3-P  Program  and  Modified 

Chordwlse  Methods  33 

12.  Comparison  of  Modified  Chordwlse  and  3~D  Calculation  Methods  . . 34 


13-  Flat  Plate  Results  - 0°  and  45'^  Sweep  3^ 

14.  Spanwise  Component  of  Shear  Stress  at  Wall 36 

15.  Spanwise  Component  of  Wall  Shear  Stress  37 

'6.  Some  Viscous- Inviscid  Interactions  in  Transonic  ilow  44 

17.  Examples  of  Viscous  Flow  Solutions  45 

18.  Types  of  Shock  - Turbulent  Boundary  Layer  Interactions  53 

19*  Pressure  Rise  Across  Shocks  54 

20.  Yoshihara's  Viscous  Wedge  Model  55 

21.  Shock  Polar 'io 

22.  Rubber  Wall  Model 57 

23.  Comparison  of  Empirical  & Theoretical  Models  5S 

24.  Implementation SO 

25.  Exajnple  of  Viscous  Weiige  Effect tO 

26.  Shock  Bounitory  Layer  Interaction  Effects  6l 

27.  Schematic  of  Boundary  Conditions  Suppiort  Surface 67 

28.  Simple  Body  Input  for  Infinite  Body  Bi-)undary  Condition  Siurfnce  . 68 

29.  Equivalent  tody  Source  Model  69 

30.  Wing-Body  Interaction  Effects  - Geometry  & Wing  Pressures  ...  70 


vi  1 


LIST  OF  ILLUSTRATIONS  (Contd) 


31.  Wing- Body  Interaction  Effects  - Comparison  of  Pressures  at 

the  Wing- Body  Junction  . . 

32.  ONERA  M6  

33.  RAE  "A"  With  Body  

31*.  L54H18A 

35.  TN  D-712  

36.  L5IPO7  

37.  A55B21  

38.  Forward  Swept  Wing  

39.  Low  Aspect  Ratio  

kO.  Transonic  Transport  

*+!•  Advanced  Fighter 

42.  Three-View  Drawing  of  TACT  Supercritical  Wing  Airplane 

43.  TACT  Wing  Section  Streamwise  Profile  Near  Midspan 

44.  TACT  Wing  Twist .' 

45.  TACT  Wing  Aeroelastic  Twist  Increments  

46.  TACT  Wing  Pressure  Orifice  Locations  

47.  TACT  Wing  Aeroelastic  Effects  . 

48.  TACT  Comparison  Case  M=.9/  a =5.2°  

49.  TACT  Comparison  Case  1^.85,  a=8°  

50.  Local  Trailing  Edge  and  Cove  Separation  

51.  Three-View  Drawing  of  F-8  Supercritical  Wing  Airplane  . 

52.  F-8  Wing  Section  Streamwise  Profile  Neeu-  Midsemispeui 

53.  F-8  Wing  Twist  

54.  Location  of  Wing  Pressure  Orifices  on  F-8  

55*  Tunnel  and  Flight  Test  Spanload  for  NASA  F-8  

56.  F-8  Comparison  Case  M=.9»  a=4°  

57-  F-8  Comparison  Case  M=.99>  ^3° 

58.  Typical  Point  Design  Data  Flow  ...  

59.  Typical  Design  Data  Flow  Requirements  

60.  QUICK  Geometry  

61.  Advanced  Airfoil  Synthesis  

62.  Low  C Tailored  Airfoil  

m 

o 

63.  Inviscid  Drag  Minimization  

64.  Wing  Drag  Minimization  

65.  Wing  Design  Example  


71 

105 

109 

113 

117 

121 

125 

127 

129 

131 

134 

141 
143 

142 

142 

143 

144 

148 

151 

153 

157 

157 

158 

159 

160 
162 
165 

174 

175 

176 

191 

192 

193 

194 
196 


I.IS'l’  OF  TABLES 


I 


I 


I 


TAlUi: 

T EXAMPLE  OiTJFIOUKATIilNS 


A 

AK 


a 


O** 


P 

F(,x'' 

i 

K 

L 

M 

M 

A' 

n 

r 

q 

Re 

s 


S(x> 


0 , 
wake 


s 

U 

e 


LIST  OF  S\W01,S 

liil'luonoo  ot'ielTlelent  jnatrlx 
Aspect  ratio 
Opeed  of  soimd 
Ooeffioient  of  drafc 
Ooofi'ioiont  of  skin  friction 
Ooeffioient  of  lift 
Ooel'i'iojont  of  pressiu*e 
’icadintr  eO^e  tr.-uisltioji  parameter 
iK.-'oal  chord 

Meim  aero<.iyniunic  eliord 
I’otal  wlra^^ 

Boiiy  siinpe  in  Mnnnan  l^ump 

Unit  vector  in  x direction 

I'riuiaonic  similarity  pai'aineter 

Total  lift 

N'ac-h  inrnher 

Frees tretun  NUvcJi  number 

ibcponent  in  Macli  number  relation  in  goveniing  equation 
lYessure 

Velocity  in  exqvuision 
Reynolds  Number 

Area  of  boundiuw  condition  suppoi't  surface 
Local  cross-section  txrea  of  body 

Regions  for  integration 
Arc  length 
Edge  velocity 


\,2,etc. 


^,2, etc. 
W., 


LIST  OF  SYMBOLS  (Continued) 

Siiock  velocities 

I’'lnite  difference  expression  in  transonic  equation 
S^'anwise  edge  velocity  in  yawed  wing  coordinate  system 
Sliock  velocities 

Crossflow  profile  in  boundary  layer 


x.y.z 


Aircraft  coordinate  system 
Yawt'd  wing  coordinate  system 


Subscripts 


ibqperiment 
Leading  etige 
Maximum 

Streamwise  coordinate 
Prediction 

Perpencli  culai'  to  streasiwlse  coordinate 
Sonic  condition 
Trailing  edge 


Creek 


Angle -of-attaek 

Effective  angle-of -attack 

Cieometric  angle-of -attack 

Ratio  of  specific  heats 

Bound.'U'y  layer  thickness 

Pisplacement  thickness 

F.xpansion  parameters 

Balliiaus  finite  difference  notation 

Convergence  criterion 

Potential 

y sptm 


X 


•/»! 


LIST  OF  SYMBOLS  (Continued'' 


Sweep  angle 

Cosine  between  yawed  wing  direction  and  streamline  direction 
Shear  stress 
Relaxation  factor 
I'ensity 

Manentum  thickness 

Momentvim  thicknes.'  in  yawed  spanwise  direction 
Planfom  grid 

Cross  flow  transition  pai'ameter 
Viscosity 


j 


J 


SECTION  I 


ANALYSIS  METJIODS 


1.  INTRODUCTION 

This  report  describes  a numerical  method  of  predicting  the  detailed 
pressure  distribution  and  integrated  forces  and  moments  for  wing-body  combina- 
tions at  transonic  Mach  numbers  less  than  one.  The  resulting  computer  code 
has  been  developed  with  the  intent  of  providing  the  user  with  a convenient, 
accurate  and  reliable  engineering  tool.  The  basic  inviscid  prediction  method 
is  the  modified  transonic  small  disturbance  theory  program  developed  by  BalUiaus, 
Bailey  and  Frick  (Reference  l).  To  provide  accurate  surface  pressure  predic- 
tions on  the  wing,  several  additional  features  of  the  typical  transonic  flow- 
field  have  been  incorporated.  These  consist  of  the  viscous  displacement  effect, 
local  strong  viscous  interaction  at  the  shock  wave  foot  and  at  the  trailing 
edge  (including  an  approximate  treatment  of  local  shallow  separations),  and 
the  interaction  effect  of  the  fuselage. 

The  resulting  code  is  a synthesis  of  the  best  methods  available  at  the 
initiation  of  the  study,  and  represents  contributions  of  a number  of  different 
individuals  througl-iout  the  U.  S.  and  British  Aerospace  community.  The  program 
has  been  applied  to  a large  number  of  cases  for  which  experimental  data  is 
available,  and  to  several  examples  of  extreme  geometric  configurations  for 
which  the  method  will  operate  successfully,  although  no  data  is  available.  The 
experience  gained  diiring  the  effort  has  been  used  to  reduce  the  complexity  of 
the  program  input  to  a simple  practical  form,  allowing  the  user  to  concentrate 
on  the  aerodynamic  aspects  of  the  analysis. ' The  program  shoirld  not  be  expected 
to  produce  results  if  there  are  significant  regions  of  separated  flow,  and  this 
fact  can  be  used  to  determine  the  limits  of  the  program  applicability  for  any 
particular  case.  Body  effects  are  incorporated  into  the  program  by  providing 
an  infinite  rectangular  cross-section  upon  which  the  fuselage  slopes  are  applied. 
Suitable  modification  of  these  slopes  is  made  via  slender  body  theory  to  take 
into  account  the  transfer  of  the  boundary  condition  from  the  fuselage  sirrface 
to  the  rectangular  cross-section  bormdary  condition  support  surface.  In  addi- 
tion, the  program  is  assembled  in  a modular  form  that  allows  for  an  easy  upgrad- 
ing as  various  improved  analysis  methods  become  available. 

Practical  numerical  methods  for  transonic  flows  became  available  in  the 
early  seventies,  as  a result  of  the  pioneering  work  by  Miuman  and  Cole  (Reference  2), 


1 


I 


who  introduotvi  tho  oonot>pt.  of  the  "mixed  differencinti."  I’eluxiitlon  metlicxl  for 
the  ijolutioti  of  ttie  triuiiioiiie  .•iiiinll  ilisturbtuice  equation.  Tlie  stundaixl  two- 
dimenslonai  tuia.Lyaia  iirogr;uu  quickly  hecfune  tlie  Jiuneaon  (Keference  3)  circle 
I'l.'Uie  i'ri>K«‘’Uiu  wtilch  ju-ovides  a numerical  solution  of  the  t\ill  potential  equa- 
tion over  arbitrary  airfoil  sections.  It  was  soon  determined  ttiat  for  airfoils 
wiUi  a slKUlfioant  dcf^i-ee  of  aft  It'adlng,  the  effects  of  the  boundoi-y  layer  liad 
to  be  Incluiled  In  order  to  obtain  gocxi  agreement  wlUi  data,  havltz  (Keference  U) 
producovl  one  of  the  first  prograjns  that  Iteratively  ctimputed  the  inviscid  and 
boimdaiw  layer  flows  until  a eonvergeii  viscous  solution  was  obtained.  I’hese 
methCK.ls  all  contain  some  empiricism  at  the  trailing  edge  to  accoiuit  for  the 
local  strong  interaction  tuul  wake  effects,  neither  of  wliich  were  included  in 
tlie  calculations.  Kecently,  Melnik,  Oliow  and  Meatl  (Reference  hi  develoiied  a 
theory  tliat  explicitly  includes  these  effects.  However,  the  present  tiiree- 
dimensional  metlnxl  employs  ;ui  empirical  model  consistent  witli  tlie  earlier  two- 
tLimensional  metluxis.  I’his  situation  exists.  In  general,  with  three-dimensional 
methods  wiiicli,  as  might  be  exi.x>cted,  are  not  developed  to  t.he  stuiie  level  of 
sophistication  as  tlie  two-dimensional  methoiis  at  any  particular  time. 

Three-dimensional  trtuisonic  flc'W  calculation  metluxis  were  initially  developed 
by  F.  K.  Btiiley  tuid  W.  F.  Ballhaus  at  tlie  NAl’A  Aim's  Research  Center,  using  the 
classical  transonic  small  disturbance  equations.  A three-dimensional  method 
that  solved  the  tlill  potential  equation  was  later  developed  by  Jiunesoii  (Reference 
oi  for  three-dimens i onal  wings.  Other  prognuiis  to  predict  three-dimensional 
t riuisonlc  flows  have  been  developed  by  Hall  (Reference  'i ) iuul  co-workers  at  the 
RAF,  iuid  Cclimidt  ^Reference  8)  fuid  co-workers  at  Hornier  in  liermany.  Fx^perlence 
with  the  classlcsU  transonic  small  disUu'btuice  theory  for  wings  with  mtxlerate 
to  large  sweep  indicated  that  some  refinement  iiad  to  be  incorporated  inti'  the 
small  disturbiuice  theory  In  order  to  properly  simulate  swept  shock  waves.  The 
first  modifications  were  proposed  by  l.omax  (Reference  O',  et  al . , and  it  is 
this  particular  equation  wtiich  has  been  implemented  in  the  baseline  inviscid 
Cixle  employed  in  the  present  method.  This  inviscid  code  is  also  the  only  avail- 
able method  wtiich  can,  at  present,  simulate  relatively  general  tXiselage  config- 
urations routinely.  The  Cixte  mtikes  use  of  the  embedded  grid  scheme  developed 
I'y  Hoppe  (Reference  10)  in  order  to  reduce  the  execution  time  and  storage  require- 
ments, wtiiie  retaining  maximum  accuracy  on  the  surface. 
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Viscous  effects  are  investigated  and  accounted  for  in  the  present  method 
through  two  different  approaches.  The  full  viscous-inviscid  iteration  is 
carried  out  assuming  a boundary  layer  of  the  infinite  swept  wing  type  at  each 
span  station,  augmented  by  local  treatments  of  the  strong  interaction  regions 
at  the  shock  foot  and  trailing  edge,  and  regions  with  shallow  separations.  In 
addition  to  this  calcxilation,  an  entirely  new  program  has  been  written  by  Nash 
and  Scruggs  (Reference  11 ) to  compute  the  fully  three-dimensional  laminar  and/ 
or  turbulent  boundary  layer  on  finite  wings.  The  Nash -Scruggs  program  is  not 
fully  coupled  with  the  inviscid  calculation,  however,  the  inviscid/strip  vis- 
cous program  will  automatically  produce  the  input  data  set  for  the  3-D  BL  pro- 
gram. The  resulting  3-D  boundary  layer  prediction  is  then  tabiilated  in  a form 
that  allows  for  the  generation  of  a data  deck  that  could  be  used  to  make  a 
fully  3-D  correction  to  the  actual  airfoil  ordinates.  The  only  other  reported 
effort  to  automatically  combine  viscous  and  inviscid  flow  calculation  at  tran- 
sonic speeds  is  the  work  by  Schmidt  and  Hedman  (Reference  12),  which  presents 
calculations  for  conventional  airfoils  only.  Kordulla  (Reference  13)  has  also 
investigated  the  coupling  of  a 3-D  boundary  layer  method  with  the  Bailey-BalUiaus 
Program. 

The  present  method  is  described  in  a three  volume  report.  Volume  I de- 
scribes the  theoretical  foundation  of  the  various  elements  of  the  method  and 
how  they  are  combined  into  a single  computer  program.  No  attempt  is  made  to 
repeat  the  excellent  and  detailed  theoretical  descriptions  contained  in  the 
references.  Instead,  the  volume  concentrates  on  providing  the  user  with  a 
practical  overview  of  the  methods  and  a working  knowledge  of  those  aspects  of 
the  methods  which  the  user  can  control.  Such  aspects  are  seldom  detailed  in  the 
theoretical  papers.  This  volume  also  contains  a large  number  of  examples  of 
the  application  of  the  program  to  a complete  range  of  aircraft  configurations. 
These  include  correlations  on  the  F-8  and  TACT  aircraft.  The  volimie  concludes 
with  some  observations  on  the  general  integration  of  computational  aerodynamic 
tools  into  a uniform  system  with  common  input  and  output  sets,  and  an  examina- 
tion of  the  potential  applications  to  the  desig;i  problem  using  the  present 
analysis  program  as  a baseline  method.  Volvime  II  contains  the  detailed  instiaic- 
tions  required  for  program  operation,  including  a review  of  the  output  results 
and  the  associated  notation/ definition. The  first  part  also  contains  a sample 
case  which  should  provide  a complete  illustration  of  the  use  of  the  program. 

The  second  part  of  this  volume  contains  a description  of  the  actual  computer 
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program  in  sufficient  detedl  so  that  a user  can  learn  the  code  well  enough  to 
make  modifications.  Hie  three-dimensxonal  boundary  layer  method  is  treated 
as  an  independent  program  and  is  entirely  described  in  Volume  III.  That  volume 
contains  both  the  theoretical  description  and  the  user's  manual  and  code 
description. 
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2.  INVISCID  CALCULATION  METHOD 


The  inviscid  transonic  flowfield  prediction  is  obtained  using  the  current 
version  of  the  small  disturbance  theory  program  developed  at  NASA  Ames,  This 
program  has  been  under  constant  development  for  a nvmber  of  years,  as  reported 
in  a series  of  publications.  The  work  has  now  reached  a relatively  mature 
stage  of  development  and  forms  the  basis  for  the  present  method,  which  is  in- 
tended to  be  used  as  a reliable  engineering  tool.  The  program  was  developed 
on  the  CDC  7600  computer  and  the  structure  of  the  code  reflects  the  requirements 
of  that  particular  machine.  It  can,  however,  be  readily  converted  to  operate 
on  IBM  machines  that  have  an  operating  system  with  the  VS  feature.  The  program 
is  big  and  long  running  by  present  engineering  standards  (although  these  are 
changing  rapidly),  so  that  users  will  require  access  to  moderately  advanced  com- 
puters. Although  the  detailed  machine  requirements  are  discussed  in  Vol\3me  II, 
the  user  should  study  the  method  with  these  considerations  in  mind. 

The  present  section  is  intended  to  provide  the  salient  features  of  the 
method,  along  with  a candid  description  of  the  assumptions  and  approximations 
that  have  been  adopted.  Recall  that  transonic  small  disturbance  theory  is 
valid  in  a double  limit  as  M - 1,  6 - 0,  where  M is  a freestream  Mach  number 
and  6 is  a measure  of  the  disturbance  size  such  as  the  thickness  to  chord  ratio 
of  the  basic  airfoil.  Ballhaus  (Reference  l4)  has  recently  reviewed  the  3-D 
small  disturbance  theory  in  the  general  context  of  contemporary  transonic  aero- 
dynamics and  the  review  is  recommended  as  an  excellent  introduction  to  the 
intricacies  of  the  3*D  small  disturbance  methods, 

A modified  form  of  the  small  disturbance  equation  (MSD)  is  included  as 
the  basic  option  in  the  code.  Its  use  is  recommended  for  the  treatment  of  flows 
about  configurations  with  swept  wings  (References  9 !)•  The  MSD  equation, 

written  in  conservation  form,  is 

M^""  0/  + i (y-3)  0y^j^ 

+ j^0y  - (y-1)  M„^  ^x  '^y  ] y ^ ^ ° 

where  0 is  the  disturbance  potential,  y ^s  the  ratio  of  specific  heats,  and 
is  the  freestream  Mach  nxmiber. 
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Figure  1 indicates  the  orientation  of  the  vehicle  with  respect  to  the 
coordinate  system.  Note  that  the  program  computes  the  flowfield  in  the  wind 
axis  system.  The  parameter  "n"  that  appears  as  an  exponent  of  the  Mach  number 
in  the  coefficient  of  the  traditional  non-linear  transonic  term  can  be  used 
to  improve  the  small  disturbance  approximation  to  the  correct  shock  jump  con- 
ditions. The  value  of  n = 1.75  has  been  found  very  satisfactory  for  the  cal- 
culations performed  thus  far  (n  = 2 being  the  classical  value). 

The  solution  of  equation  (l)  subject  to  the  boundary  conditions  specified 
below  is  obtained  by  approximating  the  terms  by  finite  differences  and  solving 
the  resulting  set  of  algebraic  equations  by  successive  line  overrelaxation 
(SLOR)  in  a computational  mesh  space. 

The  wing  boundary  conditions  are  applied  at  the  wing  reference  plane, 
which  forms  a slit  between  mesh  lines  in  the  computational  space.  The  linear- 
ized boundary  condition  is  modified  in  order  to  include  a Krupp-type  scaling 
(Reference  15),  which  has  been  found  necessary  to  provide  good  results; 


(x,y,+0)  = i i f (x,y,+0)  — o 


The  computational  mesh  is  truncated  at  a distance  several  wing  spans  from  the 
wing,  and  the  potential  is  specified  by  the  asymptotic  farfield  formula  given 
by  Klunker  (Reference  l6).  Only  that  part  of  the  potential  due  to  the  aircraft 
lift  is  specified,  the  rest  of  the  contributions  being  of  higher  order.  Recall 
that  in  three-dimensions,  the  potential  decays  as  (l/r),  while  in  2-D,  the 
potential  varies  as  log  (r),  where  r is  the  origin  of  the  win^'airfoil 
system.  This  indicates  that  the  asymptotic  matching  between  the  numerical 
and  analytical  value  of  the  potential  is  not  as  critical  in  3'D  flows  as  for 
2-D  transonic  flows.  The  trailing  vortex  sheet  extends  downstream  in  the  plane 
of  the  wing.  Althovigh  initially  the  Trefftz  plane  problem  was  sol/ed  explicitly 
at  the  downstream  boundary,  it  was  found  nvmierically  simpler  to  require  = 0 
at  the  downstream  boundary.  This  procedure  is  automatically  equivalent  to  the 
traditional  Trefftz  plane  boundary  condition. 

Once  the  computation  is  completed,  the  value  of  the  pressure  coefficient 
is  determined  fran  the  standard  linearized  relationship. 

Cp  = -2  0^  (3) 
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The  numerical  solution  scheme  devised  by  Murman  and  Cole  (Reference  2) 
is  an  elegantly  simple  method  that  manages  to  capture  the  physical  character- 
istics of  a flowfield  which  contains  regions  of  both  supersonic  and  subsonic 
flow.  When  the  local  flow  is  subsonic  and  hence,  governed  by  an  elliptic  par- 
tial differential  equation  locally,  central  difference  finite  difference  approx- 
imations are  employed.  When  the  flow  is  locally  supersonic,  the  locally  hyper- 
bolic natiire  of  the  governing  equation  is  enforced  by  the  use  of  an  upwind 
finite  difference  approximation  for  the  0^  term  in  the  classical  transonic 
small  dist\irbance  theory.  Equation  (l)  differs  from  the  classical  small  dis- 
turbance equation  (CSD)  by  the  addition  of  the  two  underlined  terms,  which  come 
from  the  full -potential  equation. 

These  two  terms  introduce  seme  additional  complexity  into  the  finite - 
difference  algorithm  for  supersonic  regions.  For  the  CSD  equation,  the  x- 
coordinate  is  the  axis  of  the  characteristic  cone.  Thus,  upwind  differencing 
the  X derivatives  and  central  differencing  the  y and  z derivatives  leads  to  a 
numerical  domain  of  dependence  that  always  includes  the  mathematical  domain  of 
dependence;  consequently,  a necessary  condition  for  stability  is  maintained. 
However,  for  the  M3D  equation,  the  characteristic  cone  aocis  lies  in  a direction 
that  corresponds  to  the  local  flow  direction.  To  maintain  stability,  domains 
of  dependence  requirements  are  satisfied  by  proper  combinations  of  central  and 
upwind  differences. 

The  proper  combination  of  central  and  upwind  differences  can  be  determined 
by  expressing  the  principal  part  of  the  governing  equation  in  am  intrinsic 
coordinate  system  (Reference  6). 


(a^-q^)  0 + a^  0 + a^  0 

' ^ ss  nn  zz 


(M 


where  q and  a are  particle  and  so\ind  speeds,  and  s and  n are  the  local  stream 

and  stream-normal  directions  in  the  x-y  plane.  The  0 , 0 and  t contribu- 

XX  xy  yy 

tions  to  0 are  upwind  differenced  and  the  remaining  terms  are  central  differ- 
s s 

enced.  The  terms  in  equation  (k)  can  be  written,  approximately,  in  the  form 
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(a'^-q^Va^^  = 1-M"„  - (6+l)  0^ 


0 =0  +2 
'^ss  ^xx 


i^/a^  = 1 - (6-1)  0 

' m ' ' 1^1  “ %i 


0 = -2 
nn 


xy  '^yy 


The  mixed  central -upwind  difference  procedure  is  implemented  in  the 
following  way.  To  begin  with,  all  terms  in  the  equation  are  centrally  differ- 
enced. Then  a test  is  performed  to  determine  if  the  flow  is  supersonic.  The 

2 2 2 

supersonic  condition  is  (approximately)  from  equation  (5),  (a  -q  )/a  ^ < 0. 

If  the  flow  is  supersonic,  the  following  terms  are  added  to  the  previously 
computed  central  difference  approximation  to  equation  (l). 


Ax  6 (V  6 ) + 2 

^ "j,k 


(v  + (1-v)  6,  ) 


] pj,k  (6] 


where  6 and  6 are  central  difference  operators  that  operate  on  0,  , to 
XX  xy  ^ j ,K 

give  0 and  0 ; 6x  is  an  operator  such  that  6xA.  . = Ax  (A.  . - A 

XX  xy  i >0  i > J ^ >3 ) 

— fc  ^ 

6 and  6 are  similar  forward  and  backward  difference  operators  in  the  y 

y y 2 n 

direction:  V is  a finite -difference  approximation  to  1-M‘'„,  - (v*^l)M^  0 ; 

^ ^ <x>  x 

and  V = (^)  for  0^.  ( $ ) zero.  Note  that  the  first  term  in  equation  (o)  is 

treated  conservatively;  the  other  term,  the  one  involving  0 , is  not.  The 

xy 

differencing  in  equation  (6)  is  ccanplicated  by  the  use  of  tl)e  wing  planform 
transformation  but  can  be  implemented  in  a direct  manner. 

The  resulting  equations  are  solved  by  a line  overrelaxation  scheme  that 
marches  throu^  the  flowfield  along  vertical  lines,  starting  at  the  upstream 
boundary  and  proceeding  to  each  downstream  mesh  pliine  by  working  outward  from 
the  center  line  to  the  outboard  limit  of  the  mesh.  The  general  mathematical 
theory  of  relaxation  applied  to  the  solution  of  partial  differential  equations 
of  mixed  type  in  transonic  flows  has  been  discussed  in  some  detail  by  Jjuneson 
in  References  17  ajid  iB. 

Although  the  mathematical  problem  has  been  IMLly  described  above,  several 
additional  refinements  are  introduced  in  order  to  obtain  a rapid,  practical 
and  reliable  numerical  solution  scheme.  These  inclvxde  the  choice  of  meslies. 


a nun\erical  transformation  to  a coordinate  system  more  closely  aligned  with 
the  planfonn,  a proper  choice  of  the  over relaxation  factor  lu,  and  a leading 
edge  boxindary  condition  mollification  via  Riegels'  Rule>^  in  order  to  improve 
the  acc\iracy  of  the  solution  in  regions  where  the  local  coordinate  system  is 
not  aligned  with  the  leading  edge.  Hie  next  several  paragraphs  describe  tlie 
purpose,  implementation  and  use  of  these  refinements  in  sufficient  detail  for 
the  user  to  gain  enough  insiglit  and  understanding  to  work  with  the  program. 

The  choice  of  computational  mesh  can  have  an  Important  effect  on  the 
accuracy  and  cost  of  the  inviscid  solution.  The  program  internally  generates 
a mesh  system  which  has  been  found  to  provide  a good  solution  at  reasonable 
computing  time.  However,  the  user  can  override  the  internal  mesh  with  an  in- 
put mesli  of  his  own  dioice.  The  program  requires  two  basic  mesii  systems  in 
order  to  obtain  the  solution  using  the  Boppe  (Reference  10)  embedded  mesh 
scheme.  An  exterior  coarse  Cai'tesian  mesh  is  used  over  most  of  the  computa- 
tional domain  and  an  embedded  fine  mesh  is  employed  in  a region  very  near  tlie 
wing  surface,  wtiere  the  flow  gradients  are  large  and  the  solution  is  of  the 
most  interest.  The  embedded  mesh  is  placed  in  a coordinate  system  locally 
aligned  with  the  planform. 

In  general,  the  e.xterior  mesh  should  extend  about  10  root  chords  above, 
below,  fore  and  aft  of  the  wing.  The  present  program  allows  for  a maximum  of 
20  vertical,  20  spanwise  and  30  streamwise  stations.  At  each  span  station, 
at  least  four  streamwise  mesh  stations  should  occur  on  the  wing  in  t)ic  cinide 
grid.  This  requirement  usually  leads  to  tiie  use  of  the  maximvmi  30  streaniwise 
mesh  stations.  The  typical  span  mesh  contains  about  two-thirds  of  the  mesh 
within  the  wing  tip,  with  the  tip  centered  between  mesh  lines  and  tiie  mesh 
widening  rapidly  as  the  distance  outboard  of  the  wing  tip  increases.  Tlie 
vertical  mesh  is  normally  placed  symmetrically  about  the  wing  plane  with  the 
first  mesh  points  about  10^.  of  the  reference  chord  above  and  below  the  wing. 
All  meshes  should  vary  smoothly.  Specifically,  the  second  difference  of  the 
mesh  spacings  should  vary  smoothly. 

* Althou^  the  relationship  of  Riegels'  Rule  to  subsonic  small  distrubtuice 
theory  is  well  understood  (see  Reference  19 )>  its  use  in  transonic  small 
disturbance  theory  is  piarely  ad  hoc  at  present.  Nevertheless,  it  has  been 
found  to  aid  in  obtaining  the  solution  for  some  cases.  Numerical  studies 
have  shown  that  Riegels'  Rule  virtually  eliminates  any  mesli  dependence  on 
the  solution. 
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'rhe  fine  mesii  is  much  denser  and  is  confined  to  a region  close  to  the 
wing.  Typically,  the  mesti  extends  1/3  ~ 1/2  chord  length  upstream  and  down- 
stream of  the  wing,  about  10)(.  of  the  semispan  past  the  tip  and  about  two  chords 
above  and  below  the  wing.  The  vertical  mesh  is  symmetrically  placed  about  tlie 
wing  with  the  first  mesh  points  about  1 ~ 2%  above  and  below  the  surface.  The 
typical  values  for  the  inner  mesli  tjre  20  vertical,  30  spanwise  and  lO  stream- 
wise  stations.  The  smooth  variation  of  the  fine  mesh  is  especially  crucial 
to  ttie  efficient  operation  of  the  program.  Uneven  meshes  not  only  degrade 
the  accuracy,  but  also  have  been  found  to  cause  numerical  iteration  cchvergence 
problems.  These  problems  can  cause  tlie  solution  to  become  much  more  expensive 
due  to  a reduced  convergence  rate,  and  in  some  cases,  even  cause  tlie  solution 
to  diverge.  It  is  well  known  that  solutions  to  the  small  disturbance  equations 
are  dependent  on  the  mesh  distribution  near  the  leading  edge.  In  addition, 
the  exact  solution  of  the  equation  contains  a singularity  at  the  leading  edge 
and  it  has  been  sliown  that  the  exact  solution  of  the  small  disturbance  tlieory 
is  not  a good  approximation  to  the  full  potential  equation  for  a distance  of 
several  nose  radii  downstream  of  the  leading  edge  (Reference  20 ).  Hence,  the 
mesh  is  usually  distributed  about  the  leading  edge  in  a mariner  that  leads  to 
an  incorrect  solution  of  the  small  disturbance  equation,  but  a better  approxi- 
mation to  the  full  potential  equation.  This  is  typical  of  small  disturbance 
solution  methods. 

The  proper  choice  of  coordinate  system  is  important  in  obtaining  accurate 
numerical  solutions.  In  three-dimensional  flows  about  wing-body  combinations, 
the  choice  of  a single  coordinate  system  that  is  appropriate  everywliere  on  the 
vehicle  is  extremely  difficult.  In  fact,  the  study  of  coordinate  systems  and 
the  restilting  numerical  meshes  is  presently  the  pacing  item  in  the  development 
of  more  advanced  inviscid  transonic  programs.  In  fact,  the  latest  Jasieson  ' 

Caughey  (Reference  21)  program  is  significant  in  that  the  numerical  procedure 
is  independent  of  the  specified  coordinate  system.  The  success  of  small  disturbance 
theory  in  large  part  is  due  to  tlie  geometric  simplification  which  arises  from 
the  transfer  of  the  boundary  condition  to  a planar  coordinate  surface.  The 
Boppe  embedded  mesh  scheme  is  also  important  to  the  success  of  the  numerical 
method  by  allowing  a dense  mesh  to  be  placed  in  the  regions  wiiere  it  is  actually 
required.  Moreover,  a "planform"  coordinate  system  has  been  employed  in  the 
small  disturbance  theory  by  introducing  a change  of  variables,  wliich  allows 
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for  the  spanwise  vai'iation  of  the  planfomi  stiapc.  The  appropriate  mapping 
is, 

5 (K,y)  = - ■H-O  <y' 

(y^  - x,^o  (i’) 

n (y^  = y 

wtiere  ^ (y)  = Xj  (y')  jukI  (y^  = Xpj,  (y''  over  the  main  part  of  tlie 

planform,  allowing  for  the  same  number  of  mosli  points  to  cover  the  airfoil 
at  eacli  sptui  station,  yleMing  a very  efficient  mesii  distribution.  A imiform 
spanwise  distribution  of  tlie  streajnwise  mesli  about  the  leading  e<.ige  is  particu- 
larly important  to  eliminate  wiggles.  Tiiere  arc,  iiowever,  some  regii'ins  (.m  tiie 
wing  vrtiere  the  ? = 0 and  ? = 1 lines  do  not  correspond  to  tlie  leading  luid 
trailing  edges.  Hgui-e  2 shows  tlie  mapping  in  gener.'il,  alivig  with  the  six-cial 
areas  vrtiere  the  ' constiuit  lines  deviate  from  lines  of  constant  percent 
cliord.  The  main  deviatitTi  occurs  due  to  the  requirement  that  the  mapping 
intersect  tlie  Idiselage  side  pei’penilicultu'ly.  Tliis  requirement  leads  to  a 
dcpai'ture  of  the  ' 0 line  frcmi  the  leading  etlge  at  the  vring  I’oot  juncture. 

For  gloved  cciifigurations , this  bending  usually  results  in  a loss  of  resolu- 
tion of  the  solution  in  tlie  vicinity  of  the  laxly.  Althougli  tlie  mapping  ;uid 
meshes  are  generated  internally,  the  user  caji  input  his  own  mapping.  This 
would  allow  for  a s^vcial  ccincentraticn  of  the  mesh  near  the  lx\l.v  mid  glove 
in  cases  vrtiere  the  user  wishes  to  concentrate  on  the  glove  results,  leading 
edge  sweep  angles  of  up  to  70^  have  been  treated  without  difficulty  diu-ing  the 
present  study. 

Another  featui'o  of  the  pi'esent  method  is  the  use  of  mi  artificial  time 
step  in  order  to  ensiu’o  a stable  numerioal  scheme.  This  approach  was  first 
adopted  by  Jameson  (Reference  o^  for  trmisonic  flow  solutions.  The  program 
assigns  a numerical  value  for  tlie  degree  of  artificial  time  dmtipliig  to  be 
used.  The  automatically  specified  value  cmi  be  overridden  by  tlie  user  if 
desired.  In  general,  tlie  larger  the  vmluc  of  Uie  dmiiping  factor,  t-he  more 
stable  the  Iteration  procedure,  at  the  expense  of  a reduced  convergence  rate. 
ITius,  the  smallest  value  of  diimping  which  allows  for  stable  cvvivergence  is 
the  proper  choice.  Normally,  tlie  automatically  specified  value  is  a good  c<.w- 
promise  and  it  will  not  have  to  be  myusted  e.xcept  in  hi  (.lily  lumsual  cases. 


Ihe  same  advice  holds  for  the  value  of  the  overrelaxation  factor  u)  (always 
<2.)  which  can  be  timed  for  particular  cases,  but  in  general,  will  not 
require  adjustment.  We  note  that  independent  values  of  u)  and  artificial  time 
damping  (Oxt  damping)  can  be  specified  for  each  mesh. 


3.  VISCaiS  CAI.CIILATION  METHOD 


Viscous-inviscid  interaction  effects  are  predicted  in  the  present  method 
by  iterating  between  the  inviscid  flow  solution  obtained  in  the  majiner  de- 
scribed in  the  previous  section  and  a boundary  layer  prediction  obtained  by 
assuming  at  each  spanwise  station  that  the  boundary  layer  is  of  the  infinite 
swept  wing  type.  The  resulting  displacement  thickness  is  then  added  to  the 
airfoil  ordinates  in  order  to  incorporate  the  weak  interaction  effects.  The 
strong  interaction  effects,  wliich  require  special  treatment,  at  shock  waves, 
the  trailing  edge  and  for  locally  siiallow  flow  separations  are  described  in 
the  next  section.  In  this  section,  we  describe  the  approach  used  to  obtain 
the  infinite  swept  wing  boundai'y  layer  prediction.  A great  deal  of  tliree- 
dimensional  boundary  layer  methodology  development  work  occurred  concurrently 
with  the  development  of  the  present  program  (including  the  contribution  of 
Volume  III  of  this  report).  The  method  adopted  for  the  simulation  of  the 
three-dimensional  viscous  effects  consisted  of  a simple  extension  to  tlie 
method  used  successiAilly  at  Giaunman  in  two-^limensions.  Its  availability 
allowed  for  immediate  coupling  with  the  inviscid  program.  Consideration  of 
the  replacement  of  this  method  with  other  emerging  codes,  such  as  the  progrrun 
contained  in  Volvuiie  III  or  the  latest  program  by  Professor  Bradshaw  (Reference 
22),  awaits  an  extensive  assessment  of  all  the  new  codes  as  well  as  a specific 
need  for  such  a revision.  This  need  was  not  demonstrated  dui’ing  the  extensive 
code  operation  that  was  conducted  during  tliis  effort. 

The  infinite  yawed  wing  coordinate  system  is  shown  in  Figure  3-  This 
coordinate  system  is  rotated  by  the  sweep  angle  A from  the  coordinate  system 
of  the  inviscid  program.  These  rotated  coordinates  are  denoted  by  x,  y,  s. 
Figure  4 shows  results  computed  in  the  inl'inite  yawed  wing  coordinated  system 
and  presented  in  the  strearawise  direction.  The  need  for  an  infinite  yawed 
wing  type  boundai-y  layer  calculation  as  opposed  to  a 2-P  streamwise  strip  is 
brou^it  out  clearly.  The  computation  of  the  boundary  layer  in  the  x,  'z  coor- 
dinate system  can  be  cemputed  as  2-P  only  if  the  flow  is  laminar  and  incom- 
pressible. This  property  is  called  the  independence  principle  and  it  does  not 
hold  for  turbulent  boundary  layers  or  compressible  laminar  boundai'y  layers. 
However,  it  has  been  found  (Reference  23)  that  for  typical  transonic  boundary 
layers,  the  laminar  boundary  layer  can  be  computed  assuming  the  independence 
principle  holds  without  making  appreciable  errors.  This  approach  has  been 


usi'J  by  Nasili  and  Scru^KB  (Hcl'crcnce  pli ) for  the  computation  of  tlio  transonic 
boundtvry  layer  development  on  the  K-H  supercritical  vrinp.  H Is  also  .jiisti - 
fled  In  view  of  the  point  transition  moilcl,  tlie  insensitivity  of  tiie  turl)- 
ulent  touniiory  layer  development  to  the  details  of  the  lamlnar-transltlon 
point  callbulation,  and  tiie  small  region  of  laminar  flow  on  swept  wings.  TJius, 
the  lamirfar  boundary  layer  is  cemputed  In  the  present  program  by  the  method 
of  Thwaitec  general!  :’.cd  for  compressible  i'low  using  tlie  inviscld  flow  normal 
to  the  leading  edge.  The  turbulent  boundary  layer  calculation  does  require 
modification  in  order  to  simulate  the  three-dimensional  effects.  We  will  use 
the  modifications  suggested  by  Nash  and  Tseng  (Reference  r‘3)  for  incompress- 
ible flow  over  an  infinite  yawed  wing.  Their  analysis  method  is  similar  to 
the  IJradshaw  (l?cfcrcnce  ffi)  imalysis  in  vise  at  Grumman  in  ttiat  ;in  equation 
for  tlie  shear  stress  is  solved  along  with  the  usual  momentum  ;ind  continuity 
equations.  The  extra  equation  in  the  Bradshaw  method,  as  compared  with  alge- 
braic eddy  viscosity  metiiods,  provides  the  flexibility  that  allows  tiic  P-O 
program  to  simulate  the  3-1'  effects. 

The  movUfied  chordw:ise  method  was  devised  based  on  the  following  features 
that  Nash  observed  as  a result  of  some  numerical  experiments  for  flow  over  nji 
infinite  yawed  wing: 

o A ?-n  cjilculation  performed  in  the  chordwi se  (x)  direction  tivkes  the 
effect  of  pressure  gradients  on  the  bounviary  layer  development  into 
aceexint  properly. 

o A P-P  calculation  performed  in  the  strcoml ine  direction  would  account 
for  the  shear  stress  direction  and  magnitude  "properly." 

By  assuming  t)iat  the  shear  acts  in  tlic  local  Inviscid  surface  streiunlinc 
direction  througinujt  Uic  sliear  layer,  Nash  has  siiown  Uiat  tlie  equation  govern- 
ing the  chordwilsc  component  of  shear  stress  is  ncai'iy  Identicfil  to  the  usual 
two-dimensional  sJicar  stress  equation.  If  we  write  the  chordwisc  cemponent  of 
the  shear  stress  as 


an<i  assume  that  X varies  slowly  in  the  x-dlroction  compared  to  the  shear  stress 
change,  then  the  equation  for  the  shear  stress  is  identical  to  tlie  2-n  form, 


1^^ 


with  sliglitly  modified  foms  of  the  Bradsliaw  tnnpirical  IVmctlonB.  'fhese 
approximations  are  the  basis  of  the  modified  chordwise  methixi.  'liie  incorpor- 
ation of  tiiis  sclieme  into  the  C-D  Bradsiuiw  prograjn  is  straightforward.  Equa- 
tion (8)  leads  to  the  requirement  that,  at  each  point  on  the  surface,  the  shear 
stress  vector  acts  in  the  sasie  direction  all  the  way  througli  tlie  boundary  layer. 
This  approximation  is  acceptable  for  most  cases  except  in  t)ic  neigjiborhood  of 
separation  lines. 

The  possible  compressibility  effects  were  not  considered  by  Nash.  In  tlie 
general  S-h  sciieme  used  for  conpressible  flow  over  finite  wings,  Nash  tmd 
Scruggs  asBvmie  that  the  shear  stress  equation  is  inaltered  by  conpressibillty 
effects.  Tills  assumption  is  used  in  the  equations  solved  by  Nash  and  Scruggs 
in  the  full  3~b  progivuti  described  in  Volume  III.  In  the  Bradshaw  analysis 
employed  in  the  C-It  calculation,  an  additional  term  arises  in  the  turbulent 
kinetic  energy  equation  due  to  compressible  effects.  This  teim  introduces  yet 
anotlier  empirical  constant  into  the  scheme.  A numerical  experiment  was  made 
in  wliicl.  this  constiuit  was  sot  to  aero  (so  that  the  extra  compressibility  term 
was  omitted  from  the  calculation).  I’or  trmisonlc  flow  speeds,  it  was  fixuid 
that  this  term  did  not  affect  the  results.  Thus,  altliougli  the  term  is  retained 
in  tile  present  calculation  mctliod  witiiout  alteration,  tlie  basic  difference  in 
the  compressible  turbuleiit  shear  stress  equations  used  by  Bradsliaw  !uid  Nasli 
does  not  appear  in  tlie  predicted  results  at  trmisonic  flow  conditions.  Tlie 
only  modification  to  tlie  methoil  to  take  into  accovuit  3-1^  compressible  effeet.s 
luviscs  in  the  use  of  the  Crocco  Integral  in  ordei'  to  obt.ain  tlie  ilensity  varia- 
tion. Tlie  density  appears  exiilicltly  in  the  equations  only  tlirougii  the  combin- 
ation (dp/d^^)/p.  I'lnploying  the  Crocco  Integral,  this  combination  csui  be 
written  as, 
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lx,  in 


Additional  factor  due  to  3-1'  flow  effects. 


for  infinit.e  yawed  wings,  we  find  no  chiuige  in  the  x-diroction  exiiression 
due  to  3-1'  effects.  The  additional  factor  in  equation  roquii-e..  a knowledgi' 
of  the  spanwlse  flowfield,  wliich  the  modified  chordwise  metlnxi  does  not  jiro- 
vide.  niyslcally,  tlie  effect  of  neglecting  tlie  extra  pieee  in  the  Croeoi-' 
Integral  is  to  kee]'  tlie  wall  from  reaeliing  the  liill  adialiatie  wall  temperature. 


! 
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Conpariaon  with  a compressible  solution  shews  that  the  modified  cliordwlse 
technique  predicts  separation  (r^  = 0)  before  the  3-D  solution.  Recall  that 
as  the  wall  becomes  warner,  the  separation  point  moves  foi-ward.  Hence,  it 
appears  that  if  the  3-D  contribution  shown  in  equation  (9)  v/ere  included  in 
the  modified  chordwlse  solution  agreement  with  tlie  3-D  solution  would  not  be 
improved.  Thus,  it  appears  that  neglecting  the  3-D  contribution  to  Uie  Orocco 
Integral  is  consistent  with  the  otlier  approximations  in  the  method. 

The  remaining  modification  of  the  2-D  program  to  simulate  3-D  effects 
requires  the  introduction  of  an  effective  Reynolds  number  which  is  found  by 
requiring  the  component  of  tlie  shear  in  the  x-direction  to  bear  the  same  rela- 
tionship to  an  effective  Reynolds  nmber  that  the  actual  shear  stress  has  to 
the  specified  Reynolds  number.  To  find  this  effective  Reynolds  number,  the 
Karman-Schoenherr  formviln  was  used  so  that  the  Reynolds  number  based  on  mcmen- 
tura  thickness  (wliich  is  specified  at  transition)  was  used.  Compressibility 
effects  were  incorporated  altliougii  Mach  number  effects  on  the  value  of  tlie 
effective  Reynolds  number  were  found  to  be  small.  Typically,  for  a sweep 
angle,  tlie  effective  Reynolds  number  is  about  h0%  of  the  actual  Reynolds 
number . 

Finally,  as  prevlimsly  noted,  it  la  importujit  to  realize  that  the  solu- 
tion in  the  spanwlse  direction  ctuinot  be  computed  using  Uie  modified  chord- 
wise  technique,  and  some  other  method  must  be  used  to  compute  this  infomiatlim. 
An  approximation  for  the  sivinwise  component  of  shear  stress  can  be  obtained  by 
making  use  of  the  observed  features  of  the  viscous  flow  over  an  Infinite  yawed 
wing  and  the  chordwise  solution.  The  first  important  obseiaration  la  that  the 
spanwlse  flow  is  basically  "flat  plate  like"  in  natui-e  since  there  is  no  span- 
wise  pressure  gradient  (Reference  27).  The  second  observation  Is  that  Uie 
spanwlse  and  chordwlse  momentum  thicknesses  are  rouglily  equal  for  the  flow 
cond'.tlons  at  which  the  modified  chordwise  method  is  useful  (Reference  28). 
Hence,  we  can  assume  that  i con  be  related  to  0 by  a tyqiical  flat  plate  skin- 
friction  fomula.  Here  we  pick  the  Dquire-Young  formula: 

1 (10) 

[8.89  log^y  fl..07‘>  R^,  r' 


where  Re  is  the  Reynolds  number  based  on  and  an  effective  viscosity  found 
using  the  same  method  that  Nash  proposed  for  the  effective  "chordwise  viscos- 
ity." The  resulting  values  of  are  in  rough  agreement  with  other  more 
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exact  prediction  methods. 

a.  Boundary  Layer  Transition  on  Swept  Wings 

When  the  wing  sweep  angle  is  greater  than  about  15°,  transition 
occurs  due  to  three-dimensional  instabilities  basically  different  than  those 
which  cause  transition  in  two-dimensional  flow.  Because  the  3-D  instabilities 
cause  transition  before  the  2-D  type  instabilities  are  fully  developed,  tran- 
sition criterion  based  on  2-D  flows  are  of  little  value  in  3-D  flows.  Adams 
(Reference  29)  has  reviewed  the  3-D  transition  criterion.  It  appears  that 
there  are  basically  two  parameters  of  engineeri^g  interest.  The  first  para- 
meter is  used  to  determine  if  the  bo-undary  layer  will  be  initially  turbulent 
at  the  stagnation  line  and  can  be  written  as, 


C*  = tan  A sin  A 


Rer 


d{V^oo)/<i{x/c)| 


= 0 


(11) 


where  x is  measured  in  the  x-direction  along  the  airfoil  and  s is  the  total 
arc  length.  If  C*  > 4 x 10^,  the  boundary  layer  may  be  turbxilent  initially. 
If  C*  > 1.4  X 10^,  the  boimdary  layer  will  remain  tiorbulent  if  tripped. 

The  second  parameter  can  be  used  to  predict  3-D  transition  if  the  flow 
is  laminar  at  the  stagnation  line.  It  is  given  as, 


X = 


’max 


V 

e 


(12) 


where 


(13) 


and  V = + vf  . 

e J e e 

the  maximum  crossflow 


6 is  the  local  boundary  layer  thickness.  W^.  is 

*'’max 

when  the  boundary  layer  profile  is  calculated  in  the 


stretwiline  coordinate  system.  iVfUisitlon  is  assumed  to  occur  wtien  X ^ \;rit 
>*i  17‘i.  In  nuuiy  cases,  is  t'oimd  to  occur  very  netu'  the  leading  edge. 

Attempts  to  estimate  X assuming  a FdiUiausen  prolUle  for  the  chordwise 
wul  sptuwiso  velocity  distributions  led  to  poor  correlation  with  the  sample 
case  presented  by  Adiuns  (Rel'erence  I’Q).  because  aji  accmaite  value  of  X 
be  found  only  if  the  crossflow  ■ ofile  is  actually  computed,  transition  based 
on  X has  not  been  Included  as  sm  option  in  the  basic  miKllfled  chordwise  bound- 
ary layer  prognun.  It  cim  be  computed  using  tlie  more  detailed  progrtuii 
described  in  Voluine  111. 

b.  Comparisons  with  Nasli  iuid  Tseng  Results  (Reference  75) 

Kigiu'es  *'  !uid  o show  tlie  comparison  l>etween  tlie  mixlified  chordwise 
methoil  results  mul  ti\e  l\illy  three-dimensional  (incompressible)  predictions 
of  Nasii  and  Tseng  (Reference  ) for  a 35^  untapered  swept  wing  wltli  miNACA 
o5^A017  airfoil.  Tills  is  tlie  case  for  wiilcli  tlie  modified  choj-dwise  metiiod 
was  originally  <iovel o>ped,  tuni  as  might  be  exiiected,  we  find  excellent  agree- 
ment. In  WgiU’e  5,  the  displacement  tliickness  agrees  closely  witli  tlie  l\ill 
3-n  calculation  tuid  the  difference  in  tlie  initial  sliear  stress  development 
may  be  due  to  the  difference  in  tlu'  t.rtuisitiivi  calculation.  Also  sliown  on 
tlie  figure  is  a civlcvilation  made  assuming  a two-dimensional  boimilary  layer 
witli  tlie  same  pressure  distribution.  The  7-n  calculations  underpredict  botJi 
displacement  tliickness  and  shear  stress  by  approximately  o-lt>X>.  Kigui*e  t>  is 
t'or  the  same  wing  at  mi  increased  lift  coefficient.  The  same  conclusions  mid 
trends  cmi  be  seen  for  this  case.  As  would  be  expected,  Uie  goixl  agreeivcut 
wrl  til  shear  stress  becomes  poor  near  the  trailing  edge.  The  deltastar  predic- 
tion remains  excellent,  however.  Wlien  the  flow  approaches  separation,  the 
simulation  of  3-n  effects  by  tlie  2-D  program  becomes  poor.  This  is  due  to 
tlie  tTuidanientol  difference  between  2-P  mid  i-P  flowl'ields  near  separation 
lines.  However,  these  results  indicate  tliat  the  growtli  of  the  cvvoy  is  con- 
fined to  regions  within  P or  3^'  of  sepai'aticui. 

c.  CiTiparl.son  witli  KxjH'rlmcnt  mid  Otlier  lYedlctlon  Methcxls 

Figures  7 and  8 provide  comparisons  with  one  of  the  few  cases  for 
whlcii  experimental  resxUts  m-e  available;  incompressible  flow  over  a swept 
wing  with  an  NACA  ti3j^-012  airfoil.  The  data  is  reproduced  fiami  NAOA  TN  2'->00 
(Reference  28).  Also  Includeil  in  the  prediction  of  Admns  (Reference  20)  tlirec- 
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dimensional  eddy -viscosity  method.  Figure  7 presents  the  zero  lift  case. 

The  agi-eement  between  both  prediction  methods  and  the  experimental  results 
are  good,  with  the  modified  chord  method  actually,  better  for  the  deltastar 
results.  Since  no  comparisons  for  skin  friction  are  available,  these  predic- 
tions of  the  modified  chordwlse  method  have  not  been  included.  Figure  8 pre- 
sents the  results  for  the  lifting  airfoil  case.  Also  included  on  this  figure 
is  the  prediction  of  the  Bradshaw  (Reference  30)  incompressible  3-B  program. 

Under  these  more  severe  conditions,  all  the  prediction  methods  underpredict 
the  displacement  and  momentvun  thickness.  The  3-U  Bradshaw  and  Adams  methods 
are  very  nearly  identical,  while  the  modified  chordwlse  method  lies  sliglitly 
below  the  predictions  of  the  other  methods.  Agreement  with  experiment  is 
better  for  the  momentum  thickness  thaji  the  displacement  tliickness.  The  luider- 
prediction  of  the  displacement  thickness  by  the  modified  chordwlse  metliod  com- 
pared to  the  3~D  Bradshaw  method  is  a trend  found  in  all  CLimparisons  to  date. 

d.  Comparisons  for  a Compressible  Flow  Case 

It  appears  that  no  experimental  results  are  available  foi'  compress- 
ible flow  over  an  infinite  swept  wing.  Thus,  Figures  9 and  10  present  a com- 
parison between  a 30°  swept  wing  compressible  case  given  by  Adams  (Reference 
29)  and  the  results  of  the  modified  chordwlse  method.  No  triuisition  location 
was  specified  and  some  numerical  experimentation  indicated  that  a transition 
point  located  at  .022c  produced  the  initial  conditions  indicated  in  the  results 
presented  by  Adams.  Tlie  displacement  thickness  predictions  are  in  good  agree- 
ment and  the  mcmentum  thickness  predictions  are  in  excellent  sigreement.  Fij^nire 
10  shows  the  ccmpai’isons  for  the  chordwlse  component  of  shear  stress.  Again 
we  see  that  the  skin  friction  prediction  deteriorates  as  separation  is  appiaxiched. 
Since  separation  actually  occurs,  the  shear  stress  I’esults  are  in  worse  agree- 
ment than  in  the  previous  cases.  Nevertheless,  the  separation  prediction  is 
within  3%  of  the  Adams  prediction.  Again,  note  that  even  near  separation  the 
deltastar  prediction  is  in  reasonably  close  agreement  witli  Adtuns.  AlthougJi 
the  30°  sweep  angle  may  seem  moderate,  it  is  tyiilcai  of  the  midciiord  sweep 
angles  encountei'ed  in  pj'actice. 
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I e.  Parametric  Study  of  the  Effect  of  Sweep,  Including  Comparison  With 

(■ 

I The  Bradshaw  3-P  I’rogram  (incompressible) 

j Figures  11  and  12  show  tlie  results  of  a parametric  study  in  wliicli  Uie 

! zero  sweep  pressure  distribution  of  the  Nash -Tseng  airfoil  and  Reynolds  number 

were  held  constant,  wliile  the  sweep  angle  A was  varied  from  0°  -*  (j0°,  and 
simple  sweep  theory  was  used  to  predict  the  actual  chordwise  pressui-e  distri- 
bution. Cases  were  run  for  botli  the  modified  chordwise  and  3-1'  Bradshaw  pro- 
grams. Figure  11a  shows  the  compai’lson  between  3-D  and  modified  chordwise 
programs  for  a case  at  zero  sweep.  At  this  sweep  angle,  the  modified  chord 
procedure  is  reduced  to  the  standard  2-D  compressible  Bradshaw  program.  This 
shows  that  the  two  different  prograjiis  protiuce  sliglitly  different  results, 
independent  of  the  modified  chordwise  procedure.  This  small  difference  shows 
the  typical  differences  between  different  mmterical  implementations  of  the 
same  basic  theory.  Figiu’e  lib  sliows  the  displacement  thickJiess  distribution 
for  A ‘j>0° . An  increase  In  the  difference  between  the  two  metliods  occurs 

with  the  modified  chord  metliod  vmderpredictlng  the  growtli  boundary  layer. 

Because  of  the  gocxl  tigreement  with  the  Nash  3-D  method,  these  calculations 
represent  a ci.imparison  between  tlie  Bradshaw  and  Nash  3-D  metliods.  Tlie  dis- 
agreement could  be  attributed  to  tlie  difference  in  the  basic  assumption  in  tlie 
methods.  Bradsliaw  adlows  tiie  siiear  stress  vector  to  act  in  a different  direc- 
tion than  the  velocity  vector's  nomial  gradient.  Nasli  assunes  tliat  botl\  act 
in  the  same  direction.  Tiius,  wo  have  obtained  va  qufuititative  measure  of  tlie 
difference  in  the  results  produced  by  the  two  different  metlitxis. 

Figure  12  shows  tlie  difference  between  the  modified  chord  method  luid  the 
incompressible  Bradshaw  3-D  method  as  a llmction  of  the  sweep  iuigle  for  the 
infinite  yawed  wing  case  of  Nash  !ind  Tseng  (Reference  The  ti’onds  are 

tlie  some  in  both  cases,  wtiile  the  quantitative  difference  grows  with  A as 
wovild  be  expected.  The  displacement  ;uid  manentum  thicknesses  are  consistently 
underpredicted  by  the  modified  chord  methods,  as  lu'e  the  shear  stresses. 

Also,  the  agreement  is  worse  at  the  trailing  c<.lge  luid  improves  forwiird  of  this 

station.  P'or  this  case,  it  ap{x?ars  that  the  two  metluxis  start  divei-ging  at.  i 

about  sweep  ,'uigle. 


f.  Flat  Plate  Study,  Including  Compressibility 


Figure  13  siiows  che  infinite  flat  plate  results  for  sweep  angles  of 
1*^  and  45^.  For  the  1°  sweep  case,  we  see  very  good  agreement  between  the 
2-D  and  3~D  programs.  We  can  see  a large  effect  on  the  displacement  thick- 
ness due  to  canpressibility  effects.  At  45°  sweep,  there  is  a larger  differ- 
ence (5%  at  the  trailing  edge)  between  the  two  methods.  Although  the  compress- 
ibility effects  are  still  quite  large,  they  are  not  quite  as  large  as  the  2-n 
case.  'Ehe  flat  plate  comparison  shows  tliat  the  shear  stress  effects  on  the 
boundary  layer  are  only  approximately  determined  in  the  modified  chordwise 


method.  Boundai'y  layers  wliich  develop  mainly  due  to  pressure  gradient  effects 
are  computed  more  accurately  than  those  wiiicli  are  dominated  by  shear  stress 
effects  wlien  the  modified  chordwise  method  is  used.  Fortunately,  the  boimdary 
layer  developnent  over  wings  is  dominated  by  the  pressui'e  gradient. 


g.  Comparison  of  Spanwise  Shear  Stress  Predictions 


The  modified  chord  method  was  developed  to  predict  the  chordwise 
properties  of  the  boundary  layer.  However,  an  estimate  of  the  shear  stress 
in  the  spanwise  direction  can  be  made  using  equation  (lOl,  wliere  O is  assumed 


to  equal  to  0^,  as  discussed  previously.  Figure  l4  sliows  a comparison  of  this 


estimate  with  the  predictions  of  the  detailed  computations  carried  out  by 
Nash  and  Tseng  (Reference  25).  The  good  agreement  is  surprising  and  perhaps 
forttiitous.  Figui’e  15  shows  tlie  canparison  with  tJie  results  presented  by 
Adams  (Reference  29).  The  3C^>  difference  seems  to  be  a more  reasonable  result 
in  view  of  tlie  approximations  used  in  obtaining  equation  (lo).  Hie  filled 
circles  show  the  results  of  a more  refined  equation  derived  using  the  siuiie 
approximations  made  in  obtaining  equation  (lO).  There  is  little  benefit  from 
the  more  refined  equation.  Finally,  note  that  the  spanwise  sheai'  stress  will 
be  sensitive  to  departvu’es  frcm  infinite  yawed  wing  conditions. 


h.  Comments  on  the  Perfomance  of  the  Method 


We  have  demonstrated  that  the  modified  chordwise  method  produces 
good  agreement  with  data  and  other  prediction  methods  for  a large  range  of 
conditions.  The  correct  qualitative  results  were  obtained  for  all  conditions 
tested.  When  assessing  the  applicability  of  the  methoti,  it  appears  tJiat  four 


basic  criteria  must  be  examined. 
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o Sweep  Aiigle 


o Separation 


o Compressibility 


o Spanwlse 

Inviscid  Flow 
Gradients 


■i? 


Results  obtained  here  indicate  that  the  method  pro- 
vides accurate  predictions  for  sweep  angles  of  40- 
45'^.  Beyond  this  sweep  angle,  the  predictions 
diverge  from  (but  maintain  the  same  trend  as ) tJie 
predictions  of  the  3-D  Bradshaw  program  (incompress- 
ible ) . 

In  the  present  scheme,  the  boundary  layer  calcula- 
tion terminates  when  t~  -•  0.  The  position  at  vrtiich 
this  occurs  is  called  the  chordwise  flow  separation 
line.  Since  the  flowfield  in  this  region  is  fully 
three-dimensional,  the  modified  chordwise  method 
must  fail  at  this  point.  This  breakdown  is  reflected 
in  the  sli^tly  premature  separation  prediction 
(as  much  as  3^  in  the  cases  studied). 

It  appears  that  for  typical  transonic  Mach  numbers 
compressibility  effects  do  not  affect  the  perfor- 
mance of  the  method.  A lack  of  experimental  data 
in  this  area  makes  a thoroxi^  evaluation  impossible. 
Tiie  difference  between  any  method  based  on  the 
infinite  yawed  wing  approximation  and  the  actual 
boundaiy  layer  on  a finite  wing  has  not  been  eval- 
uated. Tlie  basic  parameters  are  the  taper  ratio 
and  aspect  ratio  with  tip  and  root  non-uniformities 
of  interest  also.  For  wings  operating  near  a design 
point,  wliere  the  constant  isobar  sweep  is  maintained, 
the  si)anwise  gradients  are  low  and  the  prediction  of 
the  present  method  can  be  used  with  confidence.  It 
appears  that  an  experimental  program  that  investi- 
gates the  effects  of  taper,  root  and  tip  non- 
liniformlties,  and  off -design  pressure  fields  is 
required  in  order  to  assess  the  impact  that  neglect- 
ing these  effects  has  on  the  accuracy  of  the  pre- 
diction methods. 
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1.  Method  of  Computation 

The  results  presented  in  this  section  were  computed  using  pressure 
distributions  tabulated  in  Abbott  and  vcn  Doenhoff  or  taken  from  Adams'  paper 
(Reference  29)  in  the  compressible  case.  The  difference  between  the  arc  length 
and  the  distance  along  the  chord  line  x was  neglected.  The  good  results 
obtained  in  the  present  section  serve  to  validate  this  approximation  which  is 
consistent  with  the  method  developed  by  Bavitz  (Reference  4).  For  the  Nash 
and  Tseng  (Reference  24)  comparisons,  the  pressure  at  the  trailing  edge  was 
adjusted  to  avoid  flow  separation,  in  the  manner  indicated  in  their  paper 
(Reference  25). 

The  modified  chord  turbulent  b.l.  program  starts  with  a given  0/C  and  H. 
The  fV0.1y  3*D  program  starts  with  9/C  and  specified.  The  results  indicate 
that  these  different  starting  conditions  do  not  affect  the  boundary  layer 
development  after  a few  boundary  layer  thicknesses  downstream  of  the  transi- 
tion point.  In  the  application  of  the  code  to  tapered  wings,  the  character- 
istic sweep  is  taken  to  be  the  sweep  of  the  midchord  line.  Again,  a lack  of 
experimental  data  makes  the  choice  of  this  angle  diffic\ilt  to  verify.  This 
value  could  be  easily  changed  by  the  user  to  the  sweep  of  any  percent  chord 
line. 
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Figure  3.  Infinite  Yawed  Wing  Notation 
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Figure  12  Comparison  of  Modified  Chordwlse  and  3-D  Calculation  Methods 


Figure  13  Flat  Plate  Results  - 0°  and  U5°  Sweep 


Figure  15.  Spanwise  Component  of  Wall  Shear  Stress 


h.  COUPLING  OF  THE  INVISCID  AND  VISCOUS  CALCULATION  METHODS 
a.  General 


The  interaction  solution  that  combines  both  the  outer  inviscid  flow- 
field  and  the  thin  viscous  shear  layer  immediately  adjacent  to  the  body  has 
been  the  subject  of  numerous  investigations,  virtually  all  of  which  have  con- 
sidered two-dimensional  flow.  The  idea  of  iteratively  computing  the  potential 
flow  about  a body  with  an  additional  displacement  surface  due  to  viscous  effects 
was  apparently  originated  by  Prandtl  (Reference  31)*  Although  viscous  airfoil 
methods  were  continually  under  study,  the  development  of  the  aft-loaded  air- 
foil designs  for  efficient  transonic  flight  were  found  to  demand  viscous  methods 
in  order  to  predict  theoretically,  the  experimental  pressure  distributions 
accurately.  The  potential  for  increased  airfoil  performance  employing  advanced 
aerodynamic  concepts  bro\ight  about  a renewed  and  intense  interest  in  the  theor- 
etical prediction  of  transonic  flows.  The  development  of  accurate  2-D  bound- 
ary layer  programs,  followed  by  the  breakthrough  of  Murman  and  Cole  (Reference 
2)  in  inviscid  transonic  flow  predictions,  together  with  the  ever  increasing 
capabilities  of  computers  led  to  the  successful  development  of  the  first  2-D 
transonic  viscous  computer  programs  (References  U and  32,  for  example''. 

Aft -loaded  airfoils  produce  a rapidly  thickening  boundary  layer  on  the 
upper  surface  and  a thinning  boundary  layer  on  the  lower  surface.  The  result- 
ing hi^ly  asymmetrical  viscous  flowfield  acts  to  reduce  the  effective  aft- 
camber  associated  witli  the  rear-loaded  or  so-called  "supercritical"  airfoils. 

It  is  this  feattu'e  of  the  flowfield  which  causes  viscous  effects  to  become  so 
much  more  important  for  aft -loaded  airfoils  than  for  conventional  sections. 
Viscous  effects  arise  in  several  ways.  Figure  l6  shows  the  various  regions 
where  the  effects  differ  in  character.  Over  most  of  the  airfoil,  the  viscous 
effects  can  be  handled  by  adding  the  displacement  siirface  to  the  airfoil  and 
recomputing  the  potential  flow.  When  the  inviscid  flow  is  being  computed  by 
a small  disturbance  theory,  the  displacement  effect  is  treated  by  adding  the 
slope  of  the  displacement  thickness  to  the  slope  of  the  airfoil,  with  the 
resulting  boundary  condition  being  applied  at  the  axis.  For  inviscid  solutions 
that  treat  the  geometry  more  exactly,  either  the  flow  over  an  equivalent  body 
must  be  computed  or  the  flow  over  the  correct  gecmetrical  body  with  a normal 
velocity  specified  at  the  surface  is  computed  in  order  to  simulate  the  source 
flow  effect  that  the  boundary  layer  has  on  the  outer  inviscid  flow.  Tiie  latter 
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method  is  generally  preferred  because  only  a single  geometric  shape  must  be 
considered.  The  detailed  analysis  of  displacement  effect  treatment  has  been 
given  by  Preston  (Reference  33)  and  Lighthill  (Reference  34).  The  displace- 
ment surface  effect  is  a weait  interaction,  wherein  the  comparatively  small 
local  displacement  effect  of  the  boundary  layer  and  wake,  when  integrated 
over  the  whole  stirface,  leads  to  a reduction  in  lift  and  an  increase  in  the 
pressure  drag  on  the  wing. 

There  are  several  regions  on  the  airfoil  where  locally  strong  interactions 
occur.  Hie  interaction  at  the  trailing  edge  is  the  most  important.  This  is 
due  to  the  control  that  the  local  solution  at  the  trailing  edge  has  on  the 
global  circulation  and  therefore  lift.  Naturally,  the  detailed  solution  at 
the  trailing  edge  depends  both  on  the  boundary  layer  development  near  the 
trailing  edge  and  the  viscous  wake  just  aft  of  the  trailing  edge.  The  explicit 
computation  of  the  wake  and  its  effect  on  the  trailing  edge  interaction  has 
been  included  in  the  2-D  viscous  predictions  by  Hall  and  Firmin  (Reference  7) 
and  more  recently,  by  Melnik,  Chow  and  Mead  (Reference  5).  Melnik's  results 
show  explicitly  the  large  variations  in  pressure  across  the  interaction  region 
near  the  trailing  edge.  This  type  of  inviscid  pressure  distribution  was  first 
obtained  by  Pinkerton  (Reference  35 )>  who  reduced  the  circulation  to  the  mea- 
sured value  and  disregarded  the  Kutta  condition,  thereby,  obtaining  improving 
agreement  with  experimental  results.  The  wake  effects  can  be  split  into  thick- 
ness and  c\arvature  components.  In  the  present  method,  the  modification  of 
the  surface  slope  boundary  conditions  leads  to  the  implicit  modeling  of  the 
wake  as  a constant  thickness  lasrer  with  no  curvature  and  hence,  constant 
pressure  on  the  two  sides.  The  wake  curvature  appears  to  be  the  dominate 
effect  of  the  wake  contribution  to  the  trailing  edge  interaction.  In  the 
present  method,  the  wake  is  taken  into  acco\mt  by  the  empirical  treatment  of 
the  trailing  edge  problem  where  the  displacement  model  at  the  trailing  edge 
has  been  developed  by  correlation  with  experimental  data. 

The  viscous  region  at  the  foot  of  the  shock -wave  is  another  area  of  strong 
local  interaction.  As  long  as  the  shock-wave  is  not  strong  enougli  to  cause  a 
massive  separation  of  the  boundary  layer,  this  region  does  not  dominate  the 
overall  solution.  Often  the  inviscid  solution  "smears"  out  the  shock  jump 
due  to  the  lack  of  resolution  which  sirises  due  to  a coarse  mesh  (it  takes 
about  3 mesh  points  in  order  to  define  a shock  and  a typical  mesh  point  is 
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spaced  every  2^  chord  in  the  three-dimensional  methods).  For  strong  shocks, 
a "viscous  wedge"  can  he  placed  at  the  foot  of  the  shock  in  order  to  model 
the  strong  interaction  that  occurs  due  to  the  thickening  of  the  botindary  layer. 
We  note  that  in  transonic  flows,  the  direct  drag  due  to  the  presence  of  the 
shock  (wave  drag)  is  usually  small  compared  to  the  effect  of  the  thickened 
boundary  layer  on  the  pressiire  drag. 

An  additional  complication  of  the  3-D  transonic  flowfield  simulation,  in 
contrast  to  the  2-D  case, is  the  typical  existence  of  complex  3-D  shock  struc- 
txire  (Rogers  and  Hall,  Reference  36)* 

Another  important  aspect  of  the  interaction  solution  is  the  occurrence 
of  separation  on  the  airfoil.  Local  separation  vinder  shock  waves  appears  to 
be  adequately  treated  by  either  the  viscous  ramp  or  the  simple  thickening  of 
the  botindary  layer  at  the  foot  of  the  shock.  Small  regions  of  upper  surface 
separation  can  be  treated  by  a simple  empirical  extrapolation  of  the  displace- 
ment surface  slope.  A similar  treatment  is  employed  on  the  lower  surface 
when  separation  occurs  in  the  cove  region.  C'^ve  separation  can  be  termed  a 
"shallow"  separation,  where  the  region  of  separation  is  confined  to  a very 
thin  layer  such  that  it  can  continue  to  sustain  a pressure  gradient  (and  in 
fact , this  was  why  it  was  not  noticed  on  the  early  supercritical  airfoils ) . 

These  shallow  separations  usually  do  not  noticeably  degrade  the  performance  of 
the  airfoil.  Although  a rational  engineering  theory  is  needed  for  these  types 
of  shaillow  separations,  the  empirical  treatment  employed  in  the  present  pro- 
gram has  proven  reliable.  The  apparent  success  of  the  present  method  rests 
on  an  interesting  property  of  boundary  layers  on  the  verge  of  separation. 

It  can  be  readily  determined  that  when  a boundary  layer  is  about  to  separate, 
very  small  changes  in  the  imposed  pressiire  distribution  lead  to  large  changes 
in  the  displacement  thickness.  From  the  other  viewpoint,  this  fact  means  that 
large  changes  in  deltastar  result  in  only  minor  changes  in  the  pressure  distri- 
bution. This  property  is  largely  responsible  for  the  success  of  the  transonic 
viscous  flow  codes,  and  sho\ild  be  exploited  whenever  possible.  Although  a 
rational  engineering  theory  is  needed  for  these  types  of  shallow  separations, 
the  empirical  treatment  has  proven  acceptable  in  the  present  program.  There 
is  a regime  between  the  shallow  separations  that  can  be  treated  in  the 
present  method  and  massive  "constant  pressure"  separation  calculation  methods 
reported  recently  by  Maskew  and  Dvorak  (Reference  37)  and  Milgram  (Reference  38) 
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for  incompressible  flows,  where  a more  realistic  flowfield  model  for  separated 
flows  will  have  to  be  incorporated  in  order  to  make  reliable  predictions. 

For  more  extensive  discussion  of  the  physics  of  transonic  viscous- 
inviscid  interactions,  the  reader  is  directed  to  the  excellent  reviews  by 
Green  (Reference  39)  aJid  Pearcy  (Reference  UO),  et  al. 

There  have  been  very  few  investigations  of  transonic  viscous-inviscid 
interactions  on  finite  wings.  One  of  the  reasons  for  this  is  that  fully  three- 
dimensional  boundary  layer  calculation  methods  have  only  recently  become  available 
(Reference  11,  4l).  An  interesting  initial  examination  has  been  reported  by  Piers, 
Schipholt  and  van  den  Berg  (Reference  42)  for  subsonic  flow.  Hedman  (Reference 
12)  has  reported  an  interaction  method  which  treats  the  weak  displacement 
effects,  while  Kordvilla  (Reference  13)  has  also  used  the  Bailey-Ballhaus  invis- 
cid  program  to  investigate  the  effects  of  viscous-inviscid  interactions  in 
transonic  flow.  The  present  work  describes  the  first  unified  approach  to  an 
easy  to  use  and  reliable  engineering  tool. 

The  detailed  implementation  of  the  present  viscous  interaction  scheme 
can  be  described  as  follows.  First,  the  character  of  the  inviscid  solution 
must  be  established.  This  can  be  done  by  either  starting  from  a previously 
saved  solution,  or  by  i-unning  a number  of  inviscid  cycles  less  than  the  num- 
ber required  for  a converged  solution.  It  has  been  found  to  be  more  econcani- 
cal  to  recompute  the  boundary  layer  effect  every  forty  or  fifty  inviscid  cycles 
rather  than  to  reconverge  the  solution  completely  before  each  boundary  layer 
addition.  In  this  manner,  the  nvimber  of  inviscid  cycles  required  for  a vis- 
cous solution  can  be  held  to  less  than  twice  the  number  of  cycles  for  an  invis- 
cid solution,  and  more  typically  about  1-1/2  times  the  basic  number  of  itera- 
tions reqtiired  for  an  inviscid  solution. 

Starting  from  the  inviscid  solution,  the  boundary  layer  is  computed  at 
each  spanwise  station  using  the  modified  chord  method  described  in  Section  III. 

The  press\rre  distributions  are  modified  slightly  for  the  boundary  layer  cal- 
culation. Because  the  small  disturbance  theory  has  a singularity  at  the  lead- 

* 

ing  edge,  the  pressure  at  the  leading  edge  was  reset  to  the  actual  stagnation 
pressttre.  Any  pressure  cemputed  behind  the  leading  edge  above  a straight  line 
between  the  stagnation  value  and  the  tangent  to  the  pressure  distribution  with 
the  maximum  angle  is  reset  to  fall  on  this  line.  We  note  that  although  the 
two-dimensional  small  disturbance  theory  programs  often  produce  only  a single 
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pressure  coefficient  that  has  to  be  readjusted,  the  three-dimensional  solution 
sometimes  predicts  pressiire  values  above  the  stagnation  value  at  several  mesh 
points  downstream  of  the  leading  edge.  In  addition,  the  cove  pressure  distri- 
bution is  smoothed  several  times  in  order  to  reduce  the  artificially  large 
compression  of  the  flow  and  resulting  extreme  acceleration  near  the  trailing 
edge. 

The  boundary  layer  is  computed  using  the  midchord  sweep  at  each  span 
station  as  the  effective  sweep  angle.  Once  the  boundary  layer  is  computed, 
the  displacement  thickness  distribution  produced  by  the  boundary  layer  is 
interpolated  to  the  inviscid  solution  grid.  At  this  point,  the  displacement 
thickness  is  linearly  extrapolated  aft  of  the  separation  point,  and  the  separa- 
tion point  is  saved  for  futiire  use.  Once  the  basic  displacement  surface  is 
known  on  the  inviscid  grid,  the  trailing  edge  interaction  effects  are 
treated  by  linearly  extrapolating  the  slope  of  the  displacement  thickness  from 
a prescribed  point  slightly  upstream  of  the  trailing  edge,  to  the  trailing 
edge.  It  has  been  found  that  an  extrapolation  from  about  96^  or  97^  of  the 
chord  is  appropriate  for  many  airfoils.  This  factor  can  vary  slightly  between 
airfoil  sections  and  its  location  is  available  as  an  input  for  the  user  to 
adjust  if  he  desires.  No  extrapolation  is  required  on  the  lower  surface  of 
supercritical  airfoils.  The  wing  surface  slope  boundary  conditions  used  in 
the  inviscid  solution  are  altered  by  use  of  the  relaxation  formula: 


dx.  . 
input 


dz 

dx 


geometric 


liu  d6*^  + (i-u,) 

dx  dx  I ( l4 ) 


where  k denotes  the  iteration  niunber  and  u)  is  the  relaxation  factor.  Typically, 
o)  is  specified  as  .5,  although  for  a particular  case  its  optimum  value  may  be 
sli^tly  higher  or  lower.  This  procedure  forms  the  basis  of  the  viscous  inter- 
action procediire.  It  has  been  fo\md  beneficial  to  introduce  several  spanwise 
smoothings  at  the  trailing  edge  in  order  to  avoid  abrupt  changes  in  the  span- 
wise  boundary  condition.  This  is  especially  important  if  separation  is  present, 
as  described  below.  Examples  of  the  attached  flow  results  are  shown  in  Figure 
17  for  the  2-D  case  of  a supercritical  airfoil  for  which  data  is  available 
and  the  ONERA  M6  wing.  Note  that  for  attached  flow,  the  viscous  effects  are 
small  for  conventional  sections. 
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The  basic  parameters  that  the  user  can  adjust  for  a particular  case 
include  the  number  of  inviscid  cycles  between  boundary  layer  corrections,  the 
relaxation  factor,  and  the  extrapolation  location  for  the  deltastar  slope. 

This  permits  freedom  to  tune  the  method  for  a particular  case.  Along  these 
lines,  it  is  interesting  to  quote  the  recent  paper  by  Rose  and  Seginer  (Refer- 
ence U3),  who  discuss  the  coaputing  time  and  manhours  required  to  obtain  tran- 
sonic viscous  solutions: 

"Probably  the  hardest  to  estimate  and  least  talked  about 
machine  times  are  those  spent  before  the  final,  acceptable 
solutions  are  obtained.  It  has  been  our  experience  that  from 
k to  10  times  more  machine  time  is  used  to  obtain  the  final 
solution  (than  the  actual  canputing  time  of  a single  case ) . " 

Considerably  less  than  4- to  10  times  the  basic  run  time  should  be  used  in 
obtaining  results  with  the  present  code.  Nevertheless,  some  adjustments  of  the 
parameters  may  be  required  and  the  user  should  not  be  naive  about  the  typical 
results  published  by  the  research  community. 

After  the  basic  displacement  surface  modifications  are  made,  the  program 
determines  the  extent  of  separation  present.  If  the  separation  occiirs  ahead 
\ of  the  position  where  the  extrapolation  is  applied,  then  the  separation  region 

I'  . is  treated  as  described  in  Section  c.,  after  which  the  inviscid  solution  itera- 

' tion  is  reinitiated  with  the  modified  boundary  conditions.  A typical  conver- 

; ■ gence  history  is  shown  below  for  the  F-8  case,  starting  from  a previously  saved 

f ' solution: 


K 

Inviscid 

Iterations 

® Upper 

^ Lower 

1 

.412 

45 

1743 

228 

2 

.371 

45 

56 

45 

3 

.357 

45 

43 

21 

4 

.357 

45 

32 

17 

After  four  iterations,  this  case  was  considered  converged  because  the 
lift  was  effectively  constant  between  iterations. 
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Hgure  16  Some  Viscous -Invlscid  Interactions  in  Transonic  Flow 


b.  Shock  Boxmdary  Layer  Interaction  Effects* 


It  is  well  known  that  the  transonic  flow  over  airfoils  and  wings 
often  contain  embedded  shock  waves  which  have  subsonic  downstream  states. 

These  shocks  are  normal,  at  the  body  svirface  and  become  oblique  above  the  air- 
foil. For  most  aircraft  design  conditions,  a turbulent  boundary  layer  is 
present  on  the  body  s\arface  and  there  is  a resulting  shock  wave  --  turbulent 
boundary  layer  interaction.  Several  types  of  interaction  may  occur  (e.g.,  see 
References  36,  39  and  4o).  Reviewing  the  2-D  situation  first,  we  divide  the 
interaction  into  3 categories  (Figure  I8). 

o M ^ 1.05  - 1.10.  There  is  a weak  interaction  which  perturbs 
the  flow  field  only  locally  near  the  foot  of  the  shock. 

o M 1.10  - 1.30.  A strong  interaction  occurs  such  that  the 
shock  location  and  overall  strength  are  changed  from  the  inviscid  value. 

There  may  or  may  not  be  a localized  separation  bubble.  It  will  be  assumed 
that  if  a bubble  is  present,  reattachment  occurs  and  the  closed  bubble  is 
small  compared  to  the  boundary  layer  thickness.  In  general,  data  shows  that 
shocks  of  these  strengths  lie  below  incipient  separation. 

o MS  1.30.  There  is  large-scale  separation  initiated  at  the 
shock  wave  and  the  resulting  inviscid-viscous  interaction  is  strong. 

The  piirpose  of  this  section  is  to  provide  a phenomenological  method  for 
modeling  the  effects  of  the  strong  interaction  on  the  inviscid  flow  for  case  2. 
As  such,  the  method  has  been  incorporated  into  inviscid  transonic  computational 
methods  as  a first  step  in  a complete  inviscid-viscous  solution. 

(1)  A 2-D  Phenomenological  Model 

The  measTxred  pressure  increase  through  a shock  on  an  airfoil  or 
a wing  is  typically  less  than  the  normal  shock  value  for  the  corresponding  up- 
stream Mach  number.  Figure  I9  (see  Lomax,  Bailey  and  Ballhaus,  Reference  9) 
contains  a plot  of  experimentally  determined  pressure  increase  throu^  a shock 
wave  as  a function  of  the  shock  Mach  number  for  a wide  range  of  Reynolds  num- 
bers. The  curve  labeled  Rankine-Hugoniot  is  the  corresponding  normal  shock 
pressure  rise. 


* This  section  is  based  on  a consulting  report  by  Earll  M.  Murman. 
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Yoshihara  (Reference  4i+)  has  proposed  that  the  boundary  layer  thickening 
at  the  foot  of  the  shock  may  be  viewed  as  a "viscous -wedge"  as  shown  in  Figure 
20.  The  wedge  formed  by  the  displacement  surface  relieves  the  requirement  of 
a normal  shock  and  allows  the  shock  to  become  locally  an  attached  oblique  siiock. 
Although  the  viscous  wedge  model  is  not  necessarily  the  true  shape  of  the  dis- 
placement surface,  it  is  a convenient  and  consistent  explanation  for  the  effect 
of  the  boundary  layer  on  the  inviscid  flow. 

Figure  21  shows  schematically  a shock  polar  for  the  allowable  downstream 
pressiires  P2  and  viscous  wedge  angles  0.  To  adopt  Yoshihara's  model,  we  must 

require  that  the  oblique  shock  wave  be  attached,  that  is  0 s 0 . In  addi - 

max 

tion,  we  postulate  that  the  flow  downstream  of  the  shock  is  to  be  subsonic  so 

that  0 2:  0gQ^j,«  The  corresponding  pressure  jump  for  these  two  flow  angles 

is  shown  in  Figure  I9.  It  can  be  seen  that  the  measured  data  generally  fall 

between  these  two  limits  for  M 1.3.  For  larger  Mach  numbers,  it  can  be 

expected  that  the  flow  separates  and  consequently  is  not  of  interest  here. 

The  above  considerations  lead  to  the  following  phenomenological  procedure 

for  modeling  the  influence  of  the  boundary  layer  on  altering  the  shock  wave 

strength  and  location  (Figure  22).  The  airfoil  surface  is  allowed  to  have  a 

"rubber  wall"  which  deflects  at  the  foot  of  any  shock  which  impinges  on  the 

airfoil  surface.  The  initial  slope  of  the  deflected  displacement  sxirface  is 

chosen  as  0 . The  surface  is  then  bent  back  to  a required  slope  by  an  "ad 

m&Jc 

hoc"  shape  function.  The  deflecting  surface  continually  moves  with  the  sliock 

to  model  the  strong  interaction  feature.  At  this  stage,  there  is  no  attempt 

to  calculate  the  resulting  boundary  layer  displacement  surface.  This  would 

represent  a next  step  in  the  solution. 

The  choice  of  0 as  the  unique  angle  for  the  viscous  wedge  is  based  on 
max  10- 

the  following  reasoning.  The  inviscid  solution  tries  to  place  a normal  siiock 

at  the  airfoil  surface.  This  sliock  wave  is  formed  by  compression  and  expansion 

waves  emanating  from  the  surface  of  the  airfoil  and  sonic  line  upstream  of  the 

interaction.  This  wave  structure  is  basically  unaffected  by  the  interaction 

at  the  foot  of  the  shock  wave.  On  the  other  hand,  the  boundary  layer  must 

weaken  the  shock  to  obtain  a compatible  flow.  We  postulate  that  the  resulting 

shock  wave  is  the  strongest  oblique  and  attached  shock  wave  allowed;  i.e., 

0=0  . Since  the  surface  downstream  of  the  shock  must  bend  back  in  a stream- 

max 

wise  direction,  an  expansion  of  the  subsonic  flow  will  result.  It  is  not 
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inconsistent  then  that  the  measured  pressures  (Figure  19)  will  always  be 
equal  to  or  less  than 

0 

max 

The  choice  of  the  shape  function  needed  to  bend  back  the  displacement 
surface  remains  arbitrary  at  this  time.  Several  shapes  have  been  tried  with 
the  most  successful  to  date  being  a cubic  curve.  Figure  23  shows  the  com- 
parison with  an  exponential  shape  motivated  by  the  theoretical  analysis  in 
Reference  Let  the  airfoil  shape  be  given  by  z = 6F(x)  where  6 is  the 

thickness  ratio.  Then  the  modified  shape  is  given  by 

z = 6(F(x)  + Af(x)) 

where  (see  Figure  22), 


1 

0 

X < X 

s 

1 

AF(x)  = j 

0 Cli- 
max' 

Tl^  + 1 11^)? 

X^  ^ X ^ x^  . § 

— 5 9 

3 * max 

X - X 

T1  = -T— ^ 

x^  . § < X 

(15) 

where  x^  is  the  location  of  the  apex  of  the  viscous  wedge  and  5 is  a free 
parameter  defining  the  length  of  the  displacement  surface.  The  above  formula 
conforms  to  the  conditions; 

AF  = 0;  AF'  =0  at  x = x 
max  s 

AF'  -=0;  AF"  = 0 atx  = x +§. 

s 

A final  displacement  thickness  equal  to  §0_„„^3  results.  For  simplicity,  the 

nisuc 

above  fonnulation  is  termed  a "cubic-wedge"  displacement  function. 
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(2)  Implementation  of  Model  into  2-D  Inviscid  Calculations 


The  above  formulation  has  been  incorporated  into  the  inviscid 
small  disturbance  program.  At  the  end  of  each  relaxation  sweep  through  the 
field,  the  airfoil  surface  slope  is  modified  if  a "shock  point"  has  been 
detected  to  occur  at  the  airfoil  surface.  Referring  to  Figure  24,  the  point 
X.  is  taken  as  the  x-location  one  mesh  point  ahead  of  the  shock  point.  The 
formula  for  9 is  given  from  the  shock  polar  equation  for  the  small  distur- 
bance  equation  (see  Murman  and  Cole,  Reference  46)  as 


= + 


max 


Y + 1 


( 


(y  + 1)0X1  ■ ^ 


) 


3/2 


(16) 


with  the  plus  (minus)  sign  for  the  upper  (lower)  surface.  The  velocity  0^^^ 

is  calculated  from  a centered  difference  formula  for  the  location  two  mesh 

points  upstream  of  x^.  If  this  value  of  0x^  leads  to  a negative  argument  in 

the  above,  0 is  set  to  zero, 
max 

The  above  procedure  has  been  ■'■ested  on  several  airfoils,  some  presenta- 
tive  2-D  results  are  shown  in  Figure  25.  A value  of  § = .1  was  found  to  be 
adequate  to  move  the  shock  upstream  and  not  cause  too  severe  a reexpansion 
downstream  of  the  shock.  A case  was  computed  for  a shock  moving  upstream  and 
then  downstream  to  check  uniqueness  of  the  final  position.  Also,  the  unique- 
ness with  mesh  refinement  was  checked.  Some  calculations  have  been  done  using 

0 . instead  of  0 .No  major  differences  were  detected, 

sonic  max 

The  Murman  Bump  is  a model  which  participates  actively  in  each  inviscid 
iteration.  During  the  integration  of  the  method  into  the  code,  it  was  found 
that  the  original  idea  for  "hooking"  the  origin  of  the  ramp  to  a mesh  point 
inevitably  led  to  instabilities  in  the  numerical  solution.  The  cause  of  the 
instability  was  the  tendency  of  the  origin  to  alternate  between  mesh  points. 

By  interpolating  the  solution  locally  to  determine  the  actual  position  of  the 
sonic  point  in  the  smeared  shock,  the  origin  of  the  bump  could  be  attached  to 
a point  that  did  not  change  in  discrete  steps,  but  rather  continuously.  3t 
was  found  that  this  modification  uncoupled  the  physical  model  from  the  numerics 
and  led  to  a reliable  iteration  procedure.  In  addition,  it  was  found  that  the 
convergence  was  improved  if  the  bump  was  only  "turned  on"  after  the  basic 
solution  had  been  allowed  to  converge  to  the  point  where  the  character  of  the 
inviscid  solution  has  fully  emerged.  A modification  of  the  sdieme  to  use  a 
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"refined"  ramp  equation  which  more  closely  duplicated  the  strong  interaction  ! 

in  model  problems  led  to  results  essentially  the  same  as  the  origineil  poly-  , 

ncmial,  »diile  the  more  complex  formula  required  a considerable  increase  in  J 

computational  time  due  to  the  use  of  an  exponential.  Hence,  the  'polynomial 

bump'  has  been  retained  in  the  code  and  the  study  with  the  better  ptiysical  , 

model  can  be  viewed  as  verification  that  the  results  are  insensitive  to  the 
fine  scale  details  of  the  bump. 

The  results  obtained  from  these  runs  suggested  that  the  origin  of  the  \ 

ramp  shoxild  indeed  be  placed  sli^tly  upstream  of  the  sonic  location  in  the 
shock  as  the  original  method  specified  in  order  to  better  simulate  the  upstream 
influence  of  the  interaction.  Indeed,  it  appears  that  some  Improvement  can 
be  obtained  by  doing  this.  A nominal  shift  of  1 ~ 2%  chord  appears  to  improve 
correlation  and  has  been  left  as  a user  option  at  this  time. 

The  model  can  be  extended  to  three-dimensions  by  considering  the  three- 
dimensional  shock  relations.  The  shock  polar  for  a 3-D  shock  which  is  the 
weak  solution  of  the  equation 

(Ku  - ^ = 0 (17) 

is  (see  Murman  and  Cole,  Reference  46) 

(Up  + u^)^  (ug  - u^)^  + (Vg  - + (^2  - w^)^  = 0*  (18) 

Taking  y as  the  spanwise  direction  and  "z  as  the  direction  normal  to  the  air- 
foil (mean)  surface,  we  then  want  to  find  the  raaximrun  value  of  w^  for  u^, 

and  w^  = 0 given  quantities.  This  result  follows  from  setting  ^^w^/Du.,  = 0 

to  find  the  value  of  ru,  giving  maximxmi  w^.  The  result  is 
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vrtiere  the  plus  (minus)  sign  is  for  the  upper  (lower)  surface.  Thus,  the  maxi- 
mum value  of  occurs  wtien  = v^;  i.e.,  the  shock  is  normal  to  the  free- 
stream.  The  minimum  value  of  = 0 occurs  vrtien  the  deflection  in  the  span- 
wise  direction  is  the  maximum  allowable  turning  angle  to  the  given  upstream 
Mach  number. 

The  3-D  calculation  is  made  by  computing  x^  and  as  in  the  2-D  case, 

and  then  calculating  v^^  at  the  same  point  as  0^^^^  and  v^  at  the  next  point 

downstream  of  the  sliock  point.  This  then  gives  the  value  of  wo  for  the 

° “-max 

calculated  values  of  u^^,  and  v^. 

Figure  2o  sliows  some  results  for  the  3-D  case.  The  general  experience 
indicates  that  the  majority  of  the  work  done  with  the  present  computer  code 
will  be  carried  out  with  wings  at  or  neai’  design  conditions,  so  that  the  shock 
strengths  will  typically  be  minimal  and  the  explicit  use  of  the  present  mclel 
may  or  may  not  be  called  for,  as  sliown  in  the  figure,  where  the  changes  due 
to  the  ramp  were  small. 

A calculation  using  both  the  strip  boundary  layer  and  the  MiLcman  Bump 
raises  a question  as  to  the  proper  use  of  both  of  these  options  simultaneously. 
At  first  it  appears  that  this  would  produce  a "double  count"  of  the  shock 
boimdary  layer  interaction  because  the  strip  viscous  correction  also  simulates 
the  shock  bound.ary  layer  interaction.  However,  because  the  Bump  reduces  the 
strength  of  the  sliock,  the  displacement  surface  rise  at  the  shock  foot  is  sub- 
stantially reduced.  In  a similar  fashion,  the  displacement  surface  reduces 
the  strengtli  of  the  sliock  so  tluit  the  Munnan  Bi.imp  does  not  operate  with  the 
strength  wliich  it  would  have  if  the  viscous  displacement  correction  was  not 
in  effect.  Thus,  the  two  methods  tend  to  compensate  for  eacli  other,  and  the 
strong  shock  boundary  layer  interaction  is  not  dovibly  accounted  for  in  the 
present  method  when  both  options  are  in  effect.  The  precise  details  of  the 
complete  interaction  process  do  present  an  area  for  additional  study. 
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Fipxire  16  Typea  of  Shoclc  - 1‘urbulaat  Boundary  Layer  InteractiooB 


EXPERIMENTS 


O KACPRZYNSKI-OHMAN,  Re^  = 21  X lo” 

Z YOSHIIiARA-ZONARS-CARTER , Re^  = 30  X 10^' 
O PEARCEY,  Rec  = 2 X lo'' 

FLIGHT,  XI,  Re^  = 20  X lO'' 

A SINNOTT-OSBORNE,  Rec  ^ 2 X 10^' 

/ 


Figure  19  Pressure  Rise  Across  Shocks 


FORMS  WEDGE  LIKE  SURFACE  WHICH  PUSHES  SHOCK  FORWARD 
AND  WEAKENS  SHOCK. 


ngvire  20  Yosliihara's  Viscous  Wedge  Model 


RANGE  OF  POSSIBLE  0 


Figure  21  Shock  Polar 


Figure  2k  Inplementation 


oy  DATA  - STIVERS  AMES  2X2  FOOT’  TUNNEL 
THEORY,  SLOTTED  TUNMEL,  F=THIN  WALL 


Figure  25  Example  of  Viscous  Wedge  Effect 


STEVERS  DATA 


c.  Separation  Treatment 


In  the  analysis  of  transonic  flow  about  finite  wings  and  wing-body 
combinations,  it  is  important  to  be  able  to  treat  regions  of  shaJ.low  separated 
flow.  Most  airfoils  can  be  expected  to  operate  with  a region  of  separation 
occurring  on  the  upper  surface  near  the  trailing  edge,  this  is  especially  true 
of  airfoils  at  transonic  maneuver  conditions.  Current  supercritical  airfoils 
often  experience  an  additional  region  of  separation  on  the  lower  surface,  in 
the  cove.  It  is  important  to  treat  these  separated  flow  regions  in  such  a way 
as  to  prevent  the  computer  code  from  failing  when  it  detects  such  regions. 

The  procedures  used  for  treating  these  shallow  separation  regions  follow 
those  of  Bavitz  (Reference  4).  Because  of  the  equivalence  between  the  2-D 
displacement  thickness  and  the  3-D  displacement  surface  for  the  special  case 
of  the  infinite  swept  wing  approximation  employed  in  the  present  method,  the 
extension  of  the  methods  employed  in  2-D  to  3-D  is  self-consistent.  Care 
should  be  taken  in  using  these  ideas  when  a more  elaborate  3-D  boundary  layer 
calculation  is  employed.  On  the  upper  siirface,  shallow  separations  are  treated 
by  holding  the  slope  of  the  equivalent  inviscid  airfoil  constant  from  the  sep- 
aration point  to  the  trailing  edge.  On  the  lower  surface,  the  treatment  used 
in  the  presence  of  shallow  cove  separation  is  to  select  a representative  dis- 
placement thickness  shape  so  as  to  allow  the  iteration  procedure  to  continue. 

In  order  to  define  6*  for  the  separated  region,  four  points  are  chosen  along 
the  chord.  The  first  point  chosen  is  10%  chord  ahead  of  the  separation  point. 
The  second  point  is  located  8%  chord  ahead  of  the  separation  point.  The  third 
point  is  located  halfway  between  the  separation  point  and  the  trailing  edge 
and  the  last  point  is  located  at  the  trailing  edge  itself.  The  values  of  &* 
at  the  first  two  of  the  above  points  are  taken  to  be  those  previously  calcu- 
lated. The  value  of  6*  at  the  third  point  is  determined  by  incrementing  the 
value  at  the  first  point  according  to  an  empirical  equation  based  on  the 
pressure  coefficient  difference  between  the  two  chordwise  points.  The  relation- 
ship used  in  this  treatment  (which  was  developed  as  a result  of  an  extensive 
correlation  study  with  supercritical  airfoils)  is: 

\ - \)  - [ 'I’s  ■ 'I’l  ] • 


6* 


"1  + .033 


( 
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At  the  trailing  edge,  6*  is  taken  to  be  the  average  of  the  first  and  third  6* 
values.  A third  degree  polynomial  is  fit  to  these  four  values  of  6*  to  describe 
the  variation  of  6*  from  the  separation  point  to  the  trailing  edge.  As  Bavitz 
points  out,  the  most  arbitrary  values  in  the  model,  those  near  the  third  chord- 
wise  station  above,  are  defining  the  equivalent  inviscid  shape  in  the  cove 
region,  where  the  slope  is  approximately  zero.  In  this  region,  the  effect  of 
small  geonetric  perturbations  on  the  presstire  distribution  is  negligible. 

Under  these  circxmistances,  it  seems  that  qualitatively  reproducing  the  shape 
is  sufficient  to  allow  reasonable  calculations  to  proceed.  It  sho\ild  be  noted 
that  the  empirical  relationship  above,  was  determined  from  a detailed  study  of 
early  Whitcomb  airfoils. 

The  separation  treatments  described  above  are  designed  solely  for  engineer- 
ing use.  That  is,  it  is  understood  that  the  pressure  distributions  determined 
aft  of  the  separation  point  may  be  in  error.  However,  these  treatments  do 
allow  the  code  to  function  and  yield  accurate  Cp's  ahead  of  the  separated 
region  and  reasonable  engineering  approximations  beyond  the  separation  point. 
Thus,  the  code  can  continue  to  supply  reasonable  results  under  sitxiations  in 
which  the  lack  of  such  treatment  woiild  have  caused  the  code  to  halt  without 
yielding  any  information  at  all.  As  the  code  was  developed  in  modular  form, 
replacement  of  these  empirical  treatments  with  more  rigorous  treatments  as 
they  become  available  is  facilitated. 
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5.  BODY  MODEL* 

It  is  well  known  that  fuselage  effects  will  strongly  influence  the  wing 
pressure  distribution  in  transonic  flows.  In  order  to  obtain  a truly  useful 
engineering  tool  for  transonic  flow  analysis  of  actual  aircraft  configurations, 
sane  representation  of  the  body  must  be  included.  In  the  present  program,  it 

has  proven  effective  and  simple  to  model  the  body  by  providing  a constant  rec-  i| 

tangular  cross-section  boundary  condition  support  surface  upon  which  the  stream- 
wise  flow  inclination  angle  can  be  specified.  Figure  27  shows  the  way  in  which 
this  surface  appears  in  the  computational  space.  This  approach  provides  the 
user  with  a great  deal  of  flexibility  in  modeling  any  desired  geometric  effects. 

Because  the  detailed  slope  input  is  very  laborious,  the  program  has  a simple 
body  model  input  treatment  which  should  be  acceptable  in  most  cases.  The 
required  input  information  is  shown  in  Figure  28.  The  program  then  asstunes 
an  elliptic  cross-section  body  defined  by  the  input  body  lines.  Body  slopes 
are  generated  using  the  elliptic  cross-section  model,  together  with  the  input 
body  lines. 

In  order  to  provide  the  most  accurate  possible  body  effects  for  transonic 
wing  analysis  and  design,  the  slopes  specified  on  the  boundary  condition  support 
surface  should  not  be  the  actual  body  slopes,  but  rather  some  suitably  modi- 
fied set  of  slopes  that  take  into  account  the  difference  between  the  location 
of  the  actual  body  and  the  prismatic  surface  on  which  the  slope  boundary  condi- 
tion is  applied.  Slender  body  theory  can  be  used  to  determine  the  appropriate 
modifications  to  the  actual  boundary  conditions  required  when  they  are  trans- 
ferred to  the  computational  surface.  According  to  slender  body  theory,  the 
body  thickness  effects  should  be  represented  by  a source  with  strength  propor- 
tional to  rate  of  change  of  body  area.  S'  (X),  where  S(X)  is  the  cross-sectional 
area.  When  considering  the  bo\mdary  condition  to  be  specified  over  the  compu- 
tational surface  at  a particular  streamwise  plane,  we  require  the  source  strength 
over  this  surface  to  be  the  same  as  the  source  strength  for  the  actual  body  at 
this  streamwise  cross-section  location.  In  this  manner,  we  ensure  the  correct 
global  effect  of  the  body  on  the  rest  of  the  flow  field.  If  we  consider  the 

* The  authors  would  like  to  acknowledge  the  modifications  to  the  inviscid  pro- 
gram made  by  Allen  Chen  of  Boeing  and  Jack  Werner  of  the  Polytechnic  Insti- 
tute of  New  York  for  providing  the  appropriate  slender  body  corrections  to 
the  boundary  conditions. 
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body  cross-section  to  be  specified  by  the  cruve  C,  and  the  prism  to  be  speci- 
fied by  C,  then  the  requirement  that  the  source  strength  of  the  two  cross- 
sections  be  identical  becomes 


0=1  0=1 


so  that  the  simple  gradient  specified  on  the  prismatic  surface  is  related  to 
the  actual  boundary  conditions  by 


M 

9n 


prism 


Ar 


3n 


body. 


(23) 


This  procedure  is  illustrated  in  Figure  2g. 

Angle  of  attack  effects  due  to  body  incidence  can  be  treated  in  a similar 
manner.  In  this  case,  we  again  follow  slender  body  theory  by  considering 
equivalent  doublet  strengths  required  to  produce  the  correct  cross  flow  effect. 
Recall  that  the  doublet  strength  is  proportional  to  the  cross-sectional  area 
S(X). 

Since  the  body  slopes  have  already  been  specified  by  considering  the 
thickness  distribution,  the  "correction"  available  for  treating  the  body  on 
the  prism  surface  consists  of  modifying  the  actual  geometric  angle -of -attack 
to  an  appropriate  "effective"  aingle-of -attack.  The  effective  angle -of -at tack 
can  be  selected  by  requiring 


EFF 


S = a 


geo. 


S(X) 


(24) 


where  S is  the  prism  cross-sectional  area.  This  is  the  formal  statement  of 
doublet  strength  equality.  The  resulting  bovindary  condition  then  becomes  on 
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while  the  side  slopes  axe  not  modified  by  the  angle -of -attack  so  that  they 
are  given  by 

BOSS  " “ V®/  • 

SIDE 

Figures  30  and  31  show  the  agreement  for  a typical  body  case  and  also 
the  capability  to  handle  area  rule  effects.  The  results  indicate  the  strong 
influence  that  the  body  has  on  the  wing  pressures. 


Total  Input  for  Body  Simulation 

1.  Body  Top  Line,  Z^,(X^) 

2.  Body  Bottom  Line,  Zg(Xj^) 

Body  Side  Une,  Y^^Cx^) 

O 

U.  Body  Side  Line  Elevation,  Zj^(Xg) 


Figure  28  Simple  Body  Input  for  Infinite  Body 
Boundary  Condition  Surface 
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o. 


TRANSONIC  WAVE  UKAC  ON  SWEIT  WINC.S* 


a.  Oenei-al 

in  this  section  we  present  an  extension  to  the  modll’i 'd  sinali  dis- 
turbance theoiy  for  transonic  flows  by  suppl>'ing  the  farfleld  dra<,  forravila 
equivalent  to  the  farfleld  dra^  expression  for  the  classical  small  ilisturbance 
tlieory. 

All  integral  transonic  drag  relation  valid  for  unswept  wings  has  been  de- 
rived by  Murman  & Cole,  (Reference  UtO  , proceeding  directly  from  the  classical 
small  disturbance  equation, 

j^Ku  - (v  * 1)  1*]^  • ’ Vj.  = 0. 

This  small  disturbance  equation  was,  in  turn,  derived  from  an  asymptotic 
appivacli,  (Reference  U?)  , wiiich  made  the  ImiHicit  assumption  tliat  derivatives 
noimal  to  the  freestream  were  of  a different  order  of  magnitude  from  those  in 
tile  streamwise  dii'eotion.  This  assumption  is  implicit  in  the  form  adopted  for 
the  exiieuision  of  the  velocity;  i.  e. , 

“♦a  a a 

i = 1 (1  * €U  t i 6v  1^ 

U X y “ 

in  whicli  e and  6 are  assxmied  of  different  oi\ier  of  magnitude.  To  be  consistent 
with  this  assiunptlon,  any  sliock  surfaces  present  must  be  almost  noimal  to  the 
freestream.  Tliis  condition  la  generally  violated  over  swept,  wings. 

bailey  Sc  RalUiaus  (Reference  48)  have  modified  the  classical  small  distur- 
biuice  equations  by  writing, 

^Ku  - (y  * 1)  ^ • e (V  - 3)  ItJ  ^ ’ (1-«(y-1)u)  - 0 

adding  tenns  wlilcli  allow  the  resulting  equation  to  reduce  to  the  "correct  small 
disturbance  equation"  for  2-P  swept  wings.  The  proposal  of  this  equation  for 
the  colculatioj;  of  3-b  swept  wings  Implies  that  derivatives  in  Uie  sweep  direc- 
tion ai'e  sraivLl  compai'cd  to  those  normal  to  the  sweep  direction. 


• This  section  is  based  on  a Grumman  Aerospace  Corporation,  Aeroilynamlcs 
Section  Report  (390-77-0?,  12  May  1977)  by  Jack  E.  Werner,  Consultant  to 
(irumman . 
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For  the  swept  case,  an  expansion  analogous  to  that  of  References  45  and 
46  would  set; 


g = [cos  V + eu^  +...]  + [sin  v + 6V^]  + i^  6w 


where  u , v,  are  speeds  normal  and  tangentieil  respectively  to  the  sweep  direc- 
n t 

tion,  and  v is  the  sweep  angle.  This  would  lead  to  a small  disturbance  equa- 
tion for  the  swept  region  which,  unfortunately,  would  not  hold  where  shocks 
are  almost  normal  to  the  freestream  (i.e.,  over  the  center  section  of  a swept 
wing).  Rather  than  pursue  this  approach,  Bailey  and  Ballhaus  have,  in  effect, 
proposed  the  first  term  of  a "composite"  expansion  by  including  terms  which 
are  small  in  the  unswept  region  but  which  become  large  enough  to  supply  a 
correction  in  the  swept  region. 

We  now  seek  an  integral  relationship  appropriate  for  the  computation  of 
drag  from  a solution  of  the  Bailey -Ballhaus  equation.  If  we  proceed  in  a 
manner  similar  to  Murman  and  Cole  (Reference  46),  we  quickly  encounter  a serious 
difficulty  in  determining  a divergence  form  associated  with  the  drag  and  whidi 
would  lead  to  a surface  integral  for  the  drag.  To  resolve  this  problem,  we 
shall  first  explore  the  exact  procedure  of  which  Murman 's  and  Cole's  represents 
a first  order  approximation. 

b.  Exact  Expression  For  Lift  and  Drag 


An  exact  surface  integral  for  the  force  on  a wing  surface  may  be 
obtained  by  considering  the  integral, 


(pq)  dV 


over  a volume  V over  all  space  excluding  the  thin  regions  surroiinding  the 
wing,  the  wake  and  any  shock  surfaces  present.  Noting  that  the  continuity 
equation  renders  the  above  integral  identically  zero,  it  may  be  transformed 
with  the  aid  of  Gauss'  Theorem  to  obtain: 


(pq)  dV  s 0 


yp  q (q  • n)  dA  - 
A 


P (q  • V ) q dV 


(27) 
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Fi’oin  the  steady  muiicntum  equation, 


- p ( q • V ) q VP 


{VS) 


luid  a iXirther  application  of  Gauss'  Theorem,  we  arrive  at  the  equation. 


(pq")  dV 


0 


pq  (q 


dA 


(29) 


wiiere  the  surface  A * W * S ^ake  ' represents  the  closed  surface  formed 

by  siu'roundlng  shocks,  the  siu'face  W about  the  wir.g,  S w:.ike  around  ttie  air- 
foil wake  iuid  G the  remainini:^  surface  required  for  closure.  Now, 


1, 


is  the  not  force  on  tlie  wlnj;.  .Also,  we  iiave  ttie  conditions  q • n 0 on  Uie 
wing  and  over  the  w.ike  surl’ace,  wlille  p is  continuous  across  tlie  wiuke.  These 
conditions  :vllow  us  to  write  equation  (29^  in  tlie  fonii: 


(pq''  dV 


0 


(pn  • pq  (q  • lO ) ilA 


The  X component  of  tlie  first  surface  Integral  on  the  riglit  cjui  be  shown  to 
^deld  tlie  induced  drag  wliile  the  second  surface  integral  is  identified  as  the 
wave  drag,  if  a solution  of  the  classical  small  distui-biuice  equation  is  sub- 
stituted into  the  rijiiit  luuid  side  of  equation  iiO^  iuid  only  the  first  signifi- 
cant order  of  magnitude  retained,  the  results  yield  the  drag  integral  of  Muniuui 
d-  Cole.  'The  procediu'e,  in  effect,  followed  by  Mmmian  iv  Cole  was  to  exp-md 
the  left  luuid  side  of  equation  (30)  directly,  retain  the  first  significiuit 
tenns  :uid  apply  Gauss'  Theorem  to  the  result.  We  shall  proceed  in  a similar 
nuumer . 
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c.  First  and  Higher  Order  Eicpansions 


A combination  of  energy  and  continuity  equations  yields  a convenient 
form  of  the  continuity  equation  from  which  derivatives  of  the  density  p have 
been  eliminated  (Reference  *+9). 


- fi.  w.  a 
/ Pco 


^2/  o 

V(^  ^2u/) 


together  with  the  energy  equation, 


* (¥) 


Upon  introducing  the  expansions: 


q/U  = i (l  + su  ^ u^ 
^ c«  X 


> . . ) ^ 6V^ 


p''p^  = 1 + eo  ^ e Op 


1 - M*"  = eK 


and  the  following  relations  obtained  from  the  first  order  small  disturbance 
theoiy : 
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u 


+ (y+1)  -2  (y+1)  (v2)  u 

8 


Kuu2  - (y+I)  u ^2 


(40) 


-i 

The  term  in  e of  this  previous  expansion  was  obtained  by  Murman  & Cole  directly 
from  the  classical  small  disturbance  equation.  Since  the  first  two  terms  of  the 
coefficient  of  originate  with  the  additional  terms  of  Bailey  and  Ballhaus, 
the  expansion  of 


I 


corresponding  to  the  Bailey-Ballhaus  equation  would  be: 


1 V . I ^ 


3 

a)  1 e 
X 


2 3 

[K^  - (Y+1)  :i  - _I  , 

2 3 2 I X 


-e 


2u  Vt  • u + (y-1)  u“  • 


(»+l) 


Unfortunately,  the  last  term  of  equation  (4l)  cannot  be  written  in  divergence 
form,  and  the  volume  integral  of  this  last  expansion  of, 


cannot  be  converted  completely  into  a surface  integral.  As  a result,  this 
particular  expansion  of  the  left  hand  side  of  equation  (30)  cannot  match  up 
with  any  corresponding  expansion  of  the  right  side  such  as  might  be  obtained 
by  substituting  a solution  of  the  Bailey-Ballhaus  equation  into  the  right  hand 
expression. 

It  would  thus  appear  that  the  form  represented  by  equation  (4l)  cannot 
properly  be  associated  with  any  surface  integral  expression  for  the  drag.  To 
obtain  the  proper  expansion  of. 


for  swept  wings,  it  is  first  observed  that  for  unswept  wings,  the  Bailey- 
Ballhaus  terms  are  gnwn  and  may  be  included  in  the  expansion  without  affect- 
ing its  magnitude  to  within  0 (e^).  On  the  other  hand,  when  swept  shocks  are 
present,  these  terms  represented  by  the  two  leading  terms  in  the  coefficient 
of  in  equation  (40)  become  large  and  must  be  considered  along  with  the  € 3 
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f 

I 


J 


term.  But,  if  these  first  two  terms  are  now  included,  so  should  the  third  term 
be  included  since  it  is  almost  equal  in  absolute  value  to  the  second  term.  In- 
clusion of  this  third  term  now  allows  the  resulting  expansion  to  be  written  as 
a divergence;  i.  e. , 


= 0 


ih2) 


This  establishes  a vector  M analogous  to  that  of  Mumnan  & Cole,  but  appropriate 
to  the  computation  of  drag  for  swept  wings. 


d.  The  Drag  Integral 

Integrating  the  pressure  over  the  wing  surface  represented  by  z = 
6F(x,y),  an  expression  for  the  drag  is  obtained. 


D = 


dxdz 


W 

1/3 

where  [ ] ^ represent  the  jump  across  W and  z"  = 6 

condition  on  the  wing  in  the  form 


(43) 


Writing  the  boundary 


V 

1 + eu+.. 


V (1  - eu) 


(44) 
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The  drag  integral  becomes  with  the  substitution  for  p from  equation  (37) 


(45) 


Since , 


V * M d 


V 


= 0 


from  equation  (42),  an  application  of  Gauss'  Theorem  to  the  volume  bounded  by 
Sc,  W,  S„  and  substitution  of  the  results  into  equation  (45)  yields: 


V 


• n d S 
c 


S 


0 


(46) 


were  [ M]  represents  the  jump  in  M across  the  shock  surface  S^,  and  from 
equation  (42) , 


- ^ 

M = i.. 


K I - (Y*l)  f 1 
^2 


Vj  U (1-€U) 


(47) 
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The  integral  over  S may  he  interpreted  to  yield  the  induced  drag  and  when 
integrated  over  the  surface  as  x “ yields  the  same  result  as  given  by 
Murman  &Cole,  Reference  {h6). 


e.  Wave  Drag 

We  now  turn  to  the  second  integral  in  equation  (46).  But  before 
dealing  with  it  directly,  an  appropriate  form  of  the  kinematic  shock  Jump 
condition ; 


[q]  X n^  = 0 on  S(x,y,z)  = 0 


must  be  developed,  letting  n represent  the  direction  normal  to  the  sweep 
direction  and  t the  direction  tangent  to  it,  equation  (48)  becomes. 


in  ^ ^ in 


+ S 

1 


[u„]  s - [V^l  g . LV^]  X s = 


where  T represents  a direction  in  the  plane  transverse  to  n. 

1 

Noting  that  the  last  term  is  of  higher  order  than  the  proceeding  one,  we  may 
set  the  first  term  equal  to  zero  resulting  in: 


V s . ^ 

an  [un3 


A shock  jmp  condition  is  derived  by  writing  the  small,  disturbance 
equation  in  the  form: 


V • Q = 0 


and  applying  Gauss'  Theorem  to  the  region  surrounding  the  shock.  This 
results  in: 


[Q]  • n^  = 0 = 


cVi  * r»,]  • V s 


Introducing  the  kinematic  jump  condition,  equation  (50),  we  have: 


[Q^]  Cu^]  + = 0 


or,  since  the  secular  product  is  independent  of  axes  rotation, 


[Q^]  [uj  + • [v^]  = 0 


For  the  Bailey-Bal.l haus  equation,  we  have. 


2 2 \ _ 

Q = i Ku-(y+1)  ^ + e(Y-3)  + V (l-e(Y-l)u) 


Sc  S 


Introducing  this  value  of  Q into  the  shock  Jump  condition,  equation  (54),  and 
employing  results  from  the  calculus  of  Jumps  (Reference  46)  ; we  obtain: 


[^Ku-  u^]  [u]  + LV^J^  - cLv/j  [u]  = e(v-l)  [V^J^  <u>  (56) 

It  is  of  interest  to  note  that  setting  y = 2 in  the  terms  multiplied  by  e 
in  equation  (53),  yields  the  NLR  small  disturbance  equation  (Reference  50). 
Alternatively,  the  NLR  equation  is  obtained  from  equation  (39)  Ly  applying 
the  same  considerations  which  led  to  equation  (42)  for  M;  i.  e. , retaining 
the  first  three  terms  in  the  coefficient  of  instead  of  only  the  first  two. 
In  this  respect,  the  NLR  eqxiation  is  seen  to  be  wholly  consistent  with  M of 
equation  (42)  from  the  point  of  view  of  the  expansion  process. 

We  are  now  in  a position  to  evaluate  the  wave  drag  integral: 


°WAVE 


//' 


n dA  = 


// 


Thus,  with  the  aid  of  the  kinematic  condition,  equation  (50), 


J = [M]  • n = 
c 


Vs 


^S/Sn 

IW 


[M  1 [u]  + [HrJ  • [V 


9S/Bn 

IvW 


Substituting  M from  eq,uation  (U7)  and  again  using  the  calculus  of  jumps: 
/ 


This  is  considerably  reduced  by  the  introduction  of  the  shock  jump  condition 
from  equation  (56), 


J = - (Til 
12 


[u] 


3 tu] 


+ e (y-2) 


e [V^]‘ 


[ul^ 

U (60) 


In  keeping  with  previous  comments,  we  shall,  assxmie  that  disturbance 
velocities  in  the  sweep  direction  are  small  compared  to  disturbance  veloci- 
ties normal  to  the  sweep.  Thus,  if  v is  the  sweep  angle,  we  have. 


= u sin  V + V cos  v 0 

(61) 

u = u cos  V - V sin  v 
n 

(62) 
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in  conse<iuence  of  which  we  obtain : 


[u  ] = 
n 


Substituting  equations  (6o)  and  (63)  into  equation  (57)  for  'the  wave  drag. 


\ave  ^ (v+l) 

6^^737^  " 12 


[u]^  cos  V 


SS/an 


dA 


f flv.fM  . 

s I I dx  dz 


+ e (y-2) 


II 


[ui 


cos  V dA 

JWl 


^ow  projection  dS^  of  the  shock  surface  on 

tZ  or  "sweep"  plane  while  cos  v the  further  projection  of 

on  the  x7  plane.  If  we  reinterpret  the  sweep  plane  to  be  the  local  shock 
surface,  the  wave  drag  expression  may  be  written  in  the  final  form: 


^AVE  = (y+1) 
54/3p  y 12 


II 


[u]^  dx  d'z  + c(y 


^■// 


<u>^ 


dx  dz 


II 


[VJ'"  [u] 


dx  dz 


In  the  neighborhood  of  normal  shocks,  Cv^]  is  small  compared  to  [u]  and 
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the  c ntribution  from  the  last  two  integrals  in  equation  (65)  becanes  negli- 
gible and  the  Murman  and  Cole  result  is  recovered.  It  is  also  worth  noting 
that  if  Q is  obtained  from  the  NLR  equation,  the  second  tern  of  equation  (65) 
drops  out.  Near*  a swept  shock  [V^]  may  be  large  and  the  contribution  fron 
the  last  two  terms  is  correspondingly  significant. 

Eqviation  (65)  is  the  fairfield  drag  formula  which  must  be  used  when  the 
MSD  equation  is  employed.  The  integration  was  not  implemented  in  the  present 
study.  This  equation  would  have  to  be  implemented  on  the  crude  exterior  grid, 
and  the  complexity  of  the  final  equation  would  probably  lead  to  large  nmeri- 
cal  errors  when  implemented  on  such  a coarse  grid.  Thus,  the  calculation  of 
drag  in  the  present  program  has  been  based  completely  on  surface  pressure 
integration. 


7.  fULLY  THREE-DIMENSIONAL  BOUNDARY  LAYER  CALCULATION 


In  order  to  Increase  the  power  of  the  present  transonic  viscous  method, 
a new  fully  three-dimensional  finite  wing  boundary  layer  program  has  been 
written  by  John  Nash  and  Roy  Scruggs  of  Sybucon,  and  is  described  in  the  self- 
contained  document  issued  as  Volume  III  of  the  present  report.  The  data  set 
required  to  execute  the  program  is  generated  by  the  transonic  analysis  code, 
either  from  a direct  inviscid  calculation  or  after  the  interaction  solution  is 
obtained  using  the  infinite  swept  wing  strip  viscous  approach.  The  new  pro- 
gram is  a fully  implicit  numerical  solution  of  the  3-0  boundary  layer  equations. 

A compressible  S-D  turbulent  energy  equation  baaed  on  the  incompressible  equa- 
tion derived  by  Nash  (Reference  51)  is  used  to  obtain  turbulent  closure.  As 
a result,  we  expect  that  this  program  could  potentially  be  more  accurate  thaji 
the  eddy -viscosity  model  programs  that  are  now  beginning  to  appear  (Reference 
Ul).  In  addition,  input  data  requirements  are  greatly  improved  compared  to 
the  previous  explicit  program  developed  by  Nash  and  Scruggs  (Reference  2h). 

In  this  section,  we  outline  the  features  of  the  method  so  that  the  user  of  the 
programs  can  anticipate  the  complete  capability  provided  by  the  set  of  pro- 
grams described  in  the  three  volumes  of  the  present  report. 

The  3~D  BL  program  computes  the  laminar  or  turbulent  boundary  Itvyer  devel- 
opment over  a finite  swept  wing  of  arbitrary  planform  and  thickness  in  com- 
pressible steady  flow.  The  program  is  structured  in  order  to  be  iuiienable  to 
interactive  calculations,  althougli  the  full  coupling  was  not  carried  01  t 
in  the  present  effort.  The  wing  surface  is  segjnented  into  quasi -quadrilateral 
panels  and  a locally  ortliogonal  set  of  curvilinear  coordinates  is  constiucteil 
following  a conical  development  of  the  local  region.  Stresuiiwise  cuts  are  made 
over  the  pltmfom  at  arbitrary  spanwlse  locations  tuid  the  number  of  surface 
points  on  each  cut  is  the  same.  The  calculation  then  proceeds  from  tlie  pliui- 
form  leading  edge  as  turbulent  or  laminar.  The  governing  equations  are  finlte- 
differenced  with  resjiect  to  the  local  surface  normal  coordinate  y,  and  the 
locally  transverse  coordinate  z,  so  that  in  the  local  y-z  pltme,  the  difference 
equations  are  implicit.  The  solution  is  then  obtained  at  each  x (stretuiiwlse^ 
station  in  succession.  This  chordwise  forward -marching  is  applied  with  arbitrary 
step  length  proceeding  from  the  leading  eiigc  to  the  trailing  etlge  or  until 
separation  is  encountered  at  semie  spanwise  location.  Provision  is  made  to 
step  past  local  separations  by  altering  the  local  pressure  distribution,  If 


desired.  Transition  location  must  be  specified  by  the  user. 

Further  flexibility  in  the  new  program  has  been  achieved  by  adding  the 
thermal  energy  equation  to  the  set  of  governing  equations.  Althou^  the  pro- 
gram computes  adiabatic  boundary  layers  by  specifying  the  adiabatic  wall  temp- 
erature, a simple  modification  would  allow  the  program  to  compute  boundary 
layers  with  arbitrary  wall  temperature.  Another  refinement  is  the  solution 
of  the  finite  difference  equations  all  the  way  to  the  wall,  removing  the 
limitations  associated  with  the  matching  to  the  law  of  the  wall  region  as  is 
usually  done  with  turbulent  kinetic  equation  methods.  As  a result  of  this 
refinement,  the  calculation  of  both  laminar  and  turbulent  flows  can  be  carried 
out  within  the  same  basic  framework. 

Grumman  has  added  a few  additional  calculations  and  revised  the  output 
sli^tly,  in  addition  to  verifying  the  program  for  a few  simple  cases  for 
which  comparisons  can  be  made.  Checks  of  the  niimber  of  chordwise  stations 
required  to  obtain  a converged  solution  indicated  essentially  no  difference 
between  predictions  using  26  and  4l  chordwise  points.  The  program  is  dimensioned 
for  a maximm  of  60  chordwise  and  30  spanwise  stations.  It  appears  that  21 
points  across  the  boundary  layer  is  sufficient.  The  computation  of  the  stan- 
dard integral  thicknesses  was  added  as  well  as  the  incorporation  of  a summary 
of  the  results,  an  option  to  punch  some  of  the  output  for  plotting  and  an 
option  to  provide  several  levels  of  output,  depending  on  the  type  of  informa- 
tion desired.  With  the  exception  of  these  minor  additions,  the  program  has 
been  maintained  in  the  form  delivered  to  Grumman  by  Sybucon. 


SECTION  11 


APPLICATIONS 


1.  INTRODUCTION 


In  this  part  of  the  volume,  we  give  a large  number  of  examples  of  the 
application  of  the  program  to  specific  cases.  The  intent  is  to  show  typical 
accuracy  that  any  user  should  expect  to  achieve.  Although  a knowledge  of  the 
underlying  theory  is  desirable,  it  is  not  necessary.  Experience  at  Crumman 
has  indicated  that  good  results  can  be  obtained  by  users  who  do  not  know  the 
details  of  the  methodology  as  long  as  a few  rules  of  thumb  are  followed.  How- 
ever, it  is  Important  to  understand  that  the  present  programs  are  intended  to 
provide  an  experienced  aerodynamlclst  with  inform.atlon  required  for  aerodynamic 
design  studies.  We  expect  that  this  program  will  be  but  one  of  a number  of 
codes  that  form  the  aerodynaml c 1st ' s bag  of  tools.  As  such,  the  user  should 
constantly  check  the  answers  tor  consistency  and  overall  "correctness."  When 
the  program  does  not  produce  answers  that  appear  reasonable,  an  error  in  tiic 
input,  or  a simple  misunderstanding  In  the  operation  of  the  program  can  usually 
be  found  to  explain  the  anomalies  in  the  results.  We  will  warn  the  user  repeat- 
edly to  study  each  run  carefully,  applying  good  aerodynamic  judgement  to  the 
evaluation  of  the  predictions.  If  this  procedure  is  followed,  the  program 
will  consistently  provide  predictions  of  at  least  the  same  quality  as  those 
presented  in  this  part  of  the  report. 

The  ten  survey  cases  shown  here  have  evolved  over  the  course  of  the  study. 
They  constitute  a wide  range  of  configurations  which  would  be  very  typical  of 
the  actual  application  of  the  code.  Some  of  these  cases  have  formed  the  basis 
for  the  original  demonstration  of  various  methods,  while  other  cases  have  been 
used  to  demonstrate  the  use  of  the  method  to  aircraft  design  projects. 

The  F-8  and  TACT  configurations  are  two  of  the  most  advanced  aircraft 
designs  for  which  flight  test  data  is  available.  This  section  shows  the  typi- 
cal correlation  between  flight,  wind  tunnel  and  theoretical  predictions  of  the 
present  method. 


2.  EXPERIENCES  AND  LIMITATIONS  OF  THE  CODE  --  HOW  TO  SUCCESSFULLY  RUN  THE  CODE 

Although  extremely  useful  results  can  be  obtained  with  the  present  set 
of  computer  programs,  the  code  does  have  some  geometric  restrictions  and  aero- 
dynamic limits.  In  this  section,  we  review  some  of  these  limitations  and  dis- 
cuss strategies  tliat  can  be  employed  when  the  normal  solution  procedures  fail 
to  obtain  a satisfactory  result.  The  geometric  limitations  do  not  necessarily 
mean  that  the  program  cannot  be  used  for  these  cases,  but  rather  that  the  aero- 
t.iynamicist  will  have  to  supply  a more  imaginative  model  to  the  program  and  j 

judge  the  results  with  additional  care,  in  order  to  extract  the  maximum  possi- 
ble value  from  the  program. 

Cieomctric  features  wliich  the  program  does  not  explicitly  handle  are: 
o Multiple  Lifting  Surfaces  (including  canards  and  wlnglets) 
o Wings  with  Planform  Discontinuities  (such  as  snags) 
o Wing  Dihedral 
o Inlets/No2::les/Stores 
o Multi -Element  Airfoil  Sections 
o Fine  Details  of  the  Fuselage  Geometry 
Aerodynamic  limitations  include: 

o Mach  numbers  less  than  one 
o Attached  flow  over  most  of  the  surface 
o No  suction  blowing  effects  on  boundary  layer 
o No  heat  transfer  effects  on  boundai’y  layer 
o Mid-chord  wing  sweeps  less  than  about  70*^ 
o Taper  ratios  greater  than  .05 

In  addition  to  these  general  limitations,  a few  specific  configurations 
may  require  special  treatment.  Wlien  these  cases  occur,  considerable  additional 
flexibility  is  allowed  so  that  it  should  prove  possible  to  obtain  a solution. 

These  silternate  strategies  can  be  grouped  into  several  categories: 
o FYunilies  of  Solutions 
o Geometry  Simplification 
o Inviscid  Iterative  Wocedures 
o Mesli/Mapping  Alteration 
o boundary  Layer  Relaxation 

before  illscussing  the  actuiU  solution  methexis,  we  discuss  some  input 
problems.  Extreme  care  must  be  exercised  in  setting  up  the  data  sets  for 


^X) 


1 


I 


any  computational  aeroiiynsunics  type  computer  prof^rtuii.  A large  number  of 
problems  attributed  to  the  earlier  versions  of  the  code  by  I'l'ustrated  users 

have  been  traced  to  an  input  error.  In  the  iiresent  program,  the  potential  j 

for  user  errors  has  been  minimi c.ed  by  cai'etlil  design  of  the  Input  procedui’e,  j 

I 

internal  checks  on  the  cfuisistency  of  the  data  tuid  a graplilcs  package  wtiicli 
provides  grapliical  verification  of  all  the  input  geometry  as  veil  as  the  com- 
puted results.  The  graplilcs  package  CfUi  be  invaluable  and  we  urge  all  users 
to  make  the  effort  required  to  adapt  this  part  of  the  progriun  to  their  own 
system.  Although  a number  of  internal  checks  have  been  included  in  order  to 
keep  the  user  from  exceeding  the  limits  of  the  various  input  values,  undoubt- 
edly it  is  still  possible  to  make  inputs  that  are  not  corrected  or  flagj:T;ed  by 
the  prograai.  Common  errors  have  included  specifying  the  mesh  in  sane  type  of 
noi'mall^ed  units  and  the  planform  in  pliysical  dimensions  or  vice  versa,  ajid 
specifying  the  wrong  Index  for  grid  locations  sucli  as  the  first  mesh  line 
past  the  tip.  To  repeat,  experience  has  showi  that  a carefully  set  up  data 
set  will  reduce  problems  to  a minimiun  luui  reward  tlie  user  with  reliable  results. 

Families  of  solutions  can  sometimes  be  used  to  obtain  results  that  are 
proving  difficult  to  compute  by  the  standard  methods.  In  this  scheme,  a solu- 
tion that  can  be  readily  computed  is  saved  and  the  saved  solution  is  used  as 
the  starting  guess  for  another  solution  with  a minor  increase  in  difficulty. 

Typically,  this  approach  can  be  used  to  obtain  higlier  smgles -of -attack,  more 
extreme  twist  distributions,  or  extreme  NlacJi  numbers.  Several  Increments  can 
be  used  to  obtain  the  final  desired  condition.  Even  wtien  no  difficulty  is 
being  experienced  in  obtaining  a ivu’tlculai’  solution,  tills  metluxi  ciui  lead  to 
i-educed  computing  times  luid  hence,  more  econamical  computing.  This  strateg.v 
can  be  partlculsirly  ecoiuxiilcal  dui'ing  design  studies,  where  the  effects  of 
minor  geixnetrlc  clionges  are  being  studied. 

Cieometry  simplification  Ciui  correspond  to  both  the  planform  and  airfoil 
section  treatment.  Peeause  the  small  disturbance  theory  has  a singularity  at 
the  leading  edge,  the  resulting  large  flow  gradients  near  the  leading  edge 
can  cause  numerical  solution  problems  for  some  cases.  The  program  identifies 
tlie  grid  indices  that  locate  the  numerical  difficulties  in  the  grid,  so  that 
it  is  possible  to  determine  if  the  problem  occurs  at  mesti  points  about  the 
leading  edge.  Wlien  this  occurs,  the  mesh  can  be  inspected  to  Insure  that  the 
mesh  distribution  is  smooth.  The  airfoil  slopes  should  be  inspected  to  veriiy 
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that  the  proper  airfoil  ordinates  have  been  input  and  the  (smooth)  resulting 
slopes  are  reasonable  (note  that  the  input  ordinates  eire  spline  fit  in  order 
to  interpolate  to  the  mesh  locations  and  prescribe  the  slopes,  requiring  an 
accTjrate  specification  of  the  airfoil  ordinates).  If  the  ordinates  are  correct 
and  the  mesh  is  distributed  smoothly,  then  the  user  has  the  option  of  reduc- 
ing the  leading  edge  slope  in  a systematic  manner  by  employing  Riegels ' Rule . 
Riegels'  Rule  is  a method  that  renders  the  first  order  thin  airfoil  solution 
valid  near  the  leading  edge  in  incompressible  flow.  Van  Dyke  (Reference  I9) 
has  investigated  its  application  along  with  a number  of  other  methods  to 
improve  thin  airfoil  theory.  Unfortimately,  Riegels'  Rule  has  no  theoretical 
foundation  for  transonic  flows.  It  does,  however,  often  allow  the  solution  to 
proceed  without  numerical  difficulty  and  has  the  interesting  property  of  remov- 
ing much  of  the  mesh  dependence  of  the  solution  around  the  leading  edge.  How- 
ever, it  should  be  used  with  caution  because,  for  transonic  solutions,  it 
appears  that  the  resulting  solution  may  be  inacciarate,  particularly  on  the 
lower  surface. 

Another  geometric  difficulty  can  arise  if  the  user  is  trying  to  prescribe 
non-streamwise  wing  tips.  This  can  lead  to  an  extreme  bending  of  the  computa- 
tional mapping.  In  this  case,  the  simplest  improvement  would  be  to  model  the 
tip  as  streamwise,  with  an  equivalent  planform  area.  If  this  is  not  accept- 
able, the  user  may  have  to  specify  the  mapping  himself.  Attempts  to  include 
strakes  as  part  of  the  planform  will  also  lead  to  problems,  because  the  strake 
constitutes  a small  fraction  of  the  wing  area,  in  a location  where  several 
extreme  mesh  bendings  are  required.  The  resulting  loss  of  geometric  resolution 
in  the  mesh  bending  process  will  produce  a result  that  does  not  adequately 
model  the  strake  and  will  introduce  oscillation  in  the  solution.  Most  gloves 
can  be  handled  without  difficulty.  If  difficulties  are  encountered  while 
investigating  planforms  with  severe  geometry  variations,  some  numerical  experi- 
ments should  be  run  in  order  to  isolate  the  particular  geometric  feature  that 
is  causing  the  difficxxlty.  Once  the  problem  area  is  located,  a review  should 
be  made  of  the  modeling.  It  may  be  that  a slight  gecanetric  simplification  is 
perfectly  acceptable.  If  not,  then  the  span  mesh  may  have  to  be  concentrated 
about  the  region  of  interest  and  a combination  of  the  methods  discussed  in 
this  section  employed  in  order  to  obtain  the  solution. 
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The  inviscid  iterative  procedures  can  be  changed  to  allow  for  a smaller 
value  of  ID  and  a larger  value  of  the  damping  coefficient  for  Some  trials 

at  other  Mach  numbers  and  angles -of -attack  may  help  clarify  the  source  of  the 
problem . 

For  example,  if  the  inviscid  iteration  procedure  is  diverging,  a lower 
angle-of -attack  or  Mach  number  may  result  in  a converging  result,  or  it  may 
help  to  identify  the  location  in  the  mesh  at  which  the  solution  procedure  is 
failing.  Cnee  this  location  is  defined,  inspection  of  the  local  mesh  and  geom- 
etry will  probably  identify  an  irregularity  that  can  be  corrected.  In  many 
cases,  this  is  sufficient  in  order  to  obtain  a solution.  The  print  plot  of 
the  inviscid  convergence  history  should  be  examined  for  each  run.  Both  the 
exterior  and  interior  mesh  iteration  histories  are  displayed  and  the  particular 
mesh  that  the  problem  originates  on  can  be  determined.  This  can  be  used  as  a 
guide  to  determine  which  mesh  should  have  the  iteration  peirameters  changed. 

Mesh  and  mapping  alterations  have  been  mentioned  previously.  The  program 
automatically  generates  a mesli  that  is  a result  of  considerable  experience 
with  the  code.  The  user,  however,  can  input  his  own  mesh  and  mapping.  Wlien 
this  step  must  be  used,  it  should  be  done  with  considerable  care,  following 
the  instructions  contained  in  the  user's  manxial. 

A rule  of  thumb  that  works  well  for  troublesome  cases  in  all  instances 
investigated  to  date  is  that  for  regions  wtiere  the  mapping  is  bending,  such 
as  at  glove-wing  intersections,  the  spanwise  mesh  must  be  more  closely  spaced 
than  elsewhere  on  the  wing.  For  some  extreme  cases  that  are  nearly  pure  delta 
wings,  the  mapping  must  be  bent  near  or  beyond  the  tip  in  order  that  the  § = 0 
line  remains  in  front  of  the  § = 1 line  in  the  physical  space.  If  this  condi- 
tion is  violated,  the  program  will  give  a warning  and  stop.  To  attempt  to 
obtain  a solution  vuider  these  conditions  will  only  result  in  a catastrophic 
failure  of  the  iteration  procedure. 

The  boundaxy  layer  iteration  can  be  controlled  by  several  metliods.  The 
degree  of  underrelaxation  of  the  boundary  layer  can  be  controlled,  the  origin 
of  the  slope  extrapolation  can  be  varied  independently  on  the  upper  and  lower 
surface,  and  the  number  of  inviscid  iterations  can  be  atijusted.  In  general, 
however,  the  viscous  solution  will  bresik  down  when  there  is  excessive  flow 
separation  present  on  the  wing.  The  exact  degi’ee  of  separation  that  can  exist 
and  still  obtain  a convergent  iteration  will  depend  on  the  particular  case. 

On  supercritical  wings,  the  iteration  has  converged  for  separation  lines  that 
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occur  at  7<yf,  on  the  lower  surface  and  8^-90%  on  the  upper  surface.  Naturally, 
' the  larger  the  separated  region  hecomes,  the  worse  the  approximation  to  the 

actvial  solution  will  become. 

From  the  above  discussion  it  should  be  clear  that  the  more  experience 
that  is  gained  with  the  program,  the  more  useful  the  code  will  become.  While 
the  program  should  operate  satisfactorily  with  minimal  experience,  the  code 
can  be  used  to  compute  some  exceptioneLLly  challenging  cases  in  the  hands  of 
an  experienced  user. 


KXvWl.F.S  OF  CONFICURATIONS  OOMi’UTKI)  AND  TYl’lCAL  ACKEKMKNT  WITH  KXI>ERIMF;NT 

This  program  has  boon  tested  on  a luimher  ol  ooni igurat ions  (see  Table  1) 
covering  a wide  range  ot  flight  conditions  and  geometries.  These  ''ave  included 
aspect  ratios  from  2 to  over  h , wing  leading  edge  sweeps  from  JO*’  to  over  dO*’, 
glove  sweeps  in  excess  of  7s*'  and  taper  ratios  from  under  .025  to  0.6.  In 
fact,  *'ven  a forward  swept  wing  has  been  successfully  treated.  The  various 
flight  conditions  handled  have  included  Mach  numbers  from  0.81  to  0.99  and 
.»ngles-of-attack  over  a range  from  O*’  to  8.5*.  A representative  sample  of 
such  cases  has  been  included  in  tills  section.  Many  of  these  cases  are  classics, 
having  been  used  as  the  original  demonstration  cases  for  some  of  the  inviscid 
programs  in  use.  Most  of  these  cases  consist  of  relatively  simple  "rese.irch" 
type  geometric  configurations  for  which  transonic  pressure  data  is  generally 
available.  Two  of  the  cases,  the  low  aspect  ratio  and  forward  swept  wing,  .iri> 

Included  to  show  the  capabilities  of  the  program  even  though  no  experiment.il 
data  is  available  for  these  cases.  These  cases  have  been  quite  useful  even 
though  some  of  the  experimental  data  contains  wind  tunnel  wall  effects  and 
many  of  the  wings  were  not  designed  for  transonic  conditions.  The  present 
effort  shows  the  need  for  establishing  a new  set  of  baseline  experiment.il  cases 
for  use  in  aerodynamic  tool  development.  Most  of  the  good,  modern,  experimental 
pressure  data  is  closely  tied  in  with  specific  aircraft  projects.  .As  such,  the 
configuration's  geometric  complexity  is  often  of  a nature  to  make  it  difficult 
to  perform  a basic  code  development  comparison  study.  In  addition,  this  type 
of  data  Is  often  classified  .ind  can  not  be  widely  distributed.  We  have,  how- 
ever. included  two  such  cases,  the  F-8  and  TACT  aircraft,  which  will  be  exten- 
sively discussed  in  the  next  section.  These  cases  contain  the  most  widely 
avall.ible  supercritical  wing  data.  The  declassification  of  the  F-8  progr.im  w.is 
p.ir  t icularly  timely  in  this  regard. 

The  results  are  presented  in  each  case  with  some  of  the  basic  geometric 
data  for  tiie  case.  The  det.iiled  geometric  data  reqviired  to  actually  run  those 
cases  can  be  found  in  the  references  for  each  case.  Comparisons  of  experimental 
pre.ssure  distributions  witli  tlie  program  predictions  are  Included.  .All  of  the  pre- 
sented results  were  computed  with  the  final  computer  code,  using  the  autonut ica I ly 
generated  meshes  and  mapping.  The  computations  employed  the  FCR  differencing  of 
the  MSD  Equation,  with  inviscid  solutions  used  for  the  cases  with  conventional 
sections  and  viscous  flow  solutions  for  cases  with  supercritical  sections. 
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a.  ONERA  M6 

The  first  case  listed  in  Table  I is  the  ONERA  m6  (Reference  52).  The 
ONERA  m6  is  an  isolated  wing  case  with  a conventional  type  ONERA  "D"  airfoil 
section.  This  case  was  studied  at  a Mach  number  of  0.84  and  an  angle -of - 
attack  of  3°.  The  wing  has  a leading  edge  sweep  of  30°  and  a trailing  edge 
sweep  of  15.8°.  The  aspect  ratio  for  this  case  was  3-8  and  the  taper  ratio 
was  0.56.  The  basic  configuration  can  be  seen  in  Figure  32a.  Figure  32b 
shows  the  three-dimensional  pressiire  distribution  over  the  wing.  As  can  be 
expected,  the  shock  location  moves  forward  as  one  moves  outboard  on  the  wing. 

The  first  span  station  for  which  experimental  data  was  available,  was  at 
T]  = 0.2.  The  closest  span  station  at  which  the  code  gave  predictions  was  at 
T)  = 0.182.  Figure  32c  shows  the  comparison  of  experimental  data  with  the  code's 
predictions.  As  can  be  expected,  due  to  the  fact  that  the  program  prediction 
is  at  T1  = 0.182,  the  shock  location  is  predicted  farther  aft  than  the  measure- 
ments indicate.  The  lower  siorface  predictions  agree  q\iite  well  with  the 
experimental  data.  In  fact,  at  each  span  station  of  the  ONERA  solution,  the 
lower  surface  agreement  with  data  is  excellent.  The  next  span  station  (Figure 
32d),  'Hgjjp  - 0.44  and  = 0.436,  shows  excellent  agreement  with  experiment 

over  the  leading  edge.  The  upper  surface  leading  edge  pressure  peak  is  pre- 
dicted excellently  with  respect  to  both  location  and  hei^t.  The  shock  loca- 
tion is  predicted  approximately  5^  ahead  of  the  location  shown  in  the  measure- 
ments. The, upper  surface  leading  edge  pressiire  peak  is  again  accurately  pre- 
dicted at  the  next  span  station,  T]  =0.65  and  Tl  = 0.655  (Figure  32e). 

exp  pre 

The  measurements  show  evidence  of  two  shocks,  at  x/c  = O.I5  and  x/c  = 0.45. 

The  second  shock  is  predicted  quite  well,  both  in  amplitude  and  location.  The 
first  shock,  however,  is  not  predicted  as  well.  There  is  a relatively  rapid 
drop  in  pressure  from  the  leading  edge  peak,  however,  the  shock  does  not 
appear  explicitly.  This  could  probably  be  corrected  by  increasing  the  number 
of  points  along  the  chord  above  the  44  used  at  this  span  station,  and  possibly 
increasing  the  angle -of -attack  used  in  the  program  slightly.  The  next  station, 
\xp  ” \)re  ~ 32T)  again  predicts  the  leading  edge  peak  and  the 

second  shock.  However,  the  first  shock  is  still  missed  and  the  same  corrections 
could  probably  be  made.  At  the  two  outboard  span  stations  = 0'90> 

7]^^^  = 0.909  and  = 0*95j  = 0.945)  the  leading  edge  peak  is  still 

well  predicted  but  the  shock  is  not  as  sharp  as  indicated  by  the  data. 


i 


especially  at  Tl  = O.95.  The  pressure  dip  follcwiiig  the  shock  is  well  pre- 
dicted as  to  height  but  the  predicted  location  is  approximately  10^  ahead 
of  that  indicated  by  the  measurements.  Jameson  and  Caughey  (Reference  21) 
have  also  performed  inviscid  calculations  on  the  ONERA  m6  with  a namerical 
solution  of  the  exact  transonic  potential  flow  equation.  Their  predicted 
shock  location  agrees  with  the  data  somewhat  better  than  this  code;  however, 
their  code  was  run  at  a slightly  higher  angle-of -attack  and,  as  seen  prev- 
iously, a slight  increase  in  angle -of -at tack  would  seem  to  be  indicated  for 
this  code  also.  Due  to  the  concentration  of  points  near  the  leading  edge, 
the  Jameson-Cauj^iey  code  predicts  the  second  pressure  peak  better  than  this 
code,  however,  an  increase  in  the  number  of  points  along  the  chord  would 
increase  resolution  and  possibly  show  this  better. 

b.  RAE  "A" 


The  second  case  listed  in  Table  I is  the  RAE  "A"  wing-body  configura- 
tion. The  experimental  pressure  data  (Reference  53)  was  obtained  at  a Mach 
number  of  O.9O  and  an  angle-of -attack  of  1°.  Figure  33a  shows  the  configura- 
tion used  for  this  code.  The  airfoil  section  used  in  this  case  is  the  RAE 
101  conventional  type  section.  Figure  33h  shows  the  3-D  pressure  distribution 
over  the  wing.  The  predictions  were  obtained  at  an  angle-of -attack  of  1°  and 
runs  were  made  at  Mach  numbers  of  O.9O  and  0.91.  The  3-D  pressure  distribu- 
tion shown  is  for  the  Mach  O.9I  case.  Figure  33c  gives  the  pressure  distribu- 
tion at  = 0.1667  and  = 0.182.  This  span  station  for  the  predictions, 

%re  ~ 0*182,  is  the  first  span  station  outside  the  body  in  the  mesh  generated 
by  the  code.  As  can  be  seen,  the  agreement  between  prediction  and  experiment 
is  not  very  good  at  this  span  station.  However,  this  is  due  in  large  part 


to  the  difference  between  the  prediction  and  experimental  span  stations.  In 
this  region,  the  effect  of  the  body  i^  greatest,  and  the  spanwise  pressure 


gradients  ore  the  highest.  Observing  the . predicted  gradient  between  "Hp^g  = 

0.255  and  Tl  = 0.182,  it  is  apparent  that  moving  fram  H = 0.182  inboard 
pre  pre 

to  Tl  = 0.1667  would  cause  both  the  upper. and  lower  predicted  pressures  to 


drop.  Thus,  if  the  program  utilized  Tl^^g  = O.I667,  the  predictions  would 
more  closely  follow  the  experimental  results.  If  the  details  of  the  pressure 


distribution  near  the  body  are  desired,  the  user  must  input  his  own  y-mesh. 


concentrating  the  mesh  near  the  wing -body  juncture.  It  is  interesting. 


I 
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however,  to  note  that  increasing  the  Mach  mraber  fran  O.9O  to  0.91  causes 
the  shock  location  to  shift  aft  about  10^  and  brings  it  to  the  correct  loca- 


tion. Figure  33d  shows  the  pressure  distribution  at  = 0.25  and  = 

0.255.  Here,  the  agreement  between  experiment  and  prediction  is  much  better. 


Note  again  that  increasing  the  Madi  number  to  O.9I  gives  excellent  agreement 
with  data.  The  shock  location  again  moves  aft  to  the  correct  location  seen 
in  the  data.  Also,  the  lower  surface  pressures  drop  to  the  correct  values. 


The  next  span  station  sliown  (Figure  33^)  is  at  “ O.U.  Here  again. 


the  predictions  show  better  agreement  with  data  at  Macli  O.9I.  At  T] 


exp 


0.6, 


Tlpre  = 0.618,  (Figure  33f)>  the  agreement  is  again  better  at  Mach  O.9I, 


although  the  upper  surface  pressure  predictions  ore  sojnewliat  below  the  data. 
The  lower  surface  shows  excellent  agreement  with  data.  The  final  outboard 
span  stations  (Figures  33g  smd  33h)  show  similar  characteristics  with  better 
agreement  at  Mach  0,91.  These  results  clearly  demonstrate  the  Mach  number 
sensitivity  to  be  expected  in  transonic  flow  calculations.  Thus,  in  any  com- 
parison of  predictions  with  exiierlmental  data,  tlie  Macn  number  used  in  tlie 
experimental  work  must  be  known  accurately. 


c . 


L54hi8 


The  third  case  is  the  T,54lll8  (Reference  54).  This  case  consists  of 
a NACA  I'SAOO'i  wing  mid-mounted  on  a cylindrical  body  (Figure  34a).  The  wing 
is  twisted  from  0,4°  at  the  wing-body  juncture  to  -4.12°  at  the  wing  tip,  Tiie 
flight  conditions  used  for  this  case  were  a Mach  number  of  0.95  !Uid  an  angle- 
of-attack  of  4°,  Fifpare  34b  shows  the  3-0  pressure  distribution  for  the 


L54hi8  wing.  The  first  span  station  - ^*818)  (Figure  34c) 

shows  the  same  characteristics  observed  in  the  RAF  "A"  configuration.  This, 


too,  is  a situation  where,  if  the  span  station  used  for  the  predictions  had 
been  closer  to  that  used  for  the  experimental  results,  tlie  results  would  have 
been  closer.  Again,  wlien  close  to  the  body,  where  the  sponwlse  pressure 
gradients  arc  the  lai’gest,  comparisons  between  prediction  and  experiment 
require  the  span  stations  used  to  be  the  some.  Figure  34d  sliows  better  agree- 
ment with  the  data.  In  fact,  except  for  the  upper  sui'face  leading  edge,  the 
agreement  is  excellent.  The  leading  edge  peak  does  not  reach  the  heiglit  of 
that  observed  experimentally.  The  shock  location  prediction  agrees  with  the 
experimental  location.  In  Figure  34c,  we  again  see  excellent  agreement  except 
for  the  iieit^ht  of  the  pressure  peak  at  the  leading  eilge.  The  next  span 
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station  = ^pre  ~ 0.80)  (Figure  3^f)  shows  excellent  agreement;  however, 

the  upper  surface  experimental  results  shows  a second  pressure  peak  and  shock 
location  at  x/c  = 0.35.  This  peak  is  not  duplicated  in  the  predictions.  If 
; the  wing  twist  was  increased  sliglitly,  the  predictions  could  conceivably  agree 

I somewhat  better.  The  final  outboard  section  (Figure  3^6 ) shows  the  same 

characteristics  as  the  previous  span  station. 

d.  TN  0712 

The  fourth  case  investigated  was  the  TN  D712  wing-body  configuration 
(Reference  55).  The  fli^t  conditions  for  this  case  was  a Madi  number  of 
0.94  and  an  angle -of -at tack  of  2.4°.  The  airfoil  section  used  in  this  wing 
is  a MCA  65AOO4.  Figure  35a  shows  the  basic  geametry  of  this  configuration. 
The  3-D  pressure  plot  is  seen  in  Figure  35b.  The  first  span  station  shown 
^^exp  ~ \)re  ~ 0-1^2)  (Figure  35c),  indicates  excellent  agreement  between 

experiment  and  prediction.  The  only  point  at  which  agreement  is  not  canplete 
is  at  the  trailing  edge.  The  predictions  are  frati  an  invlscid  solution;  a 
viscous  solution  would  probably  show  better  agreement  at  the  trailing  edge 
due  to  separation  treatment.  The  code  fails  to  predict  the  upper  surface 
i leading  edge  pressure  peak  at  the  next  span  station  = ’'pre  ” ^•*^) 

(Figure  35d).  However,  the  rest  of  the  distribution  at  this  station  agrees 
quite  well.  In  fact,  over  the  rest  of  the  span  (Figures  35e,  35f  and  35g)j 
the  code  continues  to  fall  to  predict  the  leading  edge  pressure  peak.  However, 

I except  for  the  final  span  station,  the  shock  location  and  magnitude  is 

: accxirately  predicted.  The  outboard  span  station  (T1„,_  = 0.95>  = 0.945) 

shows  relatively  poor  agreement  over  the  first  45^  of  the  section.  The  lower 
surface  shows  a peak  at  x'c  = 0.4  which  is  not  reproduced  in  the  predictions. 

In  addition,  the  upper  surface  distribution  is  not  very  well  reproduced  until 
beyond  x/c  = 0.45,  at  wtiich  point  both  upper  and  lower  sm’face  agreement  is 
; excellent . 

e.  L51F07 

The  L5IFO7  configuration  is  shown  in  Figure  34'a.  This  consists  of 
a wing-body  combination  with  a MCA  65AOO6  airfoil  section  (Reference  50). 

The  wing  for  this  case  is  untwisted.  The  flight  conditions  for  this  case  are 
a Mach  number  of  0.93  aJid  an  angle -of -attack  of  2°.  Figure  36b  sliows  a plot 
of  the  upper  surface  isobars  predicted  by  this  program.  The  equivalent  3-b 
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upper  surface  pressure  distribution  is  shown  in  Figure  3<>c.  The  first  span 
station  shown  is  for  = 0.2  and  = 0.182,  Tlie  predictions  (Figure  3od) 

show  excellent  agreement  except  at  the  leading  edge  upper  surface  and  the 
trailing  edge.  The  leading  edge  can  be  accounted  for  by  the  fact  that  the 
span  station  used  for  the  predictions  is  inboard  of  the  span  station  used  for 
the  experimental  work.  The  trailing  edge  discrepancy  might  be  attributed  to  a 
viscous  corner  flow  displacement  effect  not  treated  in  the  present  program. 

On  conventional  airfoils,  viscous  effects  are  important  primarily  at  the  trail- 
ing edge  and  this  is  e\ idenced  by  the  excellent  agreement  between  the  data  and 
the  inviscid  solution  over  the  rest  of  the  section.  The  rest  of  the  wing 
(ilgures  36e,  36f,  3'^  and  S'^h)  show  similar  characteristics  with  the  exception 
that  as  one  moves  outboard,  the  agreement  forward  of  x/c  = .50  tends  to  decrease. 
This  could  possibly  be  improved  by  increasing  the  computational  angle -of -attack 
and  or  Mach  number  slightly. 

f.  A55B21 

The  A55B21  case  W!is  run  for  two  configurations  (Reference  57).  These 
configurations  can  be  seen  in  idgures  37a  and  37b.  As  can  be  seen,  the  wing 
planform  and  section  was  the  same  for  both  configurations.  How- 

ever, the  bodies  used  in  the  two  configurations  differed  considerably.  The 
airfoil  sections  originally  used  were  determined  by  using  the  cosine  of  the  mid- 
chord sweep  angle  since  that  was  the  specified  reference  sweep.  A second  set 
of  airfoil  ordinates  with  taper  affects  taken  into  account,  were  supplied  by 
David  Caughey.  These  ordinates,  yielding  slightly  better  results,  are  those 
for  which  the  results  are  presented.  The  fli^t  conditions  used  in  the  experi- 
mental wind  tunnel  work  were  for  a Mach  number  of  0.9^  and  an  angle-of-attack 
of  0°.  The  first  configuration  uses  the  basic  body  and  the  experimental 

results  were  obtained  at  71  = 0,471,  The  calculations  were  performed  at 

exp 

71  ^ = 0.473.  As  can  be  seen  (Figure  37c),  the  agreement  is,  in  general, 

prfe 

good  with  the  exception  that  the  shock  location  is  predicted  about  5^>  forward 

of  its  actual  location.  The  second  configuration  used  an  area  rule  body  with 

the  results  determined  at  the  same  span  locations  (Figure  37d).  Again,  the 

rapid  drop  in  pressure  is  predicted  between  5 and  10^  forward  of  the  actual 

location.  The  effect  of  area  rule  body  is  to  move  the  peak,  pressure  from 

x/c  = .75  forward  to  x/c  = O.5O.  Thus,  the  program  qualitatively 
pre  pre 
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reproduces  the  experimental  results. 

g.  Forward  Swept  Wing 

I 

The  forward  swept  wing  case  was  included  to  demonstrate  the  capahili-  j 

ties  of  the  program.  Although  no  experimental  results  exist  for  this  case,  i 

the  predictions  axe  nevertheless  Interesting.  Figure  38a  shows  the  forward  1 

swept  wing  configuraticn  (Reference  58).  The  leading  edge  sweep  of  this  con-  j 

figuration  is  -31.1°  and  the  trailing  edge  sweep  is  -42.97°.  The  wing  has  an 

aspect  ratio  of  5.22  and  a taper  ratio  of  0.4.  The  airfoil  section  used  was  one  of  jj 

the  Grumman  "K"  series  supercritical  airfoils.  The  flight  conditions  chosen 
for  this  case  were  of  a Mach  number  of  0.9  and  an  angle -of -attack  of 
8.5°.  Figure  38b  shows  the  predicted  3-D  pressure  distribution  on  the  upper 
surface  of  the  forward  swept  wing.  As  can  be  seen,  the  chordwise  pressure 
distribution  is  scMiewhat  flatter  than  that  predicted  in  most  of  the  other 
cases.  Figures  38c  and  38d  show  the  predicted  pressure  distributions  at 
T]  = 0.4  and  T]  = 0.8.  It  is  interesting  that  the  lower  surface  pressure 
distributions  are  much  flatter  than  those  usually  associated  with  supercriti- 
cal airfoils.  This  provides  an  example  of  one  of  the  non-standewd  configura- 
tions that  this  program  can  accept. 

h.  how  Aspect  Ratio 

The  next  case  calculated  consisted  of  a low  aspect  ratio  wing-body 
case  (Figure  39a).  Again,  no  data  exists  for  this  configuration.  The  config- 
uration consists  of  an  F-102  body  and  planform  (Reference  59^.  The  airfoil  section 
used  in  this  case  is  a NACA  65A00t)  section.  The  wing  planform  has  a leading 
edge  sweep  of  i'i0.1°  and  a trailing  edge  sweep  of  -5°.  The  aspect  ratio  of 
the  configuration  is  2.1  and  the  taper  ratio  is  .023.  The  flight  conditions 
used  for  this  case  are  a Mach  number  of  O.9O  and  iui  angle-of -attack  of  4°. 

Figure  39b  shows  the  3-D  upper  surface  pressure  distribution  for  this  config- 
uration. Figures  39c  and  39d  show  the  pressure  distributions  for  span  stations 
T)  = 0.4  and  T|  = 0.8.  A case  such  as  this  is  particularly  taxing  on  the 
internal  mesh  generators.  In  this  particular  case,  the  interior  mesli  had  to  be 
adjusted  beyond  the  tip  to  avoid  having  the  XI  = 0 and  XI  = 1 mesh  lines  cross 
beyond  the  tip.  This  leads  to  a large  amount  of  mesh  bending  at  tlie  wing  tip. 

It  is  apparent  that  some  care  must  be  exercised  in  setting  up  the  geanetry  of 
such  a configuration.  However,  this  case  demonstrates  that  the  progrojn  can 
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handle  such  a case. 


i . Transonic  Transport 

The  configuration  for  the  transonic  transport  is  in  figure  40a  (Reference  60) 
This  was  an  isolated  wing  calculation;  however,  the  experimental  results  are 
from  wind  tunnel  tests  of  a wing -body  configuration.  As  can  be  seen,  this 
configuration  is  an  example  of  a case  with  non-constant  leading  and  trailing 
edge  sweeps.  The  aspect  ratio  of  the  reference  trapezoidal  wing  for  this 
case  was  8.55.  The  airfoil  section  used  in  this  wing  is  a Grumman  designed 
supercritical  wing.  The  wing  was  designed  for  use  in  the  presence  of  winglets 
at  the  tip  and  an  engine  nacelle  over  the  trailing  edge  at  the  root.  The 
effect  of  the  nacelle  is  to  virtually  eliminate  the  shock  inboard.  The  effect 
of  the  winglets  is  to  straighten  the  isobars  near  the  wing  tip.  Consequently, 
the  results  shown  here  do  not  truly  represent  the  design  conditions  for  the 
wing.  The  Mach  nurabci  for  this  case  is  0.8l  and  the  angle -of-attack  is  2°. 

The  experimental  results  were  obtained  from  wind  tunnel  tests  at  an  angle-of- 
attack  of  There  was  some  uncertainty  as  to  the  wing  root  incidence  on 

the  model.  Since  a ccmputational  angle -of-attack  of  2°  seemed  to  give  better 
agreement  with  experimental  data,  this  angle  was  used  for  all  of  the  calcula- 
tions. It  sho’old  be  noted  that  the  airfoil  ordinates  used  included  a boundary 
layer.  Consequently,  only  an  inviscid  calculation  was  perfomed.  Figure  40b 
sliows  the  upper  siu*face  pressure  distribution  for  the  isolated  wing. 

figure  UOc  shows  the  pressure  distribution  at  = 0.l6  and  =.l45. 

The  agreement  for  this  station  is  not  very  good.  However,  this  is  to  be 
expected  because  the  computations  were  carried  out  on  an  isolated  wing,  wliile 
the  wind  tunnel  tests  included  the  body.  Figure  40d  begins  to  show  better 
agreement  as  the  effects  of  the  body  begin  to  decay.  At  = 0.5 

(T\^re  “ 0.509)  (Figure  40e),  the  agreement  gets  better  although  the  drop  in 
tlie  upper  stirface  press\xre  at  x/c  = 0.45  is  not  predicted.  At  = O.7O 

(Tlpre  " 0.727)  the  agreement  between  prediction  and  experiment  is  quite  good. 

The  effect  of  the  body  is  obviously  felt  quite  far  out  on  the  span.  This  case 
shows  that  pressure  distributions  on  wings  in  the  presence  of  bodies  cannot 
be  accurately  predicted  from  isolated  wing  cases  except  possibly  over  the 
outer  quarter  of  the  span.  Consequently,  it  is  important  to  include  the 
bo(.ly  in  the  geometry  model  and  the  more  accurate  the  boiiy  model  the  better 


103 


the  results  obtained  from  the  code  will  be. 


j . Advojiced  Figliter 

Hf^ui'e  4la  shows  a typical  advanced  fighter  concept  for  which  data 
is  available  Reference  61).  For  this  case,  we  see  that  a strong  shock  is 
predicted  at  the  wing  root  in  the  calculation,  as  shown  in  Figure  i+lb.  The 
experimental  data  shown  in  Figure  4lc  indicates  that  the  calculations  arc 
incorrect  at  the  root.  This  problem  continues  to  daninate  the  flowfield  until 
about  midspaji.  At  semispan,  the  predictions  begin  to  show  reasonable  a<^ree- 
ment  with  experiment.  The  striking  error  now  occurs  on  Uie  lower  surface, 
where  the  small  disturbcmce  theory  solution  caisistently  over-predicts  the 
higii  pi'essure  on  the  .lower  surface  unifomily  over  the  entire  chord.  In  addi- 
tion, the  comparison  suggests  that  the  calculation  is  being  made  at  a Macii 
number  tliat  is  too  low.  There  is  no  question  that  the  main  problem  in  this 
prediction  occurs  at  tiie  root  station.  One  source  of  the  error  in  this  case 
is  the  deciunbering  of  the  airfoil  section  in  the  glove  region,  which  could 
only  be  approximately  modeled  in  the  present  calcvilation.  However,  this 
simplification  does  not  explain  the  lower  surface  pressui’c  discrepiuicy , wliich 
occurs  consistently  across  the  span. 

k . Summary 

The  preceding  results  show  that  the  program  can  reliably  predict  the 
triuisunio  flow  over  a variety  of  configurations.  The  extreme  seiisitivity  of 
tnuasonic  flows  to  a precise  specification  of  freestrciuii  conditions  lias  been 
demonstrated  by  the  1\AF  "A"  calculatio.i.  For  wings  with  conventional  sections 
;uid  adv.anced  trjuisonic  tnuisports,  tlic  progrtuii  is  capable  of  making  good  pre- 
liictions  up  to  the  luigles-of-attack  at  wJiich  massive  shock  Induced  separation 
is  predicted.  The  advanced  fifihter  pi'edictions  require  furt.her  stiuiy.  In 
this  case,  the  flowfield  is  well  removed  from  a "small -disturb'uioe''  state  :ui  i 
the  wing-rv:)ot -glove  tiiselage  disturb'uices  arc  becoming  importiuit  ;uid  should 
be  investigated.  Nevertheless,  these  predictions  reflect  the  enonnous  advimces 
in  the  ability  to  canpute  three-dimensional  wing-body  viscous  tr-msi^mic  flows 
during  the  course  of  the  present  effort. 
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M = 0.840  RLF  = 3.000  CONV  = .327E-04  M = 0.840  HLF  = 3.000  CONV  = .327E-04 

LOCAL  CL=  0.284  LOCAL  CM=-0. 144  LOCAL  CD=0.0138  LXAL  a=  0.301  LOCAL  CM=-0.218  LOCAL  CD=0.0034 

SPAN  STATION  S r/S=0.182  44  MESH  PTS  Ofi  CHO  SPAN  STATION  13  T/50.436  44  MC5H  PTS  riN  OC 


M = 0.040  OLF  = 3.0IXJ  CONV  = . 327E-04  M = 0.040  PLF  = 3.000  CONV  = .327E-04 

LOCRL  a=  0.292  LOCAL  CM=-O.200  LOCfiL  C0=-.0032  LOCOL  CL=  0.258  LOCAL  CH=-0.312  LOCAL  CD=-.OO09 

5PfiH  STATION  19  f/S=0.655  44  MESH  PTS  ON  CHD  SPAf^  STATION  23  r/5=0.800  44  MESH  PTS  ON  CHD 


I = D.mo  flLF  = 3.000  CONV  = .327E-04  M = 0.840  FILF  = 3.000  CONV  = .327E-04 

•OCOL  CL=  0.200  LOCAL  CM=-0.28l  LOCAL  CD=-.0152  LOCAL  CL=  0. IC6  LOCAL  CM=-0.245  LOCAL  CD=-.0l89 

.PAN  STATIOtJ  26  r/5=0.g09  44  MESH  PTS  ON  CHD  SPAN  STATION  27  r/5=0.94S  44  MESH  PTS  ON  CHD 


. Hiwj  Dijur  Lijwr  louRHi  iUN  Figure  iic.  RFIE  WINt  coui  lunmujrhi  iun 

M = 0.900  flLF  = 1.000  CONV  = .991E-05  H = 0.900  OLF  = 1.000  CONV  = .991E-05 

LOCAL  CL=  0.075  LOCAL  CM=-0.035  LOCAL  COO.  0147  LOCAL  CL=  0.080  LOCAL  CM=-0.043  LOCAL  COO.  0100 

SPAN  STATION  6 r/S=0. 182  46  MESH  PT5  ON  CHO  SPAN  STATION  8 T/5=0.255  43  MESH  PT5  ON  CHO 


LOCRL  CL=  0.094  LOCRL  CM=-0.074  LOCAL  CO-.OOlO  LOCAL  CL=  0,111  LOCAL  CM=-0. 155  LXAL  C0=-.0072 

SPAN  STATION  12  T/S=0.400  44  MESH  PT5  ON  CHD  SPAN  STATION  18  T/S=0.618  44  ME5H  PT5  ON  CTC 
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M = 0.940  RLF  = 2.400  CONV  = .231E-02  N = 0.5«0  RLF  = 2.400  CtWV  = .231E-02 

LOCAL  CL=  0.120  LOCAL  CM=-0.257  LOCAL  CD=0.0100  LOCAL  a=  0.141  LOCAL  C«=-0.400  LOCAL  Cl>=O.OO0O 

SPAN  STATION  6 T/S=0. 192  48  MESH  PT5  ON  CHO  SPAN  STATION  l2  T/50.400  44  MESH  PTS  ON  DC 


LOCRL  CL=  0.168  LOCAL  CM=-0.685  LOCAL  CD=0.0058  LOCAL  CL=  0.182  LOCFfl.  CI1=-l.00H  LOCAL  CD=0.0CW1 

SPAN  STATION  18  T/S=0.618  iW  MESH  PT5  ON  CHD  SPAN  STATION  22  Y/S=0.764  “W  ME5M  PT5  ON  CHD 


LOCRL  CL=  0.153  LOCRL  CM=-0.086  LOCRL  CD=0.0153 
5PRN  5TRTI0N  6 T/5=0. 182  45  ME5H  PTS  ON  CHO 


LOCAL  CL=  0.164  LOCAL  CM=-0. 144  LOCAL  CD=0.004G  LOCAL  CL=  0.164  LOCAL  CM=-0.2ie  LOCAL  C0=-.0007 

5PAN  6TATI0N  12  r/6=0.400  44  ME5H  PTC  ON  CHD  CPAN  STATION  10  1/5=0.618  44  MESH  PT5  ON  CHD 
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(a;  Basic  model. 


lA/ing  geometry 

c/^ 

AR  = 6.0 


(b)  Area  rule  mo:iel. 


- IG.  37 

A55B^I  Configurations 


EXPERIPWTAL  DATA  . ^ EXPERIWITAL  DATA 


nc,  38  FoTWQTd  swept  wing 


I = 0.900  = 0.500  CONV  = ,9eiE-04  M = 0.900  RLF  = 8.500  CONV  = ,981E-04 

•XOL  CL=  0.915  LOCAL  CM=  0.135  LOCAL  00=0.0729  LOCAL  a=  0.870  LOCAL  CM=  1.200  LOCAL  C0=0. 1159 

;PAN  5TATI0rJ  12  r/5=0.‘400  ^C5H  PT5  ON  CHO  SPAN  STATION  23  T/5=0.800  MESH  PT5  CH) 


(^)  Upper  surface  pressure  disfributicri 


FIG.  JO  aspect  ratio 


M = 0.900  PLF  = 4.000  CONV  = .495E-02  M = 0.900  PLF  = 4.000  CONV  = .495F-02 

LOCAL  CL=  0.165  LOCAL  CM=-0.0g2  LOCAL  CD=0.0111  LOCAL  CL=  0.267  LOCAL  CM=-0.703  LOCAL  C0=0.0096 

5PAN  STATION  12  r/5=0.400  44  MESH  PTS  ON  CHO  SPAN  STATION  23  T/5=0.800  44  MESH  PT5  ON  CHD 


(b)  Upper  surface  pressure  distribution 
FIG,  41  Advanced  fighter 

m 


£ - 


L0CRL  CL=  0.753  lOCRL  CM=  0.954  LOCRL  CD=0. 1100  LOCAL  CL=  0.974  LOCAL  CM=  1.325  LXAL  CD=0.0884 

SPAN  5TBTI0r<  6 r/5=0.182  44  MESH  PTC  ON  CHO  5PRN  STATION  9 r/5=0.291  44  ME5H  PT5  ON  CHD 


F I SURE  41e,  OnVRNCED  FIGHTER  Figure  iiIf.  RDVPNCED  FIGHTER 

M = 0.^0  RLF  = 7.050  CONV  = .886E-01  M = 0.900  RLF  = 7.050  CONV  = .886E-01 

LOCRL  CL=  0.890  LOCRL  CM=  1.306  LOCAL  CD=0.0585  LOCRL  CL=  0.911  LOCRL  CM=  1.218  LOCRL  CD=0.0357 

5PflN  5TRTI0N  13  T/5=0.i*36  44  ME5H  PT5  ON  CHO  5PRN  STRTION  18  T/S=0.618  44  MESH  PT5  ON  CHO 


Figure  41s.  qcvnNCEL  FIGHTER  f^'suRE  41h,  PDVRNCED  FIGHTER 

H = 0.9Q0  RLF  = 7.050  CONV  = .8R6E-01  M = 0.900  RLF  = 7.050  CONV  = .886E-01 

LOCRL  CL=  0.910  LOCRL  CM=  1.075  LOCRL  CD=0.0206  LOCAL  CL=  0.782  LOCRL  CM=  0.766  LOCRL  CD=-.0012 

5PRN  STATION  22  r/5=0. 769  44  MESH  PTS  ON  CHO  SPAN  STATION  26  r/S=0.909  44  MESH  PT5  ON  CHD 


A.  DETAILED  DATA  COMPARISONS  - TACT  AND  F-8 


a.  TACT 

The  Transonic  Aircraft  Technology  (TACT)  program  provides  an  extensive 
set  of  wind  tunnel  and  flight  data  on  a modern  supercritical  wing  design  with  ^ 
which  the  computer  code  can  be  compared.  The  TACT  program  is  described  in 
Reference  62.  The  general  aircraft  layout  is  shown  in  Figure  42.  For  the 
transonic  flight  conditions  of  interest,  the  wing  sweep  was  26^.  A typical 
TACT  airfoil  section  is  shown  in  Figure  A3,  which  was  taken  from  Reference 
63.  The  glove  sections  were  assumed  to  be  the  same  as  the  extrapolated  wing 
sections  for  the  calculations  presented  here.  The  characteristic  wing  twist 
(which  is  contained  in  the  ordinate  definition)  is  shown  in  Figure  AA . Figure 

A 

A5  shows  the  aeroelastic  twist  increment  for  M = .91  (Reference  6A) . A 
similar  increment  was  applied  to  the  M = .86  case.  The  pressure  orifice  loca- 
tions at  ‘which  pressure  data  from  Reference  65  and  66  is  available  are  shown 
in  Figure  A6.  The  calculations  were  carried  out  for  this  case  using  the  MSD 
Theory  equation,  with  six  viscous-inviscid  iterations  required  to  obtain  a con- 
verged solution.  The  Reynolds  number  based  on  mean  aerodynamic  chord  was 
25.3  million,  which  corresponded  to  the  flight  test  condition.  The  case  was 
run  at  a second  angle-of-attack  to  demonstrate  the  sensitivity  of  the  calcula- 
tion to  angle-of-attack  change.  For  this  configuration,  NCR  differencing  was 
employed  due  to  the  FCR  prediction  of  the  shock  exactly  at  the  trailing  edge, 
which  made  the  initiation  of  the  viscous  interaction  difficult  due  to  the 
absence  of  an  adverse  pressure  gradient  at  the  trailing  edge. 

The  effect  of  the  aeroelastic  twist  increments  on  the  calculated  pressure 
distributions  can  be  seen  in  Figure  A7.  The  dashed  line  in  these  curves  repre- 
sents the  calculated  pressure  distributions  over  the  wing  with  only  the  built-in 
twist  modeled.  The  solid  lines  represent  the  calculated  pressure  distributions 
with  the  inclusion  of  the  aeroelastic  twist  increments  in  flight  article. 

Figures  A7a  through  A7b  show  the  effect,  at  a Mach  number  of  0.9  and  an  angle- 
of-attack  of  5.2°.  Figures  A7e  through  A7h  show  the  effect  at  a Mach  number  of 
0.85  and  an  angle-of-attack  of  7°.  The  effect  is  significant  and,  as  would  be 
expected  from  the  aeroelastic  twist  increments  shown  in  Figure  A5,  the  effect 


increases  at  the  outboard  span  stations. 


Figure  48  shows  the  co.nparison  ol’  predictions  ;uui  measurements  at  Much 
O.dl.  Tile  solid  lines  represent  the  calculations  at  :ui  angle-oC-attack  ol' 

;ind  tlic  daslicd  linos  represent  tin  ;uiglc-or-attacK  ol'  'I’iie  measure- 
ments were  made  at  an  angle -of -at tack  of  'liie  I'irst  sfvui  station  sliovni 

(Figure  48a)  is  at  T|  ■=  0.400.  As  can  be  seen,  tJie  leading  edgi'  pi'essiu-es  ai'e 
gener;U.ly  underpredictcil.  in  tiie  mid-chor<i  region,  tlie  lower  suid'ace  i.a  I'l'tt.er 
represented  by  a ^ An  increase  in  computationaJ  Maeli  nimiber  would  ]irob- 

ably  ePve  better  correlation  over  ttie  leailing  eilge.  Figure  481;  sliows  tlie 

correlation  at  11  0.570,  il  0.ol8.  Again,  tlie  uiipej'  surface  allows 

exp  pro 

bettor  correlation  at  o = o.P  , except  at  the  trailing  edg<’.  The  iowoj'  surJ'ai’c 

is  not  well  predicted  at  either  luigle-of-attack,  probably  due  to  the  fael.  tliat 

the  glove  and,  in  particular,  the  flowficld  around  the  engine  intiUscs  is  not 

moileled.  The  effect  of  the  flowl’ield  about  the  engine  intakes  eiui  lie  (.'xpeeted 

to  have  a large  effect  on  the  lower  surface  pressure  tiistribution.  I'iiaire  4i'e 

(n  0.'('(>1|,  ]]  O.'l'yG)  again  shows  better  correlation  on  the  upper  sair- 

Pfe  ^^exp  ^ 

face  at  a 0.2  . The  lower  surface  shows  better  correlation  at  a ''.2  . 


At  il  0.009  (il 

pre  ' 'exp 


0.921),  the  same  trcmls  hold  (Figure  48d).  iVcr  tJic 
entire  span,  the  upper  surface  generally  sliows  lietter  correlation  wlaai  the 
computvational  :uiglo-of-attack  agrees  wit.h  the  exjuvrimentai  ajig.i e-of-a l.tack , 
in  contrast  to  previous  efforts,  which  indicated  tiiat  the  lower  :uigie-ol'-at<.aek 
should  be  used . 

The  correlations  at  Mach  0.8o  are  significantly  in'tter  thaji  those  at.  Maeli 
0.91.  The  solid  line  I’epresents  calculations  at  an  ;uigle-id'-al,tack  of  V ;uid 
the  dashed  lines  represent  :ui  iingle  op  attack  of  8*^^  'I'lie  measurements  were  made 
at  un  angle  of  attack  of  8°  for  the  flight  tests  and  7-31^^  fev  the  wind  t.unnel 
tests.  The  upper  surl'ace  at  T]  = 0.400  (Figure  49a)  shows  consivlerabl e sent  l.er 
between  the  wind  tunnel  and  flight  test  results,  'llie  predietLons  again  show 
better  agreement  at  the  actual  angle  ol’  attack  on  the  upper  surface.  The  correla- 
tion here  is  much  better  with  flight  test  data  than  with  wind  tunnel  data.  I'he 
lower  surface  correlation  is  again  rather  disappointing,  for  the  same  reasons 
noted  above.  Figure  49b  again  shows  the  same  trends,  luid  in  fact,  these  trends  are 
duplicated  over  the  rest  of  the  span.  Figure  49c  ( ^ = 0.7b4,  O.'/’Mi)  shows 

excellent  agreement  as  to  tiie  location  and  strength  of  the  ndd-ctiord  sliock  at. 
or  8^  . At  a = 7°>  the  shock  is  predicted  too  far  forwaiai.  At,  tliis  span 
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station,  correlation  la  better  with  wind  tunnel  meaaureniento  than  with  rit(j;ht 
test  data.  Figure  4yd  also  shows  better  correlation  with  wind  tunnel  iiieasure- 
inents . 

It  is  apparent  that  upper  surface  agreement  is  better  when  the  coinput- 
tational  angle  of  attack  is  the  same  as  the  ex})er Imental  angle  of  attack.  The 
inboard  upper  surface  would  show  better  correlation  with  exiierlment  if  the  glove 
was  modeled  explicitly.  The  lower  surface  ae’reement  may  be  dominated  by  tiie 
effects  of  the  engine  intakes.  The  correlation  would  be  greatly  Imjiroved  if 
the  flow  field  around  the  intakes  was  modeled  by  modifying  tlie  boiiy  slopes  to 
Include  the  effect  of  this  flow  field.  In  the  absence  of  such  moilelitig,  the 
lower  surface  correlation  seems  to  be  improved  by  decreasing  the  computational 
angle  of  attack. 

Figure  ^Oa  shows  the  TACT  outer  panel  locating  t)ie  separation  regions 
occurring  on  the  panel.  On  the  upper  surface,  as  shown,  there  is  a shock 
induced  separation  region  near  the  trailing  ei.lge.  This  region  covers  i?iost  of 
the  mid-span  region.  On  the  lower  siu'face,  there  is  a large  sejiaratlon  region 
occurring  in  the  cove.  The  ext.ent  of  this  region  generally  ilecreases  on  the 
outboard  sections  of  the  span.  Figure  "jOb  shows  the  three  main  regions  requir- 
ing special  treatment  on  a super  critical  section  such  as  Is  present  on  TACT. 

The  effect  of  adding  a boundary  liiyer  on  tlie  solution,  can  be  seen  In  Figure  SOc. 
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F/G.  4-7cl  Aeroelastic  effects 


-1.6 


Isolated  wing  code Isolated  wing  code  -''9  s/iope 

Medium  grid  Medium  grid 

FIG  47e  Aeroeiast,c  effects  pfQ  Aetoelastic  effects 


LOCAL  CL=  0.669  LOCffl.  CM=-2.7i49  LOCAL  CD=0.0436  LOCAL  CL=  0.606  LOCAL  CM=-2.934  LXRL  C0=0.0152 

5PAN  STATION  12  r/5=0.400  44  MESH  PTS  ON  CHD  SPAN  STATION  18  T/5=0.61B  44  MESH  PTS  ON  CHD 


M = 0.860  RLF  = 7.000  CONV  = .2996^29  M = 0.860  RLF  = 7.000  CONV  = .299E’-29 

LOCAL  CL=  0.785  LOCRL  CM=-3. lOl  LOCAL  00=0.0429  LOCAL  CL=  0.768  LOCAL  CM=-3.503  LOCAL  00=0.0294 

5PRN  station  12  r/5=0.400  44  MESH  FT5  ON  CHO  SPAN  STATION  17  T/S=0.582  44  MESH  PTS  ON  CHD 


TACT  tact 

M = 0.860  ALF  = 7.000  CONV  = M = 0.860  ALF  = 7.000  CONV  = .299E*-29 

LOCAL  CL=  0.691  LOCAL  CM=-3.694  LOCAL  CD=0.0144  LOCAL  CL=  0.534  LOCAL  CM=-3.281  LOCAL  CD=0.0008 

5PAN  5TATI0N  22  7/5=0.764  44  MESH  PT5  ON  CHO  SPAN  ^TAFION  26  7/5=0.909  44  MESH  PTS  ON  CHO 


(a)  TACT  OUTER  PANEL 


(b) 


b.  F-8 


The  NASA  F-8  Supercritical  Wing  Program  provides  an  excellent 
opportunity  for  comparing  the  computer  program  with  an  extensive  set  of 
recently  declassified  wind  tunnel  and  flight  data.  The  details  of  the  experi- 
mental program  are  surveyed  in  Reference  6?.  The  F-8  supercritical  wing  was 
designed  for  a cruise  lift  coefficient  of  .4  at  a Mach  n\miber  of  .99-  After 
further  testing,  area  rule  fuselage  modifications  were  developed  with  the 
intention  of  producing  a weaker  shock  pattern  at  speeds  approaching  Mach  1.0. 

In  addition,  vortex  generators  were  added  to  the  lower  surface  leading  edge 
of  each  wing.  These  vortex  generators  were  designed  to  alleviate  an  unstable 
pitching  moment  problem  that  occurred  at  high  lift  coefficients.  Both  the 
model  and  flight  test  vehicle  were  tested  with  area  rule  body  additions  and 
vortex  generators  on  and  off.  The  configuration  for  which  predictions  were 
made  and  results  were  compared  had  the  area  rule  body  additions  on  and  the 
vortex  generators  off. 

The  basic  configuration  tested  is  shown  in  Figure  51.  Note  that  the 
wing  is  mounted  at  the  top  of  the  body  with  a glove  extending  far  forward 
over  the  fuselage.  Tlie  amo\ant  of  blending  between  the  wing  and  the  body  is 
extensive.  This  blending  was  not  specifically  modeled  in  the  geometry  input 
to  the  program.  Instead,  the  body  geometry  computational  input  was  based  on 
the  area  distribution  given  in  Reference  68,  which  included  the  area  rule 
fuselage  additions.  Figure  52  shows  a typical  streamwise  section  of  the 
supercritical  airfoil  used  in  the  F-8  program.  The  wing  twist  present  in 
the  model,  together  with  the  loaded  and  unloaded  flight  test  wing  twist,  is 
shown  in  Figure  53-  The  twist  distribution  input  into  the  program  is  that 
shown  for  the  model  with  the  model  airfoil  sections  employed.  The  experimental 
pressure  distributions  were  measured  at  the  locations  shown  in  Figure  54. 

The  experimental  data,  used  for  the  comparisons  to  follow,  was  obtained  from 
References  68-70. 

The  calculations  were  carried  out  with  the  viscous  interaction  option 
for  a Reynolds  nunber  based  on  the  mean  aerodynamic  chord  of  10  million  for 
the  M = .99  case  and  9 million  for  the  M = .90  case,  which  corresponds  to 
the  flight  condition.  Transition  was  specified  at  5 percent  chord.  Start- 
ing from  a saved  inviscid  solution,  with  4 inviscid-viscous  iterations  required 
in  order  to  converge,  such  that  the  wing  lift  coefficient  changed  by  less 


( 


I 


154 


than  .002  betvreen  iterations.  Figure  55a  shows  the  experimental  spanloads 
obtained  fron  flight  and  wind  tunnel  test  at  a Mach  number  of  approximately 
0.90.  From  this,  the  cases  chosen  for  comparison  were  the  flight  test  at 
an  angle-of-attack  of  ^+.03°  and  the  wind  tunnel  at  an  angle -of -attack  of 
3.56°.  The  computational  angle-of-attack  was  3.56°.  The  calculated  span- 
load appears  as  the  solid  line  in  the  figure.  Figure  55b  shows  the  measured 
spanloads  at  M = 0.99*  T'he  cases  chosen  for  comparison  at  M = O.99  are  a 
wind  tunnel  angle-of-attack  of  3*^3°  and  a model  angle-of-attack  of  3.33°. 

The  calculated  spanload  is  shown  as  the  solid  line  in  the  figure.  Calcula- 
tions are  compared  for  both  3*33  and  4.33  degrees  ajigle-of-attack. 

Figure  56a  begins  the  detailed  pressure  comparisons  as  can  be  seen,  the 

computational  span  station  was  at  = .327.  The  experimental  span  station 

was  at  T|  = O.306.  As  noted  in  some  of  the  other  cases,  this  close  to  the 
exp 

body,  where  the  spanwise  gradients  are  relatively  high,  the  agreement  would 
probably  be  better  if  the  experimental  and  computational  span  stations  were 
closer.  The  upper  surface  is  not  well  predicted,  possibly  due  to  the  rela- 
tively simple  body  model  used,  without  the  wing/ fuselage  blending  and  body 
carry-over  effects.  With  a more  detailed  body  input,  the  results  could  be 
expected  to  improve  at  the  inboard  wing  sections.  Figure  5bb  = 0.473, 

^exp  ~ 0.480)  shows  better  agreement  near  the  trailing  edge,  especially  with 
the  wind  t\innel  data.  The  forward  part  of  the  upper  sm-face  is  stil]  not 
well  predicted,  although  the  lower  surface  leading  edge  pressure  distributioi. 
is  well  predicted.  It  can  be  seen  in  Figure  56c,  that  the  agreement  improves 
moving  outboard  on  the  span.  This  span  station  = 0.655,  = O.o53) 

shows  better  agreement  over  the  entire  upper  surface.  hVen  the  leading  eilge 
peak  is  relatively  well  predicted.  The  next  spiui  station  = O.8OO, 

^exp  ~ 0.8o4)  again  shows  better  agreement  on  the  upper  svirface  (Figui'e  5('d). 

The  final  span  station  (Tl  = 0.945.  Tl  = .933)  shows  decent  overall 

pro  exi’ 

agreement  with  experimental  data  (FijTure  5l>e)  altliougii  as  in  the  previous 
span  station,  tiie  predicted  lower  surface  Ci)'s  are  somewiiat  lower  tlnin  the 
measured  values.  It  appears  tliat  some  adjustment  of  computational  angle-of- 
attack  and/or  Mach  number  would  yield  better  results,  although  the  primary 
factor  would  be  a better  body  model. 

The  Mach  = 0.90  comparisons  begin  with  Figiu’e  57a  = 0.f':'>l,  - 

0.307).  The  solid  line  in  tlie  figure  is  tlie  prediction  at  rvn  !uigle-of-attack 


o 

of  3 >33  j the  dashed  line  represents  predictions  at  an  angle-of -attack  of 

4.33°.  As  can  be  seen,  a = U.33°  tends  to  give  a more  realistic  leading 

edge  peak  although  again,  the  body  effect  may  be  incorrectly  incorporated 

and  the  overall  agreement  is  somewliat  disappointing.  Over  the  rest  of  the 

section,  a = 3*33°  yields  somewhat  better  results.  Figure  57b  shows  the 

results  at  Tl  = 0.473,  ^ = 0.480.  An  angle-of -attack  of  4.33°  again 

pre  exp  ^ 

yields  better  results  at  the  leading  edge.  The  scatter  in  the  experimental 
data  over  the  rest  of  the  section  makes  comparison  with  prediction  somewhat 
difficult.  However,  it  appears  that  the  lower  surface  leading  edge  flight 
test  data  is  better  predicted  at  o = 3-33°,  wliile  the  wind  tunnel  data  is 
better  predicted  at  a = 4.33°-  In  Figure  57c  = 0.655,  = 0.(.)53), 

it  can  be  seen  that  the  entire  upper  surface  is  better  predicted  at  a = 

4.33°*  The  lower  surface  is  not  well  predicted  at  either  suigle-of-attack . 
Again,  there  is  significant  scatter  in  the  experimental  data.  The  next 
span  station  (4p^g  = O.8OO,  11^^^  = 0.804)  again  shows  better  upper  surface 
agreement  at  a = 4.33°,  nnd  again,  the  lower  surface  data  is  reproduced  at 
neither  angle -of -attack  (Figure  57d).  It  appears  that  an  increase  in  the 
computational  Mach  number  would  lead  to  improved  correlation.  The  final 
span  station  (Tp^g  = 0.945,  ^gj^  = 0.933)  shows  better  agreement  with  flight 
test  at  o = 3.33°  a^d  better  agreement  with  wind  tunnel  data  at  o = 4.33° 
on  the  upper  surface  (Figure  57e).  Again,  the  lower  surface  is  not  well 
predicted  at  either  angle-of-attack  and  an  increase  in  computational  Mach 
number  would  improve  the  correlation. 

The  overall  correlation  between  the  computational  predictions  and  the 
flight  and  wind  turuiel  experimental  results  is  perhaps  disappointing.  However, 
the  highly  blended  fighter  type  configurations  present  the  most  challenging 
cases  for  the  present  prograiu.  The  complex  contributions  of  boUi  body  ;uui 
viscous  effects,  together  with  the  potentially  large  effect  of  acroelastic 
corrections  provide  several  sources  to  investigate  for  the  disagreement 
between  computation  and  experiment.  Supercritical  flows  are  relatively 
sensitive  to  the  various  effects  so  that  configurations  such  as  the  F-8 
provide  a good  baseline  for  ITirther  investigation  of  the  computational 
methods. 
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Figure  53.  Wing  Twist 
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Figure  55b.  Tunnel  and  Flight  Test  Spanload  for  NASA  F-8  Supercritical  Wing 
(continued) 
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5.  COMMON  INPUT/ OUTPUT  FEASIBILITY 


a.  Introduction 

Typical  aerodynamic  analysis  and  design  assignments  require  the 
aerodynamicist  to  draw  on  a variety  of  previously  obtained  aerodynamic 
data  and  ccmputational  aids,  as  well  as  experience  and  judgement.  As 
computer  capabilities  have  expanded,  correspondingly  elaborate  aerodynamic 
flowfield  prediction  computer  codes  have  been  developed.  These  codes  can 
now  provide  reliable  predictions  for  some  fairly  realistic  configurations 
and  flight  conditions.  These  advances  have  led  to  a greatly  increased 
reliance  on  computer  codes  to  provide  aerodynamic  predictions.  The 
increased  dependence  on  the  computer  leads  to  some  fund:imental  changes 
in  the  approach  to  large  engineering  jobs  in  general,  as  well  as  to  the 
specific  aerodynamic  jobs.  In  this  feasibility  study,  we  make  a brief 
survey  of  the  modern  engineering  environment  and  identify  the  procedures 
required  to  maximize  the  use  of  the  available  computing  resources. 

Clearly,  the  proper  interfacing  of  the  various  computer  codes  via  the 
input/ output  data  streams  is  one  of  the  crucial  steps  in  efficient  use 
of  computer  systems. 

b.  Common  l/O  - A Systems  Viewpoint 

Large  government  procurement  programs  involving  the  developnent 
of  technologically  advanced  systems  require  precise  planning  and  in:uiage- 
ment  in  order  to  proceed  on  schedule  and  within  budgetary  constraints. 

The  systems  approach  to  the  management  of  the  engineering  tasks  is 
usually  termed  an  "integrated  design  procedure.”  Gru/nman's  integrated 
design  procedure  is  RAVES  - Rapid  Aerospace  Vehicles  Evaluation  System 
(Reference  71).  RAVES  evolved  from  a previous  system  used  at  Grumman 
to  manage  the  F-l4  structural  design. 

In  concept,  RAVES  includes  all  major  aerospace  vehicle  analyses 
which  are  formiilated  via  computer  programs.  The  main  function  of  RAVES 
is  to  develop  a systems  framework  or  architecture  into  vriiich  computer 
programs  from  all  technical  disciplines  can  be  integrated,  resulting  in 
a data  management  system  that  can  dynamically  pass  data  from  one  program 
to  another.  RAVES  can  function  for  the  full  scope  of  the  engineering 
design  effort,  from  preliminary  design  through  detail  design.  The  system 
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uses  a time-share  computer  vfith  remote  job  entry  links  to  the  batch 
computer.  The  time-share  computer  is  used  for  user/computer  hard'ware 
communications,  to  assemble  data,  and  to  run  the  smaller  computer 
programs;  the  batch  machine  is  used  for  the  execution  of  the  larger 
computer  programs.  Users  run  programs  from  remote  terminals  equipped  with 
scope  and  digitizer  facilities. 

By  making  block  diagrams  of  the  flow  of  information  required  for  a 
design  effort,  the  input  and  output  requirements  of  individual  computer 
analysis  can  be  integrated  into  the  entire  system.  This  approach  then 
leads  to  a set  of  consistent  I/O  specifications  for  all  computer  codes 
in  the  RAVES  system.  The  resiilting  system  ensures  a smooth  flow  of  data 
from  one  program  to  the  next,  with  a minimum  of  routine  data  manipulation 
at  the  interface  between  analysis  progrtuns.  Figures  58  and  50  contain 
samples  of  the  block  diagrams  which  can  be  used  to  manage  the  data  flow. 

Aerodynamic  routines  contained  in  the  block  diagrajiis  are  labeled 
as  AXX,  wtiere  XX  refers  to  a specific  aerodynamic  program.  t)nc  aero- 
dynamic program  is  actually  used  solely  for  the  conversion  of  tiie  pre- 
liminaj-y  3-view  drawing  to  an  analytic  description  of  tlie  surface 
definition.  This  program  is  shown  in  Figure  58  and  is  desitaiated 
A9-QUICK.  Figure  ‘,i0  indicates  the  locaticn  of  the  panel  method  wing/body 
codes  in  the  design  environment,  'llie  present  code  will  supplement  these 
codes  at  the  same  location  in  tlie  diag2*am.  The  need  for  cc«nion  input/ 
output  is  evident  from  this  diagram,  since  all  tliese  cixles  occure  in  the 
same  location  in  the  functional  diagrajn  tuid  a typical  design  effort  wili 
require  predictions  fi'om  eacli  of  these  codes. 

c.  Canmon  I/O  - The  Aerodynamicist 

Aerodynamic  analysis  requires  the  ability  to  operate  a large 
number  of  different  prograjns.  Thus,  in  addition  to  a knowledge  of  aero- 
dynamics, the  job  requires  familiarity  with  the  input  peculiarities  of 
each  program  as  well  as  the  particular  notations  used  in  the  output  of 
the  solution  of  each  prograjn.  Differences  between  these  programs  often 
lead  to  painful  time  delays  iuid  confusion.  It  is,  of  course,  tlie  immense 
range  of  powerful  aerodynajiiic  simulations  available  to  the  aerodynamicist 
via  the  canputer  that  brings  about  these  problems.  Without  careful 
planning,  we  in  fact,  have  an  embarrassment  of  riciies.  Unless  Uiis 
situation  is  brought  under  control  now,  the  future  promises  to  bring 


even  more  confusion.  More  and  more  codes  are  becoming  available,  each 
with  its  unique  data  structure  and  output  format.  The  final  result  could 
be  that  the  improving  computational  capabilities  will  not  be  fully 
reflected  in  actual  engineering  practice  due  to  common  l/O  limitations. 

d.  Common  l/O  - NASF 

NASF,  the  Numerical  Aerodynamic  Simulation  Facility,  is  a NASA 
proposed  centralized  facility  that  would  be  made  available  to  the  aerospace 
community  in  much  the  same  manner  that  NASA  Wind  Tunnels  are  made  available 
(Reference  72).  The  facility  would  house  a new  special  purpose  computer 
designed  specifically  for  the  prediction  of  aerodynamic  flowfields  and  a 
hierarchy  of  computer  codes  that  would  be  dedicated  to  this  particular 
machine,  each  code  providing  the  flowfield  simulation  at  a different  level 
of  sophistication.  Any  internal  or  external  common  l/O  efforts  must  be 
closely  coordinated  with  this  effort,  since  the  coordination  of  the  data 
: set  for  NASF  with  in-house  aerodynamic  tools  will  be  by  far  the  most 

1 important  requirement  for  effective  use  of  this  facility. 

e.  QUICK  - A Baseline  Code  For  Generating  Bodj"^  Definition 

The  tool  for  generating  the  basic  body  definition  for  aerodynamic 
simulation  is  presently  available.  Designated  A9,  QUICK  is  already  part 
of  the  Grumman  RAVES  system.  Reference  73  describes  this  code.  Briefly, 
QUICK  provides  a simple  method  of  rapidly  defining  an  accurate  computer 
model  of  complex  vehicle  geometry.  The  pritiuiry  motivation  which  guided 
the  development  of  QUICK-GEOMETRY  was  the  elimination  of  the  redundant 

!i  task  of  remodeling  geometry  input  to  run  different  programs.  With  this 

I 

I system,  only  a single  effort  is  made  to  extract  a mathematical  model  from 

I the  design  drawing.  This  math  model  can  then  be  manipulated  by  computer 

I 

I subroutines  to  generate  the  required  geometry  input  data  decks.  Tlie  math 

i model  is  quite  general  and  meets  all  the  fundamental  requirements  of  a 

I 

1 geometry  model: 

(1)  It  is  independent  of  the  reference  source  of  configuration 
geometry,  hence  the  model  can  be  developed  with  equal  facility  from  a 
conceptual  description,  a lines  drawing,  or  a configuration  as  stored  in 

, the  Grumman  Master  Dimensions  System, 

(2)  It  is  adaptable  to  both  manual  input  using  ruler  and  template 

I to  specify  local  geometry,  and  to  machine  assisted  digitized  input, 

which  can  measui’e  selected  points  or  trace  out  a series  of  points  along 

fO 


an  arc. 


(i)  Ttie  model  is  sufficiently  flexible  to  adapt  to  various  types 
of  analysis.  It  will  develop  comparable  models  for  different  analyses  and 
comparable  models  for  different  configurations.  It  will  easily  adapt 
models  to  various  restrictions  imposed  by  analysis  limitations. 

(4)  Subsequent  model  refinement,  when  necessai-y  to  support  a 
more  detailed  analysis,  or  to  keep  pace  with  actual  configuration  changes 
can  be  rapidly  made  to  the  initial  quick  response  model. 

(5)  Uses  a "building  block"  approach  which  defines  each  config- 
uration component  in  its  own  natural  coordinate  system  to  simplify  its 
geometric  definition.  This  also  makes  it  possible  to  modify'  existing 
models  by  scaling,  translating  or  rotating  components  and  to  select  com- 
ponents from  a catalog  of  standard  ones. 

(6)  The  model  is  very  concise  and  can,  therefore,  be  stored 
efficiently  in  a computer  data  bank  (such  as  RAVES). 

The  QUICK-GEOMETRY  system  is  being  applied  in  a number  of  ways.  As 
part  of  RAVES,  it  is  used  as  a common  source  of  vehicle  geometry  which  is 
drawn  upon  to  generate  a variety  of  geometry  input  data  decks.  In  particular, 
data  decks  conforming  to  the  Harris  input  format  (Reference  7^),  wliich  has 
become  a standard  at  MSA  Langley,  can  be  generated.  In  another  RAVES 
application,  QUICK-GEOMETRY  has  been  adapted  to  a vehicle  lofting  system. 

The  most  sophisticated  application  to  date , has  been  to  integrate  the  QUICK- 
GEOMETRY  system  with  a numerical  flow  code  which  calculates  the  steady/ 
supersonic/hypersonic  inviscid  flow  around  real  configurations  (Reference 
75 )•  In  this  application,  vehicle  geometry  (surface  coordinates,  slopes 
and  normals)  are  generated  as  required  by  the  flow  analysis  code. 

Figure  60  shows  some  examples  of  the  typical  model  produced  by  QUICK. 
Because  the  code  described  in  this  report  contains  a very  basic  body' 
simulation,  the  QUICK  code  has  not  been  incorporated.  This  allows  the 
present  code  to  retain  its  stand  alone  configuration.  The  body  modeling 
used  in  the  present  code  has,  however,  been  strongly  influenced  by  QUICK. 

It  appears  that  the  proper  use  of  QUICK  would  be  to  continue  to  operate 
QUICK  as  a stand  alone  format,  with  capability  of  simultaneously  generating 
data  sets  for  the  various  aerodynamic  programs. 

Although  QUICK-GEOMETRY  modeling  might  not  be  required  for  all  appli- 
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cations,  this  body  oriented  system  (which  is  generally  available  through 
NASA)  would  allow  for  a universal  input  set  that  would  meet  the  require- 
ments of  all  the  various  aerodynamic  programs,  Grumman  is  presently 
considering  constructing  a processor  to  automatically  generate  body  panel 
models  for  sub-  and  supersonic  linearized  codes.  We  can  conclude  that 
common  body  input  is  feasible  and  the  basic  elements  already  exist  for 
common  subsonic,  transonic  and  supersonic  code  input.  The  main  considera- 
tion in  the  construction  of  common  input  is  the  particular  system  of  codes 
and  the  particular  time-sharing  system. 

f.  Wing  Definition 

The  planform/section  modeling  problem  is  less  complicated  than  the 
body  case.  The  representation  of  the  planform  by  (y,  x)  pairs  is  relatively 
simple  and  very  good  accuracy  can  be  achieved.  Because  QUICK  is  oriented 
toward  body  cross-sections  (which  are  natural  for  bodies),  it  is  not 
appropriate  for  most  wing  section  definitions,  although  low  aspect  I’atio 
wings  are  adequately  treated  by  QUICK.  The  wing  thickness  should  be 
defined  by  inputting  streamwise  sections.  The  method  of  input  employed  in 
the  present  code  should  form  a baseline  for  l/O  standardization.  In  fact, 
by  using  the  thickness  scaling  option,  it  should  be  possible  to  use  the 
same  airfoil  ordinate  deck  for  both  this  code  and  the  various  2-D  airfoil 
analysis  codes  in  most  cases. 

g.  Feasibility  of  Including  Graphical  Results 

Common  output  should  include  a uniform  set  of  plots  and  a uniform 
nomenclature.  This  allows  for  verification  of  the  geometry  and  should  also 
be  sufficient  for  report  and  presentation  figures.  The  graphics  package 
included  with  the  present  program  provides  wing/body  geometry  verification 
level  consistent  with  the  actual  geometric  treatment. 

h . Common  Output 

We  also  endorse  a uniform  output  of  results  for  various  codes. 

At  present  there  is  usually  a significant  confusion  factor  concerning  just 
exactly  what  each  different  code  outputs.  Subtle  changes  in  the  various 
definitions  of  output  quantities  between  various  codes  is  a continued 
source  of  frustration.  In  the  present  program,  the  output  will  be 
presented  in  as  clear  and  definitive  manner  as  possible.  In  fact,  the 
output  appears  to  be  one  of  the  best  formats  of  any  of  the  contemporary 
codes . 


17. ‘ 


i.  Recommendations  for  Common  l/O 

(1)  The  field  of  computational  aerodynamics  has  reached  a 
sufficiently  high  maturity  level  that  commonality  between  code  input/ 
output  must  be  a consideration  in  code  development. 

(2)  Practical  operation  of  the  EiSF  concept  will  require  a 
standardization  of  l/O. 

(3)  Grumman's  QUICK  Geometry  Modeling  should  form  the  basis  for 
body  and  low  aspect  ratio  wing  modeling. 

(4)  The  present  code  l/o  has  been  developed  with  an  awareness 
of  common  l/O  considerations,  so  that  it  could  easily  be  modified  to 
conform  to  any  standardization  of  l/O  for  aerodynamic  simulation  --  it 
is  feasible  to  couple  this  code  to  QUICK. 
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flGURE  56  T^.'plcal  Poir.t  Design  Data  Flow 


59  Typical  Design  Data  Flow  Req,u1reK.ents 


6.  DESIGN  FEASIBILITY 


a.  Introduction 

Although  the  use  of  the  computer  to  simulate  the  flowfield  about  a 
vehicle  vdth  a specific  geometric  configuration  is  an  extremely  useful  and 
important  capability,  it  is  an  indirect  response  to  the  aerodynamic  design 
question.  The  aerodynamic  design  question  is,  of  course,  typically  posed  at 
several  levels,  starting  with  some  vague  and  general  question  about  the  "best" 
shape  of  the  airplane  for  a particular  mission,  and  proceeds  to  more  specific 
and  detailed  questions  concerning  the  actual  wing  lines,  subject  to  a large 
variety  of  constraints.  Often,  a wing  is  designed  using  computer  programs 
submitted  and  wind  ttinnel  tested.  In  the  analysis  mode,  the  aerodynamic  com- 
puter programs  are  being  used  to  simulate  a wind  tunnel.  Of  course,  the  com- 
puter simulation  can  be  used  much  sooner  in  the  design  cycle  than  a wind 
tionnel  test  and  this  strategy  should  produce  an  improved  final  design  at  a 
reduced  cost,  in  a shorter  time  period.  This  is,  of  course,  the  proper  initial 
introduction  of  the  computer  simulations  into  the  wing  design  process.  Indeed, 
this  technique  has  been  adopted  for  subsonic  and  supersonic  wing  design  for 
some  time.  Transonic  wing  design  efforts  have  only  recently  started  using 
fully  transonic  three-dimensional  wing-body  computer  simulations  in  the 
design  cycle. 

Once  the  computer  is  introduced  into  the  design  cycle,  it  becomes  evident 
that  it  can  be  used  in  a fundamentally  different  mode  than  to  simply  supple- 
ment wind  tunnel  testing.  The  use  of  flowfield  simulation  in  this  manner  is 
naturally  referred  to  as  the  "design  mode,"  as  opposed  to  the  "analysis  mode” 
of  operation.  A "design  mode"  has  been  available  for  linearized  subsonic  and 
supersonic  flowfields  since  shortly  after  the  analysis  codes  became  available. 
The  most  extensive  use  of  a "design  mode"  appears  to  have  been  the  elaborate 
system  of  supersonic  wing  design  codes  that  evolved  from  the  work  of  Caxlson 
and  Middleton  (Reference  76).  After  a brief  review  of  the  design  problem  and 
some  of  the  presently  employed  tools,  we  attempt  in  this  study  to  identify  a 
reasonable  approach  to  the  use  of  the  present  three-dimensional  trajisonic 
methodologyr  in  a design  mode,  wtiich  is  then  illustrated  with  some  examples. 


b.  Review  of  the  Design  Process 


A variety  of  possibilities  emerge  wtten  the  problem  fomiuiation  for  a 
"design  mode"  of  operation  is  explored.  The  reason  for  this  range  of  possi- 
bilities can  be  attributed  to  the  manner  in  which  the  design  problem  is  posed, 
as  noted  above.  Ideally,  the  aircraft  designer  would  specify  the  aircraft 
mission  (or  missions)  and  a computer  program  would  provide  the  detailed  lines 
of  the  optim’ura  aircraft.  Of  course,  we  do  not  expect  that  such  a smart  com- 
puter program  will  exist  for  some  time.  However,  most  aircraft  companies  and 
governmental  agencies  routinely  employ  programs  that  predict  the  gross  fea- 
tures of  an  optimum  aircraft  for  a particular  mission  with  some  assumption 
regarding  the  rate  of  development  of  various  technologies.  These  programs 
also  include  a large  data  base  developed  from  previous  aircraft  designs.  Typi- 
cal aerodynamic  outputs  from  the  programs  are  i)  Aspect  Ratio,  AR;  ii)  Taper 
Ratio,  X;  iii)  Sweepback,  A;  and  iv)  Thickness  Ratio,  t/c.  Usually  a target 
drag  level  for  the  configuration  is  also  specified.  Typical  examples  of  this 
type  of  program  are  the  Grumman  CISE  program  (Reference  77 ) and  the  Ames 
ACSYNT  program  (Reference  78). 

Hence,  the  canputer  is  used  to  determine  the  overall  features  of  the 
required  airplane.  The  typical  aerodynamic  design  problem  thus  becomes  less 
vague  and  more  manageable,  with  the  statement  being  reduced  to  something  along 
the  following  line: 


Given : 


o AR,  X,  A,  t/c 
o M , R 

.•v>  o 


CRUISE 


MAX  ALLOWABLE 


(basic  geometric  requirement') 
(flight  regime) 


Design  Goal 


Find:*  o C for  C 
°M1N  n 


'CRUISE 


for  C.. 


-MAX  ALLfWABLE 


Design  Goal 


* Note  that  although  this  is  Ahe  typical  protlom  statement,  each  ir  actually 
a requirement  to  obtain  L In  one  case.  1.  is  .si'oci f iel , in  the 
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o Detailed  Geometry  ->  Detailed  Aerodynamic  Aircraft  Design, 
Definition, 

subject  to  geometric  constraints  on  twist,  camber,  root  bending  moment,  etc,, 
and  aerodynamic  requirements  on  performance  at  other  flight  conditions.  At 
this  point,  we  coiild  begin  to  consider  the  use  of  a computer  code  directly  to 
help  determine  the  optimum  aerodynamic  shape  and  performance  that  can  be 
obtained  for  the  specified  problem.  More  typically,  the  designer  employs  his 
experience  and  judgement  to  specify  a desired  pressure  distribution.  This 
type  of  program  is  usually  described  as  an  "inverse  method,"  while  a program 
that  attempts  to  address  the  problem  more  immediately  is  usually  termed  an 
"optimization  method,"  Each  of  these  approaches  has  its  own  strengths  and 
weaknesses . 

A contrast  between  optimization  and  inverse  methods  can  be  summarized 
as  follows: 

Optimization 

o Requires  many  analysis 
submissions  for  a single 
design  case, 

o Solution  depends  critically 
on  the  user  assmed  form  of 
tlie  answer. 

o Can  handle  a variety  of 
geometric  and  off  design 
constraints. 

o If  performed  through  large 
optimization  code,  solution 
is  not  obtained  from  "aero- 
dynamic thinking." 

Another  drawback  of  the  optimization  approach  is  that  the  path  taken  to 
the  final  result  is  often  rather  obscure  and  the  relative  importance  of  the 
various  aspects  of  the  final  design  produced  in  this  manner  are  not  readily 
apparent . 

I>ir  orM’T' ; .-at  liir.  ‘eotiri^ues  presently  rmrli'yed  in  Un-  lest.-r  -.••.hi'i: 
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Inverse 

o Generally  as  fast  as  a single 
analysis . 

o The  geometry  may  not  always  exist 
for  a given  pressure  distribution. 

o Difficult  to  treat  off  design  and 
geometric  constraints. 

o Solution  as  a direct  resvilt  of  best 
current  "aerodynamic  thinking." 


vaxiations  in  order  to  obtain  the  maxima.  Although  an  entire  book  (Reference 
79)  has  been  devoted  to  aerodynamic  optimization  using  calculus  of  variations, 
these  concepts  have  not  yet  proven  to  be  valuable  for  the  class  of  optimiza- 
tion problems  under  study  here.  It  would  not,  however,  be  wise  to  automati- 
cally dismiss  these  methods  from  consideration  in  general,  althoiigh  we  will 
not  discuss  these  methods  further  in  the  present  study. 

c.  Review  of  2-D  Transonic  Design  Methods 

A variety  of  numerical  approaches  have  been  used  to  design  transonic 
airfoils.  The  book  by  Thwaites  (Reference  80)  discusses  the  classical  approaches 
to  the  incompressible  inverse  methods  and  points  out  that  some  judgement  must 
be  used  by  the  designer  in  specifying  the  desired  pressure  distribution;  a 
solution  does  not  necessarily  exist.  Inverse  methods  for  transonic  speed  air- 
foil design  have  to  contend  with  this  same  problem.  However,  in  practice  the 
aerodynamicist  has  been  able  to  use  inverse  methods  without  any  undo  hardship. 
Typical  inverse  methods  are  due  to  Tranen  (Reference  81),  Volpe  (Reference  82 1 
and  Carlson  (Reference  83).  The  programs  have  proven  to  be  very  useful  and 
Figures  6l  and  62  show  examples  of  the  application  of  Volpe 's  program  to  high 
performance  and  cruise  section  design,  respectively.  The  cruise  section  is 
an  example  of  a mixed  design,  where  the  top  of  the  original  design  airfoil  was 
retained  and  the  bottom  of  a section  designed  to  have  the  same  upper  surface 
pressure  distribution  as  the  original  and  a lower  surface  pressure  distribu- 
tion with  reduced  pitching  moment  was  employed.  This  example  shows  how  the 
inverse  methods  are  in  practice  quite  versatile. 

The  "nxmierical  optimization"  approach  to  airfoil  design  appeai-s  to  be  a 
very  recent  concept,  unlike  the  inverse  methods,  wliich  were  available  in  the 
forties  for  subsonic  flows  (like  many  of  the  currently  used  aerodynamics 
methods,  inverse  methods  for  detailed  aircraft  work  were  not  routine  engineer- 
ing tools  until  the  widespread  availability  of  computers''.  An  initial  study 
of  numerical  optimization  was  presented  in  1^7^  by  Hicks,  Murman  and  Vanderplaat 
(Reference  84).  The  underlying  idea  in  this  approach  is  to  couple  a modern 
optimization  co>ie  with  an  aerixlyniimic  analysis  cole.  The  airfoil  dc."i,ni 
problem  is  then  cast  as  an  opt  ' ’it  i iT.  p’-.Mor  and  t nr  ent  ri  arpanfis 
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I j the  optimization  method  is  its  ability  to  handle  design  constraints.  These 


constraints  include  both  off  design  performance  requirements  and  design  point 
f geometry  restrictions.  Ballhaus  (Reference  l4)  has  recently  reviewed  the 

I procedure  and  the  report  by  Vanderplaats  and  Hicks  (Reference  85)  provides  the 

I most  detailed  description  of  the  techniques  used  to  formulate  the  design  prob- 

I lem  as  an  optimization  problem. 

Another  approach  to  transonic  airfoil  design  must  be  mentioned  in  any 
review.  Hodograph  methods  have  been  used  to  design  some  very  good  airfoil 
sections.  The  method  has  proven  to  be  very  productive  in  hands  of  the  profi- 
cient operators  at  the  Courant  Institute  (Reference  86).  However,  because  the 
method  seems  to  have  no  application  in  three -dimensions,  it  will  not  be  dis- 
cussed further  in  the  present  study. 

The  optimization  method  of  aerodynamic  design  is  one  of  the  more  promis- 
ing methodologies  under  development;  there  are  some  drawbacks  at  present  which 

t 

need  to  be  examined.  These  drawbacks  are  mainly  related  to  computer  run  times. 
In  optimization  methods  Jargon,  optimization  methods  minimize  some  "objective 
function,"  which  is  a function  of  a set  of  "design  variables"  subject  to  a 
set  of  constraints.  The  "objective  function"  could  be  dra^,  for  example, 
while  the  "design  variables"  are  typically  the  variables  used  to  specify  the 
shape  of  the  airfoil.  The  constraints  could  be  airfoil  thickness,  off-design 
drag  values,  or  virtually  any  other  requirement  that  might  arise  in  practice. 

The  selection  of  the  appropriate  objective  function  and  design  vsuriables 
are  crucial  to  the  success  of  the  optimization  methods.  Although  drag  is  the 
most  direct  measure  of  aerodynamic  performance,  the  transonic  methods  do  not 
presently  predict  drag  within  the  extreme  accuracies  required  for  aerodynamic 
design.  Indeed,  the  computation  of  2-D  transonic  drag  is  only  now  becoming 
a possibility,  with  accurate  engineering  codes  for  general  application  still 
not  available.  Although  incremental  drag  changes  are  quite  useful,  an  optimi- 
zation procediire  that  compares  absolute  values  of  wave  drag,  induced  drag  and 
profile  drag  demands  an  accurate  prediction  of  absolute  drag  values  of  each 
of  these  components.  The  pilot  calculations  performed  by  Hicks  and  co-workers 
have  concentrated  on  minimizing  the  wave  drag  of  transonic  airfoils.  Although 
wave  drag  can  be  computed  over  an  inviscid  foil,  the  shock  which  generates 
wave  drag  also  increases  the  form  draig  due  to  the  thickening  of  the  boundary 
layer  at  the  shock  wave.  Experience  suggests  that  for  typical  modern  transonic 
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airfoils,  the  wave  ilrag  is  small  compareiJ  to  the  form  drag  and  reliable  can- 
putatlonal  methods  which  could  be  used  to  accurately  investigate  the  trade- 
off in  various  drag  components  are  only  now  being  developed.  This  is  an 
example  of  the  present  difficulty  In  specifying  a valid  objective  function. 

As  an  example  of  the  application  of  the  optimization  method  to  transonic 
airfoil  design,  we  exaiiiine  a pressirre  distribution  from  a recent  publication 
by  Hicks  and  Vanderplaats  (Keference  8‘>)  presented  here  as  Figure  03.  Tn  this 
case,  the  wave  drag  is  reduced  from  that  of  the  baseline  foil  by  altering  the 
pressure  distribution  from  one  whicli  accelerated  continuously  into  the  shock 
wave  to  a constant  rooftop  type  pressuie  distrlbuti  on  wtiich  maintains  the  sec- 
tion lift  while  decreasing  the  shock  Mach  number.  This  is  entirely  consistent 
with  current  practice  in  airfoil  design.  It  might  be  more  instructive  to  begin 
with  a section  designed  by  the  hodograili  method  or  Inverse  metiicKl  mid  investi- 
gate improvements  to  these  foils  using  optimization  techniques.  This  would 
appear  to  be  the  appropriate  method  by  wiiicii  the  relative  effectiveness  of 
optimization  methods  could  be  illustrated.  A particularly  interesting  example 
would  be  a case  in  wtiich  the  off  design  characteristics  of  a foil  designed  by 
the  hodograph  raethexi  are  improved  using  optimization  techniques. 

The  "design  variable"  specification  is  perhaps  the  biggest  challenge  in 
the  application  of  optimizatiom  methods.  In  principle,  the  number  of  airfoil 
ordinates  used  to  specify  the  sliape  could  each  be  used  as  design  variable.^, 
however,  if  i>0  upper  surface  and  40  lower  surface  points  (a  typical  iuunl>er  of 
ordinates)  ai’e  used,  then  there  are  100  design  variables,  in  practice,  iu> 
more  than  about  10  Independent  design  variables  can  be  treated  reliably,  iiius, 
the  airfoil  shape  must  be  constructed  from  siiape  liinctions  that  describe  more 
than  a single  ordinate;  i.e.,  coefficients  of  polynomials  used  to  approximate 
airfoil  shapes.  Kx]ierience  led  to  the  reiUization  that  polyniwiials  were  not 
aj)proprlate  "shape  iTinctions",  and  sch»?mes  that  use  linear  civiibi nations  of 
present  supercritical  slmpes  and  local  geometric  perturbations  to  tiiese 
shapes  appear  to  be  the  most  practical  method  to  obtain  u.seful  resiLlts  witii  a 

small  number  of  design  variables.  It  is  important  to  realize  that  the  op-  jj 

timization  method  will  only  identify  the  best  of  a particular  set  of  possible  ^ 

airfoil  shapes  arising  from  the  shape  liinctlons.  If  the  actual  optlmvun  air-  1 

foil  is  not  among  this  set  of  shapes,  the  methoil  cannot  find  this  shape.  !j 

Hence,  the  optimization  methods  require  the  user  to  apply  as  much  (and  pos-  I 
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sibly  more)  insigiit  into  the  problem  then  in  the  use  of  inverse  methods. 

In  summary,  optimization  methods  would  seem  to  deserve  more  attention 
than  they  are  presently  receiving  outside  of  NASA  Ames.  Attention  should 
be  concentrated  on  developing  optimization  algorithmus  that  Incorporate 
aerodynamic  concepts  more  directly,  while  also  providing  information  on  the 
relative  importance  of  the  various  design  variables.  Advanced  airfoil  sec- 
tion designs  will  remain  critically  dependent  on  the  aerodynamic  insiglits 
that  are  not  directly  a part  of  the  optimization  theory;  i.e.,  it  appears 
that  good  airfoils  are  not  described  by  simple  analytical  functions,  lin- 
ally,  the  direct  relationship  between  viscous  effects  and  airfoil  design 
has  been  difficult  to  include  in  a reasonable  optimization  approach, 
d.  Review  of  3-n  Transonic  llesign  Methods 

No  general  three-dimensional  transonic  wing  design  programs  are 
presently  available  in  this  country.  Normally,  transonic  wings  are  designed 
by  procedures  that  employ  2-0  transonic  tools,  panel -method  wing-bexiy  programs 
in  the  design  mode,  and  3-1'*  transonic  wing -body  programs  in  the  tuiulysis  mode. 
The  detailed  design  procedures  are  usually  proprietary,  however,  insight  into 
typical  industry  procedures  can  be  gained  by  studying  Reference  88  find  80. 

One  pilot  3-D  transonic  inverse  method  calculation  has  been  reported  by  Steger 
and  Klineberg  (Reference  90).  Moreover,  Schmidt  and  Hedman  (Reference  oi)  have 
reported  results  obtained  using  an  inverse  method  available  at  Oornler  (Germany^ 
ajid  the  FfA  (Sweden).  In  addition.  Kali  and  Firmin  (Reference  yl  and  Haines 
(Reference  92)  have  reported  the  unpublished  work  by  I^angicy  anil  Forsey  at 
ARA,  Bedford,  England.  Pilot  calculations  applying  optimization  methods  to 
3-D  wing  design  have  been  reported  in  Reference  93* 

The  comments  concerning  inverse  and  optimization  methods  in  2-D  in  the 
previous  section  carry  over  to  the  3-D  design  case.  Because  there  has  been 
very  little  work  reported,  the  review  of  3-D  transonic  design  methotis  is 
necessarily  short.  The  work  on  inverse  methods  by  Schmidt  and  Hedman  (Refer- 
ence 89)  indicates  that  a stable  iterative  procedure  can  be  found  and  a con- 
verged design  obtained.  One  disturbing  aspect  of  the  3-D  inverse  methoil  la 
the  possibility  that  the  solution  may  be  non-unique  near  the  wing  root,  a 
result  which  was  recently  reported  by  Sloof  (Reference  9t().  The  extension  of 
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the  optimization  methods  to  3-D  is  straightforward.  In  fact,  shape  functions 
for  spanwise  variations  would  presiimably  be  less  critical  and  easier  to  con- 
struct than  the  airfoil  shape  functions. 

In  the  next  section,  we  illustrate  the  possible  use  of  the  three- 
dimensional  transonic  methodology  in  a design  environment  by  application  to 
two  model  problems. 
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e.  Application  of  the  S-f  Transonic  IVogram  to  Wing  Design  Problems 
The  feasibility  of  using  the  present  computer  program  in  transonic 
wing  design  as  more  than  a straightforward  analysis  tool  was  Investigated 
througii  two  model  problems.  Hie  first  model  problem^*  was  conducted  without 
using  the  transonic  code,  but  making  use  of  the  NASA  optimization  program 
CONMIN.  The  main  purpose  of  the  excercise  was  to  gain  familiarity  with  the 
use  of  optimization  codes  in  aerodynamic  applications.  The  second  model 
problem  was  undertaken  in  order  to  assess  the  effort  required  to  introduce 
an  autanatic  geometry  alteration  loop  driven  by  the  results  of  a previous 
iteration  into  the  code.  The  stability  of  this  type  of  iterative  procedure 
was  also  of  interest.  An  inverse  method  was  not  investigated  because  the  work 
of  Schmidt  and  Hedman  (Reference  89)  demonstrated  clearly  that  such  a scheme 
could  be  implemented. 

The  first  mcKlel  problem  provided  an  opportunity  to  obtain  a great  deal 
of  experience  using  CONMIN.  The  problem  was  specified  simply  as  follows: 

Dsing  lifting  line  theory  for  the  aci-odynamic  representation  of  the  finite 
wing,  have  CONMIN  detemine  the  twist  distribution  required  to  minimize  the 


induced  dra*.'.  In  this  case  the  exact  solution  can  be  found  to  be 
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for  straight  tapered  wings.  For  an  mitapered  wing,  equation  (oo)  for  ag 
shows  i.hat  the  basic  incidence  variation  along  the  span  is  elliptic.  Ob- 
serving tiie  ’"unctional  form  of  the  exact  solution,  we  note  that  this  parti- 
cular ratio  of  the  root  of  a second  order  polynonial  to  a first  order  poly- 
nomial would  have  beei  an  unlikely  selection  for  the  assumed  variation  of 
spanwise  twist.  To  repeat.,  unless  Equation  (60)  was  contained  as  a subset 
of  the  tdmctiohal  forms  selected  for  the  optim’zation  study  the  true  op- 
timii.iun  twist  distribution  would  not  iiave  been  found.  This  fact  serves  to 
demonstrate  the  importance  of  using  the  insight  gained  from  analytical 
theories  in  order  to  .n;iximize  the  benefits  of  numerical  solutions. 

* The  effort  describeil  iiere  concerning  the  first  model  problem  was 
conducted  witii  internal  vlrumman  iTinding,  but  is  included  in  this  report 
for  completeness  and  perspective. 
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In'leed,  initial  efforts  to  obtain  the  minimum  solution  using  a cubic  poly- 
naiiial  for  the  span  variation  of  twist  were  not  particularly  satisl^yitig. 

The  results  never  approached  the  true  minimum,  and  apparently  tiiere  were 
several  combinations  of  coefficient  values  that  were  equally  close  to  true 
minimum,  such  that  several  substantially  different  answers  for  the  twist 
variation  could  be  obtained,  depending  on  the  initial  guess  supplied  to  the 
program.  These  calculations  ti’pieally  took  on  the  order  of  ten  iterations, 
each  of  wtiich  required  a nunber  of  function  evaluations  in  order  to  obtain 
the  local  gradient  of  the  objective  function.  In  aerodynajiiic  terms  this  memis 
that  there  were  ten  main  executions  of  the  aerodynamic  program,  and  a number 
of  "small"  executions  wiiich  were  required  to  be  run  long  enougli  to  provide  the 
local  gradient  of  the  solution  with  respect  to  each  design  variable.  It  is 
clear  that  this  can  quickly  lead  to  an  immense  tunoimt  of  computational  effort. 

Finally,  the  optimication  scheme  was  run  with  the  design  variables  con- 
sisting of  a coefficient  to  equation  (66)  and  tiie  coefficient  of  ;ui  additional 
term  added  to  equation  (66).  Figure  (J|  shows  the  patli  tiirough  design  space 
for  this  two  parameter  optimisation  run.  Note  that  the  minimum  occurs  wlien 
^ -=  0,  and  ■«  1 (P.?*l  exactly  because  a lift  cvu’ve  slope  slightly  dif- 
ferent then  2 n was  employed'',  iiie  run  temiinated  after  elgjit  iterations,  with 
the  numerical  solution  predicting  that  the  optimim  had  been  achieved.  The  re- 
sult demonstrated  tliat  the  program  could  in  fact  select  the  ti-ue  optimim  if  it 
wjis  embedded  in  tiie  design  variable  space.  This  effort  demonstrated  both  tlic 
difficulties  and  possibilities  associated  with  the  use  of  optimisation  metliods. 

'Ihe  second  model  problem  was  considerably  different  in  concept.  For 
ttiis  problem  the  question  posed  was  simply:  For  a given  planform  mid  spfuiload, 
determine  the  twist  required  to  produce  the  spanload.  Initially  lifting  line 
thc-'ry  was  employed  to  verify  that  the  basic  iteration  sciieme  adopted  would 
converge  for  a simple  aerodynamic  model  before  attempting  to  incorporate  tlie 
iteration  into  tlie  main  program.  The  twist  was  determined  by  adjusting  the 
seecion  incidence  at  the  finite  set  of  span  stations  at  vrtiich  the  conputation 
provided  results,  without  making  any  assumption  concerning  the  functional  re- 
lationship between  the  incidence  at  adjacent  span  stations.  The  basic  itera- 
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where  ^1  denotes  the  particular  span  station,  P denotes  the  desikai  condition  luid 
K indicates  the  iteration  tuunber,  is  approximated  by 


por  the  lil*tini^  line  simulati^ai  this  iteration  procedure  converfied  to  tlie 
exact  solution  »?iven  by  Kquatlon  (i>6)  in  about  four  or  five  iterations.  This 
result  was  obtained  wltliout  difficulty  even  tliougii  Uie  approximation 

given  in  Kquatliin  is  very  pinir  for  munerlcal  computation  due  to  tlie  pro- 

gressively smaller  differences  between  tlie  values  as  the  iteration  converges. 

Kquation  (o7^  is  enuivaJent  to  a more  general  form: 


1^  i 


[» 


Wliere  [AJ  is  sui  approximaticui  to  the  actual  influence  coefficient-  mat.rix 
wliich  relates  and  a: 

fC^]  --  [Al  fol. 


In  the  present  method  [ 'K ] has  been  given  by  the  extremely  simplified  relation- 
ship in  Equation  (till  for  the  diagonal  terms,  with  the  off-diagonal  tei-ms  as- 
sumed to  be  cero.  This  result  shows  that  [A]  can  be  crudely  approximated  if  :ui 
iteration  to  determine  the  final  result  is  allowed.  Naturally,  as  the  appi’oxi- 
mution  to  [a1  improves,  the  number  of  iterations  required  is  reduced. 

Moiiifi  cat  ions  to  the  basic  invlscid  program  to  include  tliis  type  of  iter- 
ation scheme  were  ineoriwrated  without  difficulty.  It  was  found  that  a rela- 
tively fine  grid  was  required  in  order  to  obtain  tiie  straiglit  wijig  result  com- 
puted previously  using  the  lliXing  line  aeri.xiyn!unic  mode.  Kefinements  to  the 
iteration  included  tlie  use  of  underrelaxation  of  the  twist  increment  ;uid  tlie 
use  of  the  Initial  value  for  all  iterations.  These  refinements  led  to  a smooth- 
ly converging  solution  that  took  about  fifty  percent  l».mger  Unut  the  basic  solu- 
tion. Ilie  metluxl  wjis  then  applied  to  tv  swept  luitapered  wing.  The  refined 
procedure  led  to  a solution  with  the  results  obtained  shown  in  Elgure  I'p,  wliich 
also  contains  the  str.aiglit  wing  results.  In  this  case,  attempts  to  caiquite  Uie 
result  wliileCj^  changed  from  iteration  to  iteration  led  to  a diverging  result 
at  the  point  wliere  no  rlilft  in  tingle  was  required  (about  4V){>  semispan tmd  shows 
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Uint  in  an  aotnal  piwluol  ioi\  pin't-nnuii  :ui  impi-ovtni  approx iiiiat. itin  bo  [A  1 r.lu'aJii 
bo  ino.lnioJ.  iiowovoin  this  improvod  [A'  oould  \uuloubtodly  bo  oonstniotod  witli- 
v'Ut  dibt'iouity  no  that  a doaii’in  optivin  v't'  tho  tyiu'  dosoribod  abovo  oou.Li  bo 
iiK'lvidod  in  tho  baaLi-  ptnv'inuit  without  -ii  t't'ioult;,  . 

in  this  soot  ion  wo  havo  vioiv.oitst  rat  od  tlio  varioty  of  possibiii  l ion  that  aviso 
wluni  inoorporation  of  dosif^n  options  is  su^^t-rostovi . I'no  of  ttio  options  would 
pi'ovido  i’lniodiato  bonofits  to  ttio  dosij^nor,  wtiilo  itio  optiiiiiiiation  approaoh,  al - 
tlK'U^'lt  'ptito  pivvm isiti;',  woulvl  tnko  tnuoii  tiun’o  study'  in  oivior  to  produot'  a rolia- 
I'lo  dosi^n  tool. 

A braotioal  Approaoh  to  ;u!  intonnoii iat o i'osign  i>\xlo 

I'ho  i-osults  of  tlio  offorts  roportod  thus  far  provido  a basis  for  oon- 
sidoriiv  tiow  host  to  appvi\\ol\  tlio  quostion  ot'  actually  imploniontint.":  a dosip.n 
option  in  tlio  prosont  codo.  i'ho  procodiiut  discussion  ajid  oxiuni'ios  havo  showii 
that  a varioty  of  dosi;:;n  inothods  could  bo  adopted,  i’ho  f^oal  of  tho  approach 
dosoribod  in  t-his  soctiivi  is  to  provide  a rationai  foundation  for  a cixlo  wtiioh 
would  provido  till'  dosiftnor  with  a roliabio  tool  ^juickly.  while  also  doinonst ra- 
ting: basic  tochniduos  tluvt  could  bo  adopted  in  more  sophisticated  triuisonic 
ra'tluxiolvvti^'-’  as  they'  are  dovolv’pod.  .-Iny  resulting:  codo  must,  bo  easy  to  use. 

I'iio  of  tho  keys  to  a succosstdil  result  is  to  incorporate  a number  of  features 
one  at  a time  in  order  to  reduce  tho  risk.  Kach  i.vio  of  those  items  couLl  bo 
tested,  sop.aratoly  in  order  to  assess  its  practicality  before  including  it  in 
the  actual  design  optiiiii  prognuu.  i'ho  approach  described  here  moots  t.hese 
requirements  by  proposing  a "bleiuioii"  metluxioliigy  wliich  c>.imbines  inverse  tuui 
I'pt  imittalion  metluxls  in  a iiuuiner  that  eliminates  the  weaknesses  of  eaeli  of 
tlie  metluxls  wlien  used  seyviratoly. 

I’ho  tr;uiso;iic  wing  vlesign  procedure  wv'  envision  proceeds  as  follows; 

tl'  A "target"  pressure  di  st  ribut  iiin  sliould  be  specified  in  the 
manner  of  present  wing  design  work.  I’liis  requires  Svxiie  oxjH'rience  on  the  p.art 
of  the  fierotyauutiicist . but  this  step  is  the  vxie  in  which  he  c!Ui  contribute  his 
insights  tuid  knc'wUxlge.  Headers  not  fajiiiliar  with  the  tyq^ical  isobar  ^iesign 
pix-cedurer.  employed  In  industry  should  consult  the  papers  by  iiaines  ^Kefei-ences 
OS  and  o,A,  it  will  be  svwe  time  before  the  cimiputer  opt  i mi  •.'.at  ion  procedures 
ciui  ciVTipetent ly  eliminate  tliis  step. 


(2''  The  related  baseline  P-h  airfoil  section  (or  sections'  aliould  he 
ileslf^ned  with  inverse  metluMs  in  the  manner  presently  used.  The  mmiericaJ 
resolution  availahle  in  f'-p  piaY.rtms  exceeds  juiy  3-h  piavriun,  so  that  baseline 
sections  should  always  be  desif;ned  usinc:  the  f'-P  methods.  The  l-P  desii^n 
et'Torts  will  tiien  considei’  perturbations  to  tJiese  baseline  shapes. 

(3^  Ttie  basic  wine';  twist  required  to  obtain  the  vlesii-ed  spfuiiwtd 
usinei  tl\e  P-P  sectlvins  fv.'iu\d  in  stei'  should  be  deteiaiiined  by  ;ui  itei-ative 
procedure  th-it  w>.'uld  be  ;ut  extension  of'  tfie  ituxlel  problem  pj-ocedure  described 
above.  At  this  point,  tin  d.esLj'ner  shoulii  have  a f^ood  first  winy  desiyn.  wtiicli 
Cvuild  be  coji  si  dent'd  tiie  baseline  winy  Tor  ilu'Uie!'  inqiivvement  s. 

(4'  In  reyicMis  oji  the  winy  wtiere  the  taryet  pressurs'  di stri but  i t'n  has 
not  beei\  aeliieved,  the  neroi.vauimie i st  should  speeiiy  tiie  desircvi  pressure 
distribution  tuid  ;ui  inverse  method  should  bo  employed  in  order  to  determine 
ttie  required  yeomt'try  modifi eaticai  necessary  to  acb.ieve  tb.e  tar.yet  pi-essure 
di  st  ’•  ii'utiivf . 


I'he  key  step  in  the  pri'oovifire  is  to  ox-uuine  t'.ie  yeomotry  modiPiea- 
tlonr  Tt'Utni  in  stt'p  'i . ilu'  result  s found  in  that  stop  probably  violat  ed  so.ne 
basic  ctaistraints  associated  with  the  desiyn,  such  as  section  thiekut'ss  t'r 
closure,  howt'ver.  the  inverse  methoi  rJuniid  !>ave  provitievi  the  desirt-ti  shape 
tlinction  whieli  e:m  tiien  be  used  i)\  an  I'pt  imi ratiei-i  pri'cedure,  wliieh  can  then 
be  used  to  mittimi’ce  tlie  .lifferi'iiee  I't'tween  ttie  taryet  pressvire  distrii'vdion 
ajtd  the  best  possilvU'  pressure  distribution  wnicli  e;ui  i'c  ol'taincvi  while'  also 
meet  in, y ci'iistraints  such  as  yi'o.a'trie  limitations  and  off-dosiyn  per  foriuan.ee. 

Tlu'  *'  step  proei'diirt'  I'Ut  lined  ab,'ve  requi  n's  tlu'  nu' ii  I'ieat  ion  I'f  tin' 
small  vi i st \u'banee  iheoey  proyriuii  to  inelndt’  a 'leval'  inverse  methoi  wnoi'e'in 
yt'.'metry  -uid  pi'i'sr.nrt's  e:ui  both  he'  speoil’ie.i  I'ver  dit'feri'nl  portions  oi'  the 
pbuil'orm.  I'he  other  reqn i re'mi'nt  is  aj^  opt  Imi  . ation  class  v'f  alyor  i thm  whieii 
e*!Ui  tutumie.e  the  »lifferi'nee  bi'twi'i',  two  pressure  distribntiens,  subje'ct  te  -i 
varie'ty  ef  eenst  ra  int.s . Peth  ot'  tiu'so  capabilities  >'!Ui  be'  ineerporat  <'d  with- 


out difficulty. 

The  uevel  idea  eei\t.ai!ud  in  this  appreaeli  is  tlu'  use  ot'  the  it^verse 
metiu'et  te'  suyyi'St  Ih*'  r.hape  t'lriet  ious  whieli  e:ui  bo  vised  as  desiyn  variable's 
in  tJie'  e'ptimlrat  ii'ti  prevovinre',  i'his  prv''ee'eiure  remewos  the  main  d.e'f i e ie-ney 
in  t.'u'  I'lptimi rat  K'n  me'tluxls  by  int  re^eUiciny  a ’’iritural''  tihape  tlinetie'n  which 


sluxilvi  ullow  the  optiitii;*.ati(.in  procedure  to  eonverf-te  quickly  to  the  final 
result.  lliis  fast  convergence  is  IXirtlier  sided  by  the  fact  that  the  initial 
starting  point  for  the  optimi::ation  progrjuii  sliould  be  "close"  to  tlie  final 
result,  because  the  inverse  method  is  not  required  to  produce  a desigii  meet- 
ing the  geametrical  constraints,  its  ni.ajor  drawback  is  removed,  by  blending, 
the  two  meUKxlologies  in  the  proposed  nrumer,  lui  outstiuuling  traiuuvnc  wing 
design  tool  coxild  be  developed. 

g.  «.\'i;iolusions  v'ii  resign  Feasibility 

It  is  feasible  to  construct  ;ui  outstajiding  triuisonic  wing  desig.n 
tool  based  on  the  present  trtuisonic  ;uialysis  metlu>iol(.\:y.  This  coile  would  be 
a blenvi  of  inverse  ;uid  optimisation  mettiods  ttiat  would  remove  t.he  restrictions 
ii\herent  in  eitlier  approacii  wlien  considered  separately.  Hie  resulting  program 
wo.ild  make  only  sliglit  additiv'inal  demtuids  on  ci.imputing  system  resoui-ces. 

having  outlined  in  the  last  section  Uie  approac);  t.hat  'would  provide  this 
tool,  we  rooivimend  tluit  tne  developnent  of  a blended  inverse  optimisation 
method  be  undertaken.  I'his  effort  would  bo  based  on  metliods  wnioh  have  alrea.ly  -a 
been  demonstrated  separately,  insuring  tlie  suecess  of  tlie  new  effort.  Farliou- 
la!'  atteiitii,'!!  sliould  be  pJaet'd  ivi  t.he  Aiudiuiiental  aspects  of  tiie  melhodolog.,v 
in  order  to  prv'ivide  a solivl  t'oiuulation  t'or  oontiirued  development  along  tho 
ines  outlined,  i.lie  of  the  more  signifieiuit  features  of  tlie  blended  approach 
is  tliat  the  optimisation  method  is  basically  independent  of  tiie  aero.i.mam ic 
motliods  .'ui.i  will  be  partioularly  insensitive  to  tlie  rei'inoment  of  tlie  aero- 
d>aiam i e met lio.iv' Log.,v . 
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Figure  65.  Wing  Design  Example 


SECTION  III 


CONCLUSIONS 


Major  advances  are  beinc  made  in  the  computation  of  transonic  flows 
over  wing-body  combinations.  In  the  present  work  on  inviscid  three- 
dimensional  small  disturbance  theory  program  has  been  coupled  with  a 
boundary  layer  program  and  the  resulting  code  re-designed  to  be  extrem- 
ely easy  to  use. 

In  the  new  progran  the  user  has  the  ability  to  examine  a number  of 
effects; 

o Viscous  effects  via  a fully  coupled  strip  boundary  layer  of  infinite 
swept  wing  type. 

o Local  strong  shock  interation  at  sliock  waves 

o Shallow  separations  at  the  trailing  edge  and  in  the  cove  region 
of  supercritical  airfoils 

o Reynolds  number  and  boundary  layer  transition  location 
o Infinite  yawed  wing  results  including  2-D  cases 
o Fully  conservative  (FCR)  and  non-conservative  (NCR)  inviscid  so- 
lution procedures 

o Modified  small  disturbance  theory  (MSI')  or  classical  small  distur- 
bance theory  governing  equations  (CSD) 
o Riegels'  Rule  effects 

o Body  effects  on  wing  pressures,  including  simple  area  ruling 
effects  via  body  lines  input,  or  detailed  fuselage  effects  via 
body  slope  input 

o Ptilly  three-dimensional  boundary  layer  predictions  via  automatic 
data  set-up  for  the  new  boundary  layer  program  written  during  the 
present  effort. 

In  order  to  provide  a user-oriented  program  the  following  special  pro- 
visions were  included: 

o A complete  visualization  of  the  program  input  and  results  are  avail- 
able via  the  automated  grapliics  package. 


IQY 


o 


An  internal  check  is  made  of  input  data  sets  for  consistency,  and 
the  program  stops  when  the  solution  is  either  converged  or  encoun- 
tering difficulties  in  order  to  conserve  computing  time. 

The  resulting  computer  code  was  applied  to  numerous  configurations , 
wherein  the  identical  code  was  used  for  each  case  without  any  special 
"tuning" . 

The  conclusions  resulting  from  these  studies  can  be  summarized  as 
follows : 

o For  relatively  simple  conceptual  type  aircraft  configurations,  the 
program  provides  excellent  results  for  a wide  range  of  Mach  numbers, 
sweep  angles  and  taper  ratios. 

o The  agreement  with  the  advanced  technology  transonic  transport  re- 
sults provide  confidence  in  the  prograjn  for  transport  type  aircraft 
o Advanced  technology  fighter  configurations  present  special  problems 
wliich  still  require  some  resolution,  such  as  the  detailed  treatment 
of  wing -glove-body  blending.  Nevertheless,  encouraging  results  have 
been  obtained  in  several  instances  although  in  many  cases  a nvmiber  of  com- 
plicating factors  prevent  a final  determination  of  the  source  of 
the  disagreement. 

Perhaps  the  single  most  important  result  of  this  effort  has  been 
tlie  development  of  a code  which  allows  the  routine  use  of  the  new  3-P 
transonic  methods  outside  of  a few  government  and  industrial  research 
labs.  The  present  code  and  its  documentation  provides  a baseline  upon 
which  an  aerodynamicist  has  extreme  flexibility  to  investigate  both 
aerodynamic  configurations  and  program  modifications.  The  resulting 
program  should  have  major  impact  on  the  design  of  advanced  aircraft. 


SECTION  IV 


RECOMMENl^ATIONS  FOR  FUTURE  WORK 

IXie  to  the  extremely  broad  range  of  aerodynamic  effects  investi- 
gated in  the  present  work , the  recommendation  for  future  work  tends 
to  be  all  encompassing.  The  various  requirements  can  be  divided  be- 
tween basic  research,  general  requirements  for  aerodynamic  code  develop- 
ment and  specific  extensions  and  investigations  related  to  the  present 
program . 

Basic  research  programs  which  would  have  a direct  application  to 
the  computation  of  transonic  flowfields  over  wing-body  configurations 
are: 

o A baseline  experimental  three-dimensional  boundary  layer  study 
for  a typical  supercritical  wing  body  combination.  This  study 
is  required  in  order  to  validate  tlie  fiilly  three-dimensional  boiui- 
dary  layer  prediction  methods  and  to  assess  the  range  of  applicabil- 
ity of  the  more  approximate  methods  such  as  tlie  infinite  swept 
wing  strip  boundary  layer  approach. 

o Basic  theoretical  work  in  the  area  of  separated  flovc  is  required, 
witli  the  specific  goal  of  producing  rational  engineering:  methods 
tiiat  could  be  incorporated  into  codes  designed  along  the  lines  of 
the  present  program.  This  work  should  include  tlie  calculation  of 
the  three-dimensional  viscous  wake  downstivam  of  finite  wings, 
o Basic  nimterical  research  into  tlie  acceleration  of  the  inviscid 

solution  process  is  required  in  order  to  reduce  the  computationaJ 
cost  of  tlie  calculation. 

dieneral  requirements  for  I'uture  computational  aerodynamic  code 
development  should  include; 

o A uniform  input  fomat  for  body  and  wing  geometry  should  be 
adopted. 

o Standard  output  format  and  nomenclature  should  be  adopted. 

Specific  extensions  of  the  present  code  that  would  prove  in- 
valuable in  the  future  use  of  the  metluxi  include: 
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o A general  re -examination  of  the  body-glove^wlng  .lunoturc,  and 

the  mapping  bending  near  the  body  is  needed  in  order  to  improve 
the  quality  of  Uie  local  solution, 
o djeivnetrlc  extensions  to  the  configurations  that  caji  be  hmidled 
should  be  considered.  Many  advanced  technology  configurations 
Will  require  the  explicit  treatment  of  canards  find  winglets.  In 
addition,  nacelle  inlet  effects  are  usually  required  in  order  to 
perform  detailed  transonic  wing  designs.  Ibese  features  will  have 
to  be  incorporated  if  the  program  is  to  Ailly  supplant  the  present 
day  panel  methods  as  a design  tool  at  transonic  ai^eeds. 
o Steady  state  aeroelastic  effects  must  be  included  wlien  computing 
tile  transonic  Qi'w  over  mo.iern  aircraft.  The  wing  deflection 
could  be  incorporated  at  the  same  time  the  boundary  layer  correc- 
tions are  made  wi  tii  a minimu.n  of  additional  computati onal  effort. 

This  would  greatly  enliance  tlie  general  utility  of  ttie  code, 
o The  restriction  to  Mach  mmibers  less  than  one  is  not.  a basic  limit- 
ation of  any  of  the  methods,  and  tlie  code  could  be  extended  to  tlie 
upper  transonic  regime  (M=.1.2‘5)  witliout  difficulty.  This  extension 
woulvl  complete  the  bridging  of  the  linear  subsonic  :uid  supersonic 
regimes  so  that  there  would  be  no  "gaps"  in  the  range  of  Miicli  numbers 
over  wliich  nxitine  aerodynamic  predictions  can  bo  made, 
o Substfuitial  reductions  in  the  real  tuid  computational  ti  e requireii 
for  wing  design  Wv->rk  could  be  achieveil  by  ttic  incorporation  of  sune 
desigji  options  into  tlie  program.  Tlie  options  sliovild  include  local 
Inverse  solution  capability,  and  a simple  application  of  tlie  geo- 
metry alterati«.in  iteration  procedures  luialogous  to  those  in  op- 
timlzatiivi  techniques. 

Althougii  the  prognuii  outlined  above  is  si-xnewliat  imibitious  it  does  in 
fact  reflect  tiic  broad  scope  of  the  effoid  luidertiiken  in  tlie  present  j'rogriuii. 
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